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PREFACE 

Initially,  my  thesis  project  was  to  set  up  Gulf  Gener¬ 
al  Atomic’s  multigroup  neutron  cross  section  code,  GGC-4,  on 
an  IBM  7094  digital  computer.  I  was  to  understand  the  theory 
used  in  the  code,  be  able  to  operate  the  code  well  enough  to 
document  its  operation  for  others  (produce  a  local  user’s 
manual)  and  compile  a  source  book  of  neutron  group  cross 
sections  for  weapons  problems. 

I  obtained  the  GGC-4  code  and  the  GGC-4  cross  section 
library  (data  for  45  nuclides)  from  the  Air  Force  Weapons 
Laboratory  where  a  CDC  6600  digital  computer  was  used  to 
copy  the  code  and  the  cross  section  library  onto  magnetic 
tapes.  I  was  unable  to  program  the  GGC-4  code  on  the  IBM 
7094  due  to  what  appeared  to  be  insufficient  computer  memory, 
but  with  the  possibility  of  it  being  due  to  a” compatibility" 
problem  with  the  CDC  6600  output  tape  being  used  on  the 
IBM  7094.  I  was  eventually  successful  in  obtaining  data 
for  the  45  nuclides  from  the  GGC-4  data  tapes. 

My  thesis  advisor.  Dr.  C.  J.  Bridgman,  and  I  decided 
to  write  our  own  code  for  the  IBM  7094»  This  code  calculates 
the  zero  moment  of  the  neutron  flux  which  is  used  to  flux 
weight  the  GGC-4  cross  section  data.  This  code  was  the  bulk 
of  my  thesis  work.  In  the  course  of  this  work,  I  was  able  to 

obtain  two  new  ENDF/B  neutron  cross  section  data  tapes  from 
Oak  Ridge  National  Laboratory.  I  was  able  to  use  these 
ENDF/B  data  tapes  with  minor  modifications  to  my  original 
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moments  code. 

A  few  words  should  be  said  about  the  title  of  my  code, 
New  Barnyard.  The  Resident  School  of  Engineering  at  AFIT 
has  a  cross  section  code  called  Old  Barnyard  that  calculates 
group  cross  sections  for  the  energy  range  10  MeV  to  thermal. 

Old  Barnyard's  cross  section  library  was  updated  by  using 
the  cross  section  results  from  my  code  as  the  "new"  input 
library  to  Old  Barnyard.  Thus,  the  idea  occurred  to  me  of 

calling  my  code  New  Barnyard. 

I  am  indebted  to  several  people  in  connection  with  this 
thesis:  to  my  advisor,  Dr.  C.  J.  Bridgman,  for  guidance  and 

counseling;  to  Lieutenant  Robert  Barry  of  the  Air  Force 
Weapons  Laboratory,  Kirtland,  AFB,  for  help  with  the  GGC-4 
library  tapes;  to  my  classmates,  Captains  Jim  Fisk  and 
Robert  Winchester,  and  Lieutenants  Gary  Knutson  and 
Fred  Damm,  for  many  interesting  discussions  on  transport 
theory,  and  to  my  typist,  Mrs.  Bobbie  Thompson.  Finally, 

I  express  special  thanks  to  my  wife,  Cely,  for  her  encour¬ 
agement  and  good  humor  throughout  my  thesis  study. 
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ABSTRACT 

New  Barnyard,  a  moments  code,  was  written  to  calculate 
multigroup  neutron  macroscopic  cross  sections.  Total,  ab¬ 
sorption,  fission,  and  scattering  transfer  (group  to  group) 
cross  sections  can  be  calculated.  The  transport  cross  sec¬ 
tion,  diffusion  coefficient,  and  the  average  cosine  of  the 
scattering  angle  can  also  be  calculated  for  each  group.  The 
energy  range  of  these  cross  sections  extends  from  14.91B  MeV 
to  .4139  ev.  Two  versions  of  the  moments  code  were  written 
so  that  two  different  data  sources  could  be  used.  With  one 
version  PO,  PI,  P2,  and  P3  elastic  scattering  transfer  cross 
sections  can  be  calculated  for  22  broad  groups,  and  with  the 
other  version  PO  through  P$  elastic  scattering  transfer  cross 
sections  can  be  calculated  for  20  broad  groups. 

New  Barnyard  was  written  in  the  Fortran  IV  language 
for  use  on  an  IBM  7094  digital  computer.  The  code  calcu¬ 
lates  the  zero  moment  of  the  neutron  flux  which  is  then  used 
to  flux  weight  basic  neutron  cross  section  data  over  the 
energy  limits  of  each  group.  Unlike  some  other  moments 
codes,  no  first  and  second  moments  of  neutron  flux  are 
calculated,  and  no  B^,  Pl»  ag e,  or  resonance  calculations 
are  performed.  These  calculations  were  excluded  to  achieve 
simplicity  and  speed  of  calculation.  New  Barnyard  uses  the 
GGC-4  cross  section  library  as  well  as  the  new  ENDF/B  data 
from  ORNL  (Oak  Ridge  National  Laboratory). 

Excellent  agreement  was  found  when  the  results  of  New 
Barnyard  were  compared  with  results  from  other  cross 
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ABSTRACT  (Contdf) 

section  codes.  Results  from  the  two  versions  of  New  Barn¬ 
yard  showed  relative  differences  generally  much  less  than 
10^  for  comparisons  between  absorption  and  between  total 
cross  sections.  Larger  relative  differences  were  found  in 
the  comparisons  of  the  PN  scattering  transfer  cross  sections 
due  mainly  to  the  fact  that  the  number  of  terms  used  in  the 
Legendre  polynomial  expansion  of  the  differential  scattering 
cross  sections  for  each  set  of  data  on  the  two  data  tape 
sources  were  different. 
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NEW  BARNYARD:  A  MULTIGROUP 
NEUTRON  CROSS  SECTION  CODE 
I.  Introduction 

A  need  exists  for  multigroup  neutron  cross  sections 
for  weapons  problems.  Although  many  sets  of  group  cross 
sections  exist,  most  are  based  on  reactor  spectra  and  con¬ 
sequently  are  of  limited  use  from  about  2  to  15  MeV,  the 
principle  region  of  interest  for  weapon  physics. 

Purpose  and  Method 

The  purpose  of  this  study  is  to  produce  a  code 
that  will  calculate  group  neutron  cross  sections  in  this 
higher  neutron  energy  range.  The  computer  code  described 
in  this  study  calculates  group  neutron  cross  sections  and 
related  constants  for  the  energy  range  .4139  ev  to  14.91$ 
MeV.  It  is  written  in  Fortran  IV  and  has  been  executed  on 
an  IBM  7094  computer.  The  cross  sections  are  determined  by 
first  calculating  the  energy  dependent  flux  in  an  infinite 
homogeneous  mixture  of  the  isotopes  or  compounds  specified. 
This  flux  is  then  used  to  flux  weight  cross  section  data. 

The  input  cross  section  data  is  obtained  from  99  group 
libraries  of  fast  cross  sections.  The  energy  dependent  flux 
is  calculated  by  the  method  of  moments. 

Sequence  of  Development 

The  basic  theory  for  the  numerical  equations  used  in 
New  Barnyard  is  given  in  Chapter  II.  The  cross  section 
data  tapes,  the  numerical  equations,  and  the  special  output 
features  of  New  Barnyard  are  presented  in  Chapter  III.  The 
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operating  instructions  for  using  New  Barnyard  on  an  IBM  7094 
are  given  in  Chapter  IV.  In  Chapter  V,  results  from  other 
cross  section  codes  are  given  and  then  the  conclusions  reached 
on  the  validity  and  accuracy  of  the  code  are  stated. 
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II.  Theory 

In  this  chapter  a  brief  review  of  the  theory  of  flux 
weighting  is  given  followed  by  a  theoretical  discussion  of 
the  calculation  of  the  zero  moment  of  the  neutron  flux  which 
is  used  to  flux  weight  cross  sections  in  New  Barnyard. 

Flux  Weighting 

3 

The  reaction  rate  R  (events/cm  sec)  in  the  pre¬ 
sence  of  polyenergetic  neutrons  is 


R  = 


(1) 


where  E  is  the  energy, 

X  2 

(J*  (E)  is  the  energy  dependent  flux  in  Qn  /cm  -sec, 

E(E)  is  the  energy  dependent  macroscopic  cross  section 


-1 

in  cm 


The  total  flux  is,  by  definition 

p  =JFl(E)dE  (2) 

In  terms  of  this  total  flux,  (J>,  the  reaction  rate  can  be 
written 

R  =E  <j>  O) 

where  £  is  some  average  macroscopic  cross  section  over  the 
energy  range  of  interest.  It  follows,  by  equating  the 
reaction  rates,  i.e.,  the  right  hand  sides  of  equations  (l) 
and  (3)  and  substituting  from  (2),  that  this  average  cross 
section  must  be  given  by 
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(E)  dE 

r(XJ - 

j>  (E)  dE 


(A) 


-'o 

or,  in  other  words,  the  average  cross  section  is  a  flux 
weighted  average  of  the  energy  dependent  cross  section. 

Equation  (4)  can  be  generalized  to  any  range  (limits 
on  the  integral),  say  from  E^  to  E2,  in  order  to  produce  an 
average  cross  section  applicable  to  that  energy  range.  Such 
averages  are  called  energy  group  cross  sections,  and  may  be 
expressed  as 


(5) 


where  £  is  the  group  cross  section  for  the  n-th  energy  grotp, 

E”  is  the  lower  energy  boundary  of  the  n-th  group, 

E+  is  the  upper  energy  boundary  of  the  n-th  group, 
n 

From  inspection  of  equation  (5)  it  is  seen  that  the  varia¬ 
tion  of  neutron  flux  with  energy  must  be  known  in  order  to 
determine  group  cross  sections. 

Zero  Moment  of  the  Neutron  Flux 

The  flux  calculations  performed  in  this  code  are 
in  solution  to  the  energy  dependent  Boltzmann  transport 
equation  for  an  above-thermal  energy  region  (14.918  MeV  to 
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.4139  ev) .  The  energy  dependent  flux  for  this  energy  range 
is  calculated  by  the  method  of  moments.  A  brief  review  of 
the  zero  moment  of  the  neutron  flux  will  be  discussed  here. 

A  more  detailed  analysis  of  the  method  of  moments  can  be 
found  in  AFIT  Technical  Report  67-17  (Ref  1). 

The  steady  state  Boltzmann  equation  written  for  an  in¬ 
finite,  non-multiplying  homogeneous  medium  which  scatters 
and  absorbs  neutrons  is 

-  S(X,E ,fi) 

r  r°° 

+  J&’  Jd E  Es(e’^e,M0)^(x,e ',M') 

O 

(6) 

Where  X  is  the  spatial  positions  of  the  neutrons 
E  is  the  neutron's  energy, 

jUQ is  the  cosine  of  the  scalar  angle,  cosine  Qg, 

through  which  a  neutron  is  scattered, 

U-OT  fl'  is  the  cosine  of  the  scalar  angle,  cosine  9,(0*) 

between  the  neutron's  direction  and  the  X  axis, 

is  the  neutron  flux  in 
i  ,  2 

0n  /cm  -sec-steradian-unit  energy 

(E)  is  the  energy  dependent  macroscopic  total 
cross  section  in  cm"\ 

Y  (E' — *E >fX  )  is  the  energy  dependent  macroscopic 
s  o 

differential  scattering  transfer  cross  section 
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for  a  neutron  with  an  initial  energy  E'  and  an 
initial  direction M'  that  scatters  into  a  unit 
energy  interval  about  E  and  within  a  unit  solid 
angle  about  /X  in  cm"*  /steradian. 

In  the  above- thermal  energy  region  the  target  nuclei 
motion  may  be  neglected  with  respect  to  the  neutron  ener¬ 
gies.  The  moments  method  assumes  that  an  isotropic  source 
consisting  of  a  plane  of  infinite  area  is  located  at  the 
coordinate  position  x  =  0,  i.e., 

Me) 

S  (X,E ,jl)  -  <?(X)  (7) 

Where  (f (X)  is  the  Dirac  data  function  at  X  =  0.  The  flux 
in  equation  (6)  is  expanded  in  terms  of  Legendre  polynom¬ 
ial  as 

oo 

(JKX.E.M)  -  Y,  2^—  lL(X’E)PmW  (8) 

m  =  o  4  /( 

Similarily,  the  differential  scattering  transfer  cross  sec¬ 
tion  is  expressed  as 

CO 

2 n  +i 

Es(EWE,M0)  -  L  — =-  Esn(E'-E)Pn0io)  (9) 

n=o  4/( 

The  coefficients,  ^]  (E * — >E)  ,  in  the  polynomial  expansion 

sn 

in  equation  (9)  are  referred  to  as  PN  scattering  transfer 
cross  sections  where  PN  is  given  by 
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PN  (E ' — >E)  =  (2n  +1)  V  (E »E) 

sn 


(10) 


where  N  =  n. 

The  substitution  of  the  source  given  by  equation  (7) 

and, of  the  expansions  (8)  and  (9)  into  equation  (6)  yields 
00  ' 


2mfl  a(t»(X,E).  2m+lc_,  S(E) 

kU  TT~UPm (M)+4TT  £c<E>i(x>E)pm(^)=  + 


DO  co 


fdQ'J  d^T 


n^o*- 


2n+l 


47 T 


£  (E'->E)P  (U) 

^sn  no 


2m+l 


_  47T 


(|)  (x,e')p  (/XO 

Tm  m 


1(11) 


Equation  (11)  is  simplified  by  using  some  of  the  special 
properties  of  Legendre  polynomials  (Ref  2:  115)  to  obtain 

(X,E) 


'  m 


m=o 


|3X 


(nri-l)P  UU)  +  mjp  (AX) 
m-H  m-| 


(2mfl)Zt(E)j)m(X,E)Pm(/X)  =  S(E)«f(X)  + 

"  f° 

(2m+l)P  JfM)  JdE'  £  Sm  (E '  ^E)  (j>m  (X ,  E '  ) 


(12) 


Equation  (12)  ,  a  single  equation  with  an  infinite  number  of 
terms,  is  transformed  into  an  infinite  number  of  coupled 
equations,  each  with  a  finite  number  of  terms,  by  operating 

on  equation  (12) ,  term  by  term,  with 
■±1 


J <eq.(12))<fM  1-0 

The  results  for  JL=  0 


1  2 

,  1.  ,  ,  •  •  • 


are: 
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9(1)  (X,E)  r 

- 1 - +  £  ♦  (X,E)  =  V  (Ef — >E)(|)  (X ,  E 1 )  dE '  +S  (E)  (f  (X) 

9  X  to  0J  so  o  o 


(13) 


Equation  (13)  is  the  basis  for  the  calculation  of  the  zero 
moment  of  the  neutron  flux. 

The  basic  energy  dependence  of  the  flux  is  expressed 
as  the  volume— angle  integral  of  the  flux,  i.e.. 


-fOO 

®  (0)  (E)  =fdxfdn  <j>(X,E,M) 
-00 


(14) 


or  integrating  over  d/2  to  give 

+°° 

®  (0)(E)  =  JdX  <^(X,E)  (15) 

-oo 

where  (J)  (X,E)  is  the  all  angle  flux.  Equation  (15)  is  the 
o 

zero  moment  of  the  spatial  distribution  of  the  total  (all 
angle)  flux.  This  "zero  moment  flux"  can  be  calculated 
exactly  by  using  some  properties  of  the  Fourier  integral 
transform  (Ref  3)  on  the  spatial  variable  in  equation  (13) , 


i.e. 


-v  +00 

fy(X,E)  j  =  Jt£  (X,E)e‘1Px  dxS  6j,(P , E) 


(16) 


It  can  be  seen  that  if  equation  (16)  is  written  forj^=0 
and  the  special  case  P  =  0,  then  it  is  equivalent  to  the 
zero  moment  of  equation  (15),  i.e., 
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®  (E)  = 


e  (p,E) 

O 


(17) 


p  =  0 


Therefore,  equation  (13)  is  Fourier  transformed  using  the 
properties  (Ref  3) 


3<t>  (X, E) 
o 


=  iP  6  (P,E) 
o  ’ 


(18) 


^ X 

and 

F  (X)>  =1 

To  obtain 

iP  6-i  (P,E)=£  (E)0  (P,E)=S(E)+/£  (E'-+E)0  (P,E')dE' 

i  t  o  q  so  o 


CO 


(19) 


Equation  (19)  is  valid  for  any  arbitrary  P  but  from  equa¬ 
tion  (17)  the  zero  moment  flux  is  obtained  for  P=0;  there¬ 
fore,  Maclaurin  expansions  for  0Q  (P,E)  and  0^(P,E)  are 
chosen  as 


0  (P,E)=0° 

o x  *  '  o 

0i(P,E)=0° 


(E)  +  (-iP)01(E)  +  (-iP)2  02(E)+... 

o  2!  ° 

2 

(E)+(-iP)  0^ (E)+ ~lP^—  02(E)+... 

1  2!  1 


(20) 

(21) 


Where  the  argument  of  the  Maclaurin  expansion  is  (-iP) 
rather  than  P.  The  primes  denoting  differentiation  in  the 
expansion  have  been  replaced  by  superscripts.  That  is. 


e"  (E)=02  (E)  (22) 

o  o 
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Under  this  notation  the  zero  moment  flux  is  symbolized 

®  °(E)  =  9°  (E)  (23) 

Which  is  the  energy  dependent  flux  of  interest. 

Davison  and  Sykes,  (Ref  4:343)  show  that  the  n-th 
moment  of  the  neutron  flux  can  only  involve  spherical  har¬ 
monics  of  order  n  or  less.  Further  they  show  that  due  to 
the  odd -even  nature  of  the  functions  involved 


oo 

/ Xn  L  (X,E)dx/0 
-oo  ^ 


(24) 


only  when  both  n  andXare  both  even  or  both  odd.  Thus 


the  Maclaurin  expansions  become 

9  (P,E)  =  9° (E)  +  (ziEl2  62(E)  +...  (25) 

o  o  2!  ° 

9  (P,E)  =  (-iP)01(E)  +  (-iP)3  03 (E)  +...  (26) 

1  3! 


By  substituting  equations  (25)  and  (26)  into  equation  (19) 

o 

and  equating  like  powers  of  (-iP)  the  zero  moment  equation 
is  obtained  qq 

£  (E)9°(E)  =  S  (E)  +J  £  (E >  E)  0°  (E '  )  dE '  (27) 

to  o  so 

where  0°(E)  is  the  zero  moment  of  the  neutron  flux  and  is 
o 

the  energy  dependent  flux  necessary  to  calculate  the  group 
cross  sections. 
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Solution  of  the  Zero  Moment  Equation 

Equation  (27)  is  solved  numerically  by  expressing  it 
in  multigroup  notation.  In  group  form  it  is 


+ 

-Eti 


-E 


.+ 


Zt(E)0°(E)dE=  S(E)dE  + 


'dE  /^so(E'-^E)0O(E,)dE' 


n 


(28) 


From  the  definition  of  a  flux  weighted  group  cross  section, 
equation  (5) ,  it  is  seen  that  the  first  term  becomes 
E+ 

i  Zt(E)9°(E)dE  -£”  0°n  (29) 

ill 

n 

<-,n 

Where is  the  macroscopic  total  cross  section  of  the  n-th 

energy  group,  and  0  is  the  total  group  flux.  Similarly 

on 

the  second  term  is 


S(E)dE  =X 
n 


(30) 


Where  V  is  the  fraction  of  the  source  emitted  in  the  n-th 
n 

group,  providing  S (E)  is  a  normalized  source.  The  last 
term 


(E ' — ^E)  0  °  (E '  )  dE ' 
o 


(31) 
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is  considered  as  a  group  of  n  double  integrals  where  the 
integral  from  0  toaoon  dE*  has  been  broken  into  n  finite 
intervals;  /\E|,  ^E^,  /\E^,  ...  over  the  energy  range  .4139 
ev  to  14.918  MeV,  It  is  important  that  the  integration  on 
the  variable  E  be  performed  first  since  the  limits  of  in¬ 
tegration  on  E,  i.e.,  E^  to  E+,  are  actually  functions  of 
E'  and  the  maximum Ae  of  the  scattering  nucleus.  The  last 

term  becomes 

+ 

E 

fn  dE f  y  (E' — *E)  0  (E')dE'  + 

JE“  J&E'  Lso  ° 

n  1 


E  K 

f  "<“[  Z  so(E,->E)0°(E,)dE'...rn_dEr  £  (E’-^E)0°(E'  )  dE' 

Je;  JAe-  J  En 

(32) 


Each  double  integral,  say  the  j-th,  in  equation  (32)  is  a 
measure  of  the  neutron  scatter  transport  from  the  j-th  to 
the  n-th  group.  Recalling  the  definition  of  a  flux  weight¬ 
ed  cross  section  in  a  slightly  different  form 


=  L-  dE  (E-  Zso(E'-^E)eo<E'>dE' 


(33) 


E+ 

r  j 


_  0  (E' ) dE' 
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It  is  seen  that  equation  (32)  can  be  expressed  as 


(34) 


Substituting  equations  (29) ,  (30) ,  and  (34)  into  question 
(28)  and  rearranging  terms  yields  the  numerical  zero  mo¬ 
ment  equation, 


£  -£So<nH>n) 

t 


(35) 


The  total  group  flux,  0°  ,  is  determined  by  solving  equa 

on 

tion  (35) . 


New  Barnyard  reads  in  a  99  group  cross  section 
set  (or  sets)  from  a  cross  section  library  tape,  which  is 
discussed  in  Chapter  III,  calculates  a  99  group  macrosco¬ 
pic  cross  section  set,  then  flux  calculation  is  begun  with 
group  one  (highest  energy  group)  where  the  scatter  in  term 
is  zero.  Calculation  proceeds  consecutively  through  all 
remaining  groups  down  to  .4139  ev.  It  is  assumed  that  no 
"scatter  up"  in  the  neutron's  energy  occurs  in  this  above- 
thermal  energy  region.  Once  the  total  flux  has  been  cal¬ 
culated  for  each  of  the  99  groups,  it  is  used  to  flux 
weight  the  99  group  macroscopic  cross  section  set  over 
each  of  the  desired  broad  group  limits. 
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III.  The  Code 

In  this  chapter  a  few  comments  about  New  Barnyard  are 
made  followed  by  the  contents  and  structure  of  the  two  data 
tape  sources.  Next,  the  numerical  equations  used  to  calcu¬ 
late  the  group  cross  sections  and  related  constants  are 
given  and  finally  the  special  features  of  the  output  are 
mentioned. 

New  Barnyard  is  written  in  Fortran  IV  for  an  IBM  7094 
digital  computer.  Two  source  decks  are  available  in  order 
to  use  the  GGC-4  cross  section  data  tape  (Ref  5)  and  the 
two  ENDF/B  cross  section  data  tapes  (Ref  6) .  The  GGC-4 
source  deck  (i.e.,  the  deck  that  uses  the  GGC-4  data  tape) 
and  the  ENDF/B  source  deck  are  listed  in  appendices  A  and 
B.  Extra  comment  cards  have  been  positioned  throughout  the 
source  decks  to  make  them  easier  to  read.  Appendix  C  lists 
the  variables  and  the  meaning  or  use  of  each  variable  used 
in  the  two  source  decks.  The  code  reads  in  cross  section 
data  and  then  determines  the  zero  moment  of  the  neutron 
flux  which  the  code  uses  to  calculate  broad  group  cross 
sections  and  other  related  constants  over  the  energy  range 
.4139  ev  to  14.918  MeV. 
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The  Cross  Section  Data  Tapes 

As  mentioned  earlier,  either  the  GGC-4  data  tape 
or  the  ENDF/B  data  tapes  can  be  used  in  New  Barnyard.  Both 
data  tape  sources  have  been  compiled  recently  (1966-1967) 
and  the  results  for  absorption  or  total  cross  sections  ob¬ 
tained  from  the  two  source  decks  for  the  same  problem  are 
within  approximately  5  to  10%  of  each  other.  The  ENDF/B 
data  is  updated  from  time  to  time  by  R.S.I.C.  (The  Radiation 
Shielding  Information  Center)  at  Oak  Ridge  National  Labora¬ 
tory.  At  this  writing  they  are  also  preparing  data  for 
nuclides  other  than  those  that  are  contained  on  the  ENDF/B 
data  tapes  described  here.  New  data  tapes  can  be  obtained 
through  R.S.I.C.* 

Both  data  tape  sources  contain  total,  absorption, 
fission**  and  scattering  transfer  microscopic  cross  sections 
for  99  energy  groups.  The  99  "fine"  group  structure  is 
shown  in  Appendix  D.  The  P0  transfer  cross  sections  on  the 
GGC-4  tape  are  listed  separately  as  elastic,  inelastic, 
n-2n,  and  total  transfer.  The  ENDF/B  tapes  list  only  a  P0 
transfer  cross  section  which  is  the  sum  of  the  elastic,  in¬ 
elastic  and  2(n-2n)  transfer  cross  sections.  The  n-2n 

*  ENDF/B  data  tapes  are  obtained  from  R.S.I.C,  Oak  Ridge 
National  Laboratory,  Post  Office  Box  X,  Oak  Ridge, 
Tennessee,  37830 

**  The  ENDF/B  data  tapes,  list VOf  for  99  groups.  The  GGC-4 
data  tape  list!/  and  Of  separately. 
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transfer  cross  section  is  multiplied  by  2  to  conserve  neu¬ 
trons.  It  is  assumed  that  both  neutrons  are  emitted  in 
the  same  energy  group.  The  GGC-4  tape  has  PO  through  P3 
elastic  scattering  transfer  cross  sections.  The  ENDF/B 
tapes  have  PO  through  P8  elastic  scattering  transfer  cross 
sections;  however,  the  PO,  as  mentioned  above,  also  in¬ 
cludes  inelastic  and  n-2n  transfer  cross  sections.  Both 
data  sources  include  only  PO  transfer  cross  sections  for 
the  n-2n  and  inelastic  reactions.  The  GGC-4  tape  also 
contains  other  information  such  as  the  99  energy  group 
boundaries,  fission  source  spectrums,  one  dimensional  cross 
section  arrays  (for  (n,0t)  ,  (n,tf)  and  other  similar  reac¬ 
tions)  and  resonance  data.  The  99  group  fission  source 
spectrums  for  several  nuclides  are  listed  in  Appendix  F. 

The  GGC-4  tape  has  data  for  the  following  nuclides*: 


lc 

Hydrogen 

7. 

Boron (natural) 

2. 

Deuterium 

8. 

Boron- 10 

3. 

Helium 

9. 

Carbon 

4. 

Lithium-6 

10. 

Nitrogen 

5. 

Lithium-7 

11. 

Oxygen 

6. 

Beryllium 

12. 

Sodium 

*  It  was  mentioned  that  data  for  45  nuclides  was  on  the 
GGC-4  data  tape;  this  is  true  but  data  for  some  nuclides 
is  repeated.  This  repeated  data  was  obtained  by  Gulf 
General  Atomic  from  other  sources. 
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13. 

Magne  s ium 

14. 

Aluminum 

15. 

Silicon 

16. 

Sulfur 

17. 

Calcium 

18. 

Titanium 

19. 

Chromium 

20. 

Manganese 

21. 

Iron 

22. 

Cobalt 

23. 

Nickel 

24. 

Copper 

25. 

Molybdenum 

The  ENDF/B  tapes  have  data 


1. 

Hydrogen 

2. 

Deuterium 

3. 

Lithium-6 

4. 

Lithium-7 

5. 

Beryllium 

6. 

Boron- 10 

7. 

Carbon 

8. 

Nitrogen 

9. 

Oxygen 

10. 

Sodium 

11. 

Magnesium 

12. 

Aluminum 

13. 

Titanium 

14. 

Chromium 

15. 

Manganese 

26. 

Cadmium 

27. 

Tungsten (natural) 

28. 

Tungsten-180 

29. 

Tungsten -18 2 

30. 

Tungsten-183 

31. 

Tungsten-184 

32. 

Tungsten-186 

33. 

Lead 

34. 

Uranium-233 

35. 

Uranium-235 

36. 

Uranium-238 

37. 

Plutonium- 241 

for  the  following  nuclides: 

16.  Iron 

17.  Nickel 

18.  Tungsten-182 

19.  Tungsten-183 

20.  Tungsten-184 

21.  Tungsten-186 

22.  Vanadium 

23.  Uranium-235 

24.  Uranium-238 

25.  Plutonium-238 

26.  Plutonium-239 

27.  Plutonium-240 

28.  Pliitonium-241 

29.  Plutonium- 24 2 


Tables  XIII  and  XIV  show  the  structures  of  the  data  tapes 
and  Table  XV  lists  comments  about  these  tapes.  Appendix  E 
contains  these  three  tables. 

Each  of  the  99  group  cross  section  sets  contained  in 
the  two  data  tape  sources  was  calculated  by  flux  weighting 
energy  dependent  cross  sections  with  a  1/E  flux  dependence 
over  the  limits  of  each  fine  group.  The  specific 
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calculations  of  these  99  group  cross  section  sets  can  be 
found  in  the  description  of  the  GGC-4  code  (Ref  5) . 

Since  repeated  reference  is  made  to  PN  (PO,  Pi,  etc.) 
scattering  transfer  cross  sections,  a  review  of  what  is 
meant  by  the  PN  elastic  scattering  transfer  cross  sections 
is  included  here.  The  differential  scattering  transfer 
cross  sections  for  elastic  scattering  of  neutrons, 

Cfc  (E *E,U  )  are  experimentally  determined  and  are  us- 

u  O 

ually  tabulated  as  a  plot  of  (J  (E'—>E,llQ)  versus where 

IX  is  the  cosine  of  the  scattering  angle.  BNL  400  (Ref  7) 
o 

shows  typical  differential  scattering  transfer  cross  sec¬ 
tion  curves  for  nuclides  with  Z  numbers  from  1  to  22.  In 
order  to  cut  down  on  the  amount  of  data  that  would  have  to 
be  retained  for  each  of  these  curves,  a  Legendre  polynomial 
expansion  of  the  differential  scattering  transfer  cross 
section  is  performed  as  follows: 

cf  (E'— )E,4d)  =  Y1  2n±l  Cfsk  (E<— >E)P  (Mo)  (36) 

3  fci  n 

P  (fj.  )  is  thek-th  Legendre  polynomial  and  (E'— ^E)  is 

k  ^  sk 

the  k-th  scattering  transfer  cross  section  coefficient. 

These  coefficients,  Cfsk(E'— ^E),  can  be  calculated  equat¬ 
ing  equation  (36)  to  the  experimental  values  of  Cf^i E'-^E,^). 
The  PN  fine  group  scattering  transfer  cross  sections  are  then 
calculated  by  weighting  coefficients  Cfsk(E’-^E)  with  1/E' 
over  the  limits  of  each  fine  group;  that  is, 
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PN(i->j)  = 


i-1 


(37) 


k  =  0,  1,  2,  . . .K 
i  =  2,  3,  ...  99 
j  -  2,  3,  ...100 


where  N  =  k.  This  polynomial  expansion  of  (J  (E* — >E,/4  ) 

u  O 


can  lead  to  negative  cross  section  values  when  the  number 

of  terms  (K+l)  used  in  the  expansion  is  small.  Thus,  there 

is  an  advantage  of  having  the  ENDF/B  data  since  it  has  P0 

through  P8  elastic  scattering  transfer  cross  section 
coefficients. 

The  Numerical  Equations  Used  in  the  Gode 

After  the  code  reads  the  99  group  set  (or  sets) 
of  microscopic  cross  sections  specified  for  a  particular 
problem,  it  calculates  a  99  group  set  of  macroscopic  cross 
sections  using  the  relation 


NNUK 


(36) 


Where  ($  is  the  microscopic  cross  section  in  bams  for  the 


i-th  nuclide  and  the  n-th  fine  group,  N^is  the  number  den¬ 
sity  in  nuclei/cm-bam  for  the  i-th  nuclide,  NNUK  is  the 
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number  of  nuclides  for  the  problem.  The  code  then  uses 
the  fine  group  total  macroscopic  cross  sections  and  the 
fine  group  PO  total  transfer  macroscopic  cross  sections 
and  calculates  the  total  flux  for  each  of  the  99  groups 
from  the  zero  moment  equation  (eq.  (35))  which  is  repeated 
here  for  convenience: 


o 

0  = 
on 


1 

-  Z  fr-*1)) 

t  so 


n-1 


£so(j^n)ej 


n=l,  2,  3,  ...99 


(35) 


Next,  the  code  flux  weights  the  99  group  macroscopic  cross 
section  set  over  the  limits  of  each  of  the  broad  groups. 
The  specific  cross  sections  calculated  along  with  the  nu¬ 
merical  equations  used  to  calculate  them  are  listed  below: 

(In  the  following  equations,  the  subscripts  n  and  j  denote 
fine  groups  (any  one  of  the  99  groups)  and  K  and  L  denote 
broad  groups.) 


1. 


Absorption  cross  section 

V  '  o 
K  ) 


z  =  ^ 


K  =  1,  2,  3, 


* - 1 


a 


0^ 


n  £K 


. . .NBBG 
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2.  Total  cross  section 


3.  Fission  cross  section 


K=l,  2,  3,  ...NBBG 


K=l,  2,  3,  ...NBBG 


4.  Nu*  (Fission  cross  section) 

Cl/1  Ztt 

_  K  '  f 

<i/£c>  - 


'f  ^  o 

^  c 

n£K 


K-l,  2,  3,... NBBG 


5.  Scattering  transfer  cross  section 


PN  (K— *-L)  = 


0  /  PN  (n— »j) 

. — on  4,  ‘ 
n£K 


K-l,  2,  3,... NBBG 


e 


ng* 


07? 


Where  NBBG  is  the  number  of  broad  groups, 


0o^  is  the  total  flux  for  the  n-th  group, 


is  the  summation  symbol  denoting  that  the  sum  is 
over  all  the  fine  groups  in  the  broad  group  K, 
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y  is  the  fraction  of  the  source  neutrons  emitted 
n 

in  the  n-th  group  (source  normalized  to  1) , 
is  the  PO  total  scattering  transfer  macroscopic 
cross  section  for  a  neutron  which  scatters  from 
group  j  to  group  n, 

is  the  macroscopic  absorption  cross  section  for 
the  n-th  group, 

is  the  macroscopic  total  cross  section  for  the 
n-th  group, 

is  the  macroscopic  fission  cross  section  for  the 
n-th  group, 

is  the  average  number  of  fission  neutrons  pro¬ 
duced  per  fission  event  in  the  n-th  group, 

PN(n— *j)  is  equal  to  (2N+1)  V  (n-^>j)  where  N=0,  1,  2,  3 

SN 

for  the  GGC-4  data  tape  and  N=0,  1,  2,  ...8  for 
the  ENDF/B  data  tapes, 

LSN(n^J)  is  the  scattering  transfer  macroscopic  cross 

section  coefficient  for  a  neutron  which  scatters 
from  group  n  to  group  j .  The  subscript  N  denotes 
the  N-th  scattering  transfer  cross  section  co¬ 
efficient  (term)  in  the  Legendre  polynomial  ex¬ 
pansion  of  the  differential  scattering  transfer 
cross  section.  For  N=0,  this  can  be  the  elastic, 


Z 


n 


,n 


n 


\)n 
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inelastic,  n-2n,  or  the  total  scattering  transfer 
cross  section  when  using  the  GGC-4  data  tape  but 
only  the  total  scattering  transfer  cross  section 
when  using  the  ENDF/B  data  tapes. 

The  code  also  calculates  the  macroscopic  transport  cross 
section,  the  diffusion  coefficient,  and  the  average  cosine 
of  the  scattering  angle  for  each  broad  group.  These  three 
constants  are  defined  as  follows: 

transport  cross  section,  T  =  T  -  T 

tr  t  si 


1 

diffusion  coefficient,  Ds  ~ — ^ - 

3*£ 

tr 

average  cosine  of  the  scattering  angle, 

so 

The  multigroup  calculations  performed  in  New  Barnyard  for 
these  constants  to  obtain  broad  group  constants  are: 

1.  Transport  cross  section 


K  K 

£  =  £-Pl(K*->K)/3  K-l,  2,  3,  ...NBBG 

tr  t 

(Recall  that  the  PN  coefficients,  defined  in 
both  data  tapes,  include  the  factor  2n+1.) 

2.  Diffusion  coefficient 


K  1 

D  -  -  K=1 ,  2,  3,  ...NBBG 

tr 


3.  Average  cosine  of  the  scattering  angle 
_  K 

jJL  -  (Pi  (K— >K) / 3)  K=l,  2,  3,... NBBG 

PO  (K— ^K) 
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Output  Features  of  the  Two  Source  Decks 

Both  source  decks  calculate  total,  absorption,  and  fis¬ 
sion  cross  sections  for  each  broad  group  as  well  as  the 
transport  cross  section,  diffusion  coefficient,  and  the  av¬ 
erage  cosine  of  the  scattering  angle  for  each  broad  group. 

The  ENDF/B  source  deck  can  also  provide  PO  through  P8 
macroscopic  elastic  scattering  transfer  broad  group  cross 
sections;  however,  the  PO  transfer  cross  sections  include 
inelastic  and  n-2n  transfer  cross  sections  as  mentioned 
earlier.  The  order  of  the  PN  transfer  cross  sections  de¬ 
sired  is  specified  in  the  input  data.  The  maximum  number 
of  broad  groups  is  20  due  to  the  fact  that  insufficient 
core  storage  occurrs  when  more  than  20  broad  groups  are 
used. 

The  GGC-4  source  deck  also  provides  a  list  of  the  nu¬ 
clides  on  the  data  tape,  the  99  fine  group  structure,  the 
broad  group  structure,  the  PO  elastic,  inelastic,  n-2n, 
and  total  transfer  broad  group  cross  sections,  and  it  lists 
PO  through  P3  elastic  scattering  transfer  broad  group  cross 
sections.  This  deck  can  calculate  group  cross  sections  for 
22  broad  groups. 
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IV.  Operating  Instructions 

This  chapter  contains  the  information  necessary  to  use 
New  Barnyard  (either  source  deck)  on  an  IBM  7094  computer. 
The  input  data  cards  required  for  each  source  deck  are  giv¬ 
en.  Next,  the  control  cards  required  by  the  IBM  7094  are 
shown,  and  finally  the  composite  deck  (control  cards,  source 
deck,  data  cards)  is  described. 

The  user  should  consider  the  following  factors  when 
selecting  which  source  deck  to  use: 

1.  The  GGC-4  source  deck  (i.e.,  the  deck  that  uses 
the  GGC-4  data  tape)  runs  approximately  twice  as 
fast  on  the  computer  as  the  ENDF/B  source  deck. 

2.  Results  for  absorption  or  total  cross  sections 
from  the  two  source  decks  for  the  same  problem 
are  approximately  within  5  to  10%  of  each  other. 

3.  P0  through  P3  elastic  scattering  transfer  cross 
sections  can  be  obtained  from  the  GGC-4  source  deck 

PO  through  P8  elastic  scattering  transfer  cross 
sections  can  be  obtained  from  the  ENDF/B  source  deck. 

4.  The  PO  output  of  the  GGC-4  source  deck  lists  sep¬ 
arately  the  elastic,  inelastic,  2(n-2n),  and  total 
transfer  cross  sections.  The  PO  output  of  the 
ENDF/B  source  deck  lists  only  the  total  transfer 
cross  sections. 
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As  a  convenience  to  the  user  of  New  Barnyard,  sample  pro¬ 
blems  including  input  data  cards  and  computer  output  from 
both  source  decks  are  presented  in  Appendix  G. 


CARD  NUMBER  1 

Format  (18A4) 

Columns  1-72 

Problem  description  -  Anything  can 

be  punched  on  this  card.  This  in¬ 
formation  will  appear  at  the  top  of 

each  page  of  results. 

Symbol* 

BXCX(T),  1=1,  18 

CARD  NUMBER  2: 

Format  (313) 

Columns  1-3 

The  number  of  broad  energy  groups 

Symbol 

NBBG,  NBBG&22 

Columns  4-6 

The  number  of  nuclides  in  the  pro¬ 
blem 

Symbol 

NNUK 

Columns  7-9 

This  value  is  either  1  or  0.  If  it 

is  1,  the  code  will  calculate  the 

flux.  If  it  is  0,  the  code  will 

have  a  flux  input. 

Symbol 

KKK 

*  The  symbol  (variable)  used  in  the  code  for  this  speci¬ 
fied  data. 
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CARD  NUMBER  3; 

Columns  1-3 
Columns  4-6 


etc. 


Format  (2213) 

Lower  broad  group  boundaries -The 
boundaries  are  the  numbers  of  the 
lowest  fine  group  in  each  of  the 
selected  broad  groups.  For  example, 
if  one  of  the  lower  energy  bounaries 
was  10.00  MeV,  the  value  4  would  be 
specified.  The  number  of  values 
punched  on  this  card  is  equal  to  the 
number  of  broad  groups,  NBBG. 


Symbol  LBGB  (I),  I  =  1,  NBBG 

Each  of  the  actual  energy  boundaries  selected  for  a  parti¬ 
cular  problem  must  correspond  to  one  of  the  99  energy  groups 
therefore,  the  broad  group  boundaries  should  be  selected 
from  the  99  fine  group  structure  (Appendix  D) . 

Sometimes  it  is  desirable  to  know  the  broad  group  scatter¬ 
ing  transfer  cross  section  for  scatter  into  a  "thermal 
dump"  group.  This  can  be  done  in  New  Barnyard  by  specify¬ 
ing  a  broad  group  for  the  energy  range  .4139  ev  to  0.  The 
lower  broad  group  boundary  value  that  would  be  input  for 
this  group  would  be  100.  This  procedure  is  valid  because 
both  data  sources  give  scattering  transfer  cross  section 
values  for  scatter  into  this  "100th"  group;  however,  one 
should  be  aware  of  the  fact  that  no  100th  group  total 
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or  absorption  cross  section  values  are  included  on  the  data 
tapes . 


CARD  NUMBER  4: 

Format  (1F12.7,  1E13.6) 

Columns  1-12 

Identification  number  of  the  nuclide 

on  the  data  tape  -  The  I.D.  number 

for  the  nuclides  are  shown  in  Table  I. 

Symbol 

AID  (IXL) ,  IXL  =  1,  NNUK 

Columns  13-25 

The  number  density  (nuclei/cm-barn) 

Symbol 

of  the  nuclide 

24  2 

(nuclei/cm-barn)  =  nuclei  *  10”  cm 

cm3  barn 

DENT (IXL),  IXL  =  1,  NNUK 

Card  4  is  repeated  for  each  nuclide  in  the  problem.  All 
card  4's  precede  card  5.  It  is  necessary  to  start  with  the 
card  with  the  lowest  I.D.  number,  followed  by  the  cards  with 
increasing  I.D.  numbers,  since  the  nuclides  on  the  data  tape 
occur  in  the  order  of  increasing  I.D.  numbers. 


CARD  NUMBER  5: 

Format  (6E12.6) 

Columns  1-12 

Either  a  flux  spectrum  or  a  source 

Columns  13-24 

• 

• 

spectrum  is  specified  depending  on 

the  value  on  card  number  1  for  KKK 

• 

etc. 

(columns  7-9) .  99  values  (6  per 

card)  must  be  specified. 

Symbol 

FLUX  (I)  or  SSSS(I),  1=1,99 
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If  a  flux  spectrum  is  included  as  input  data,  none  of  the 
values  can  be  zero.  Appendix  F  contains  99  group  source 
spectrums  for  fissionable  nuclides.  The  source  spectrum 
used  must  be  normalized  to  1  as  is  the  case  for  the  spectrums 
in  Appendix  F. 

Preparation  of  Input  Data,  Source  Deck  for  ENDF/B  Data  Tapes 

The  first  three  input  data  cards  are  the  same  as  the 
data  cards  for  the  GGC-4  source  deck  with  the  exception 
that  the  maximum  allowable  number  of  broad  groups  is  20 
(i.e.,  NBBG*20). 

Card  4  is  repeated  for  each  nuclide.  All  Card  4's  precede 
Card  5.  It  is  necessary  to  start  for  the  card  with  the 
lowest  material  number  followed  by  the  cards  with  increas¬ 
ing  material  numbers,  since  the  nuclides  on  each  of  the  two 
data  tapes  occur  in  the  order  of  increasing  material  num¬ 
ber.  If  both  data  tapes  ar$  being  used  for  a  particular 
problem  it  is  only  necessary  to  have  the  cards  arranged  such 
that  the  material  number  (nuclides)  that  are  on  the  same 
tape  be  placed  in  increasing  order. 

CARD  NUMBER  4:  Format  (1x,1A4,2I2,1I6,1E13.4) 

Columns  2-5  Material  number  of  the  nuclide  on 

the  data  tape-The  material  numbers 

for  the  nuclides  are  shown  in  Tables 

II  and  III. 
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Symbol 

MATN0 

Columns  6-7 

Order  of  the  PN  scattering  transfer 

cross  sections  desired-The  values 

allowed  are  0  through  8.  This  value 

must  be  the  same  for  each  nuclide 

used  in  the  problem. 

Symbol 

LORDER 

Columns  8-9 

Logical  unit  number  of  the  data  tape- 

If  using  data  tape  B  (Table  jj  .)  the 

input  value  will  be  1.  If  using 

data  tape  C  (Table  III )  the  input 

Symbol 

value  will  be  2.  (No  decimals  are  used 
with  these  values  since  I  format) 

N 

Columns  10-15 

The  number  of  data  records  before 

reaching  the  nuclide  of  interest. 

Symbol 

NOR 

Columns  16-28 

The  number  density  (nuclei/cm-barn) 

of  the  nuclide 

nuclei  =  nuclei  *  10"  4cm  | 

■i  "3  t  1 

l  cm- bam  cm  barn 

Symbol 

DENT 

The  number  of  data  records  before  reaching  the  nuclide  of 
interest  is  obtained  with  the  aid  of  the  "No.  of  data  Re¬ 
cords"  columns  in  Tables  II  and  III.  For  example,  if  C-12 
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TABLE  I 


Data  Tape  A,  GGC-4  Data 

Master  Duplicate 

Tape  No. 

2602 

0S122 

Block  Size* 

256 

25 6 

Track 

7 

7 

Density 

556BPI  556BPI 

Nuclide 

Identification 

Nuclide 

Number 

Description 

1.0000000 

Hydrogen 

1.2000000 

Deuterium 

2.0000000 

Helium 

3.0062000 

Lithium-6 

3.0072000 

Lithium- 7 

4.0000000 

Beryllium 

5.0000000 

Boron 

5.0099999 

Boron- 10 

6.0200000 

Carbon 

7.0000000 

Nitrogen 

8.0200000 

Oxygen 

11.0000000 

Sodium 

12.0000000 

Magnesium 

13.0000000 

Aluminum 

14.0000000 

Silicon 

16.0000000 

Sulfur 

20.0000000 

Calcium 

22.0000000 

Titanium 

23.9999990 

Chromium 

25.0000000 

Manganese 

26.0000000 

Iron 

27.0000000 

Cobalt 

27.9999990 

Nickel 

29.0000000 

Copper 

41.9999990 

Molybdenium 

48.0000000 

Cadmium 

73.9999990 

Tungsten 

74.1799990 

Tungsten  -  180 

74.1820000 

Tungsten  -  182 

*  Binary  tape 
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Nuclide 

Identification 

Number 


74.1821000 

74.1829990 

74.1831000 

74.1840000 

74.1841000 

74.1859990 

74.1861000 

82.0000000 

92.2334990 

92.2350000 

92.2351990 

92.2379990 

92.2381000 

92.2381990 

92.2382990 

94.2411990 


Table  I  (Contd') 


Nuclide 

Description  P-N 


Tungsten-182  (Resonance  3 

data  included) 

Tungsten-183  3 

Tungsten-183  (Resonance  3 

data  included 

Tungsten- 184  3 

Tungsten-184  (Resonance  3 

data  included) 

Tungsten-186  3 

Tungsten-186  (Resonance  3 


data  included) 

Lead 

Uranium-233 

Uranium-235  (NASA  Data) 
Uranium-235 

Uranium-238  (NASA  Data) 
Uranium-238  (Resonance 
data  included) 

Uranium-238  (ORNL  data)  3 

Uranium-238  (ORNL  Re so-  3 

nance  data  included) 
Plutonium- 241  3 
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TABLE  II 

Data  Tape  B,  ENDF/B  Data 


Master 

Duplicate 

Tape  No. 

1152 

2601 

Block  size* 

80 

80 

Track 

7 

_ 7_ 

Density 

800 BP I 

800 BP I 

Material 

Number  Material  P-N 


No.  of  Data 
Records 


1003 

H-2 

8 

1005 

Li-6 

8 

1006 

Li-7 

8 

1007 

Be 

8 

1009 

B-10 

8 

1010 

C-12 

8 

1012 

N-14 

8 

1013 

0-16 

8 

1014 

Mg 

8 

1015 

Al-27 

8 

1016 

Ti 

8 

1017 

V 

8 

1018 

Cr 

8 

1019 

Mn 

8 

1020 

Fe 

8 

1021 

Ni 

8 

1044 

U-235 

8 

1047 

U-238 

8 

2892 

1537 

1531 

1479 

1447 

1392 

1072 

991 

1048 

1079 

1081 

1136 

963 

1110 

1079 

975 

1197 

1172 


*BCD  Tape 


Total 

Records 

2892 

4429 

5960 

7439 

8886 

10278 

11350 

12341 

13389 

14468 

15549 

16685 

17648 

18758 

19837 

20812 

22009 

23181 
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TABLE  HI 

Data  Tape  C,  ENDF/B  Data 


Master 

Duplicate 

Tape  No. 

1148 

0S578 

Block  size* 

80 

80 

Track 

7 

7 

Density 

800BPI 

800BPI 

Material 

Number 

Material 

P-N 

No.  of  Data 
Records 

Total 

Recor< 

1001 

H-l 

8 

8286 

8286 

1050 

Pu-238 

8 

1133 

9419 

1051 

Pu-239 

8 

1160 

10579 

1053 

Pu-240 

8 

1157 

11736 

1054 

Pu-241 

8 

1215 

12951 

1055 

Pu-242 

8 

1151 

14102 

1059 

Na-23 

8 

1073 

15175 

1060 

W-182 

8 

1139 

16314 

1061 

W-183 

8 

1148 

17462 

1062 

W-184 

8 

1151 

18613 

1063 

W-186 

8 

1173 

19786 

*BCD  Tape 


34 


GNE/PHYS  69-8 


and  N-14  were  the  nuclides  needed  for  a  particular  problem, 
the  value  8886  would  be  punched  on  the  data  card  for  C-12, 
and  the  value  0  would  be  punched  on  the  data  card  for  N-14. 
When  the  computer  reads  the  value  8886  it  "skips"  8886  data 
records;  that  is,  the  data  tape  is  advanced  until  the  nu¬ 
clide  of  interest  is  reached  and  then  the  computer  begins 
to  read  the  data  from  the  tape.  When  the  computer  reads  0, 
no  data  records  are  skipped.  If  two  nuclides  of  interest 
are  separated  by  other  nuclides,  the  number  of  data  records 
that  have  to  be  skipped  is  the  sum  of  the  records  for  the 
nuclides  separating  the  two.  F0r  example,  if  H-2,  and  Be  were 
the  nuclides  of  interest,  3068  (1537+1530  would  be  specified 
on  the  data  card  for  H— 2  and  0  would  be  specified  on  the  data 
card  for  Be. 

CARD  NUMBER  5:  Format  (6E12.6) 

Columns  1-12  Either  a  flux  or  source  spec,  is  speci¬ 

fied.  (Same  as  the  GGC— 4  source  deck) 
Symbol  FLUX ( I )  orSSSS(I),  1=1,99 

Control  Cards  for  the  IBM  7094 

A  brief  explanation  of  the  control  cards  used  in  the 
two  source  decks  is  given  here  to  insure  proper  usage  of 
these  source  decks  on  an  IBM  7094  digital  computer  located 
at  the  Digital  Computation  Division  (ASNCD) ,  Wright-Patter- 
son  AFB,  Ohio.  As  a  convenience  to  the  user,  the  only  con¬ 
trol  card  that  will  have  to  be  prepared  is  the  $J0B  card, 
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the  first  card  of  the  composite  deck;  all  other  control 
cards  have  been  prepared  and  positioned  in  the  source  decks. 
The  user  should  be  aware  of  the  fact  that  control  cards  are 
subject  to  change  and  that  each  computer  facility  has  its 
own  characteristic  control  cards.  The  user  should  always 
have  on  hand  the  current  literature  for  the  computer  being 
used. 


Control  Cards,  GGC-4  Source  Deck 

The  control  cards  used  with  the  GGC-4  source  deck  are 


as  follows: 

Columns  1 

8 

16  31 

$JOB 

Priority,  Time,  Lines  Job  I.D. 

$SETUP 

1 

OS 12 2, NORING 

SlBJOB 

MAP, F IOCS 

$IBMAP 

FILES 

ENTRY 

.UN01. 

. UN01. 

PZE 

UNI TO 1 

. UNIT01 

FILE 

, , READY , INPUT , BIN , BLK= 256 

ENTRY 

.UN02 . 

.UN02. 

PZE 

UNI TO 2 

UNI  TO  2 

FILE 

, , READY , INPUT , BIN , BLK= 256 

ENTRY 

.UN03. 

.UN03 . 

PZE 

UNI TO 3 

UNIT03 

FILE 

, , READY , INOUT , BIN , BLK= 256 

ENTRY 

.UN04. 

. UN04 . 

PZE 

UNIT04 

UNIT04 

FILE 

, ,  READY , INOUT, BIN , BLK=256 

ENTRY 

.UN07 . 

. UN07 . 

PZE 

UNI TO 7 

UNI TO 7 

FILE 

, , READY , INOUT , BIN , BLK= 256 

ENTRY 

.UN08 . 

.UN08 

PZE 

UNI TO 8 

UNI TO 8 

FILE 

END 

, , READY , INOUT , BIN , BLK- 2  5  6 
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Columns  1  8 

$IBFTC  MAIN 
$IBFTC  REWE 
$IBFTC  ONEE 
$IBFTC  CSAVE 
$DATA 
$EOF 

The  $JOB  card  is  an  orange  card  supplied  by  the  computer 
facility.  The  priority  number  should  always  be  0  unless 
the  user  is  authorized  some  other  number.  The  time  number 
is  the  estimated  7094  time*  needed  in  minutes.  The  line 
number  is  the  estimated  number  of  lines**  that  will  be  out 
put.  This  number  should  include  all  printing  and  card 
punching  that  will  be  done  by  the  7094.  If  either  the 
time  estimate  or  the  line  estimate  is  exceeded,  the  com¬ 
puter  run  will  be  terminated.  The  job  identification 
should  include  the  user's  account  number,  office  symbol, 
and  name.  After  the  user's  name,  there  should  appear  a 
slash  followed  by  the  name  Bridgman.  The  data  tapes  are 
filed  under  this  name  and  it  has  to  be  present  on  the  $J0B 
card  to  use  the  GGC-4  data  tape.  An  example  of  the  $J0B 
card  is  shown: 

$J0B  0,3,2000  68-564-00  AFIT-SE  JONES /BRIDGMAN 


*  An  estimate  of  3  minutes  is  usually  adequate  when 
using  the  GGC-4  source  deck. 

**  An  estimate  of  2000  lines  is  adequate  for  the  GGC-4 
source  deck. 
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Control  Cards,  ENDF/B  Source  Deck 

The  control  cards  for  the  ENDF/B  source  deck  are  al¬ 
most  the  same  as  the  control  cards  for  the  GGC-4  source 
deck.  The  exceptions  are  the  following: 

(1)  No  file  is  defined  for  logical  unit  7 
(only  3  scratch  tapes  are  used) 

(2)  2  $SETUP  cards  are  used  for  the  two 
input  data  tapes. 

(3)  Usually  an  adequate  time  estimate  for  the 
$JOB  card  is  5  minutes. 

Composite  Deck  Setup 

The  composite  deck  setup  for  either  source  deck  is 
shown  below.  The  user  has  to  supply  the  $JOB  card  and 
the  input  data  cards. 

$JOB 

$SETUP  (2  cards  for  the  ENDF/B  source  deck) 

$IBJOB 
$IBMAP  FILES 

ENTRY  . UN01. 


END 

$IBFTC  MAIN 

(Source  Deck  -  includes  all  subroutines) 
$DATA 

(Input  data  cards) 

$EOF 
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The  $SETUP  Card  is  used  to  indicate  the  use  of  a  data  tape. 

The  number  1  on  this  card  refers  to  the  logical  number  of 

the  tape  used  in  the  program  (GGC-4  Source  deck) .  The  num¬ 
ber  0S122  is  the  number  of  the  duplicate  copy  of  the  GGC-4 
data  tape.  The  word  NORING  that  follows  the  tape  number  is 
used  to  indicate  that  the  tape  is  not  to  be  written  on. 

The  $IBJOB  card  as  used  here  is  to  have  a  print  out  (map) 
of  the  storage  location  of  each  variable  used  in  the  pro¬ 
gram.  The  word  FIOCS  is  used  to  cut  down  on  the  core  stor¬ 
age  requirement  of  the  program.  The  $IBMAP  card  and  all 

the  cards  that  follow  up  to  and  including  the  END  card  are 
used  to  define  files  for  the  input  data  tape  and  the  5 
scratch  tapes  used  by  the  program.  The  $IBFTC  cards  are 
required  for  the  main  program  and  the  subroutines  used. 

The  $IBFTC  MAIN  card  precedes  the  GGC-4  source  deck  and 
each  of  the  other  $IBFTC  cards  precede  a  subroutine  in  the 
source  deck.  The  names  REWE,  ONEE,  and  CSAVE  that  appear 
on  these  cards  are  arbitrary;  any  name  (6  letters  or  less) 
other  than  the  actual  name  of  the  subroutine  can  be  used. 
The  $DATA  card  is  required  when  input  data  cards  are  used. 
The  $EOF  card  is  the  end-of-file  card  and  it  is  the  last 
card  of  the  composite  deck. 
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V.  Some  Sample  Results  and  Conclusions 
In  this  chapter  results  of  the  GGC-4  version  of  New 
Barnyard  are  compared  with  results  from  other  cross  section 
codes;  also,  group  cross  sections  calculated  from  the 
ENDF/B  version  of  New  Barnyard  are  compared  with  group 
cross  sections  calculated  from  the  GGC-4  version.  Through¬ 
out  this  chapter  repeated  reference  will  be  made  to  the 
zero  moment  of  the  neutron  flux,  zero  moment  flux,  energy 
dependent  flux,  and  flux.  These  terms  all  mean  the  same 
thing.  The  reader  is  reminded  that  the  GGC-4  version  of 
New  Barnyard  refers  to  the  source  deck  that  uses  the  GGC-4 
data  tape  and,  the  ENDF/B  version  of  New  Barnyard  refers 
to  the  source  deck  that  uses  the  two  ENDF/B  data  tapes. 

The  Zero  Moment  of  the  Neutron  Flux  for  Aluminum 

It  was  pointed  out  in  Chapter  II  that  in  order  to  cal¬ 
culate  group  cross  sections,  the  energy  dependent  flux  had 
to  be  known.  This  flux  is  obtained  in  New  Barnyard  by  solv¬ 
ing  the  zero  moment  equation  (35) .  The  validity  of  this 
flux  calculation  was  checked  by  comparing  New  Barnyard's 
calculation  of  the  zero  moment  flux  for  aluminum  with  Gulf 
General  Atomic's  GAM-1  (Ref  9)  calculation  of  the  zero 
moment  flux  for  the  same  material.  GAM-1  calculates  group 
neutron  cross  sections  and  other  related  constants  such  as 
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age  for  the  energy  range  10  MeV  to  .4139  ev.  GAM-1  uses  a 
68  group  cross  section  library  data  tape.  In  order  to  have 
a  better  basis  for  comparison,  the  GGC-4  version  of  New 
Barnyard  was  modified  slightly  to  use  GAM-1' s  library  data 
tape.  Table  IV  shows  the  results  of  these  two  flux  cal¬ 
culations.  The  ratio  of  the  GAM-1  flux  to  the  New  Barnyard 
flux  for  each  energy  group  is  4.000,  as  shown  in  Table  IV 
This  indicates  that  the  shape  of  the  two  flux  spectrums  are 
identical  but  with  the  GAM-1  flux  having  a  magnitude  4 
times  larger  than  the  New  Barnyard  flux.  This  factor  of 
4  is  unimportant,  however,  since  in  the  flux  weighting 
process  for  the  calculation  of  group  cross  sections  this 
factor  will  always  cancel  and  either  flux  will  produce  the 
same  group  cross  sections. 

Seven  Group  Cross  sections  for  Carbon 

The  validity  of  the  flux  weighting  process  to  calcu¬ 
late  group  cross  sections  was  checked  by  comparing  a  7 
group  macroscopic  cross  section  set  for  carbon  calculated 
by  the  GGC-4  version  of  New  Barnyard  with  a  7  group  macro¬ 
scopic  cross  section  set  for  carbon  calculated  by  the  GGC*4 
code.  The  GGC-4  code  was  run  on  a  CDC  6600  computer  at 
the  A.F.  Weapons  Lab.  The  same  flux  spectrum  was  input  to 
each  of  the  codes.  The  group  cross  sections  obtained  from 
each  of  the  codes  were  exactly  identical.  Table  V  shows 
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the  broad  group  structure  for  this  cross  section  set.  Table 
VI  shows  the  results  from  both  codes  for  the  total  cross 
sections  and  Table  VTI  lists  the  PO  (total  transfer) 
scattering  cross  sections  from  both  codes. 
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TABLE  IV 

Flux  Spectrum  for  Aluminum  Calculated 
in  New  Barnyard  and  GGC-4. 


New  Barnyard  Flux  GAM-1  Flux  Ratio  of  GAM-1 
12  12  Flux  to  New 

Group  (°n  /cm  -group)  (on  /cm  group) Barnyard  Flux 


1 

9.2597 

X 

102 

3.6959 

X 

10-1 

4.000 

5 

2.4437 

X 

10° 

9.7748 

X 

10° 

4.000 

10 

8.9699 

X 

10° 

3.5880 

X 

101 

4.000 

15 

1.8522 

X 

101 

7.4089 

X 

101 

4.000 

20 

1.3666 

X 

101 

5.4664 

X 

101 

4.000 

25 

6,7206 

X 

101 

2.6882 

X 

102 

4.000 

30 

2.3125 

X 

101 

9.2500 

X 

101 

4.000 

35 

2.4929 

X 

101 

9.9715 

X 

101 

4.000 

40 

2.4753 

X 

101 

9.9013 

X 

101 

4.000 

45 

2.4753 

X 

101 

9.9013 

X 

101 

4.000 

50 

2.5750 

X 

101 

9.4999 

X 

101 

4.000 

55 

2.1040 

X 

101 

8.4162 

X 

101 

4.000 

60 

1.7247 

X 

101 

6.8988 

X 

101 

4.000 

65 

1.1993 

X 

101 

4.7972 

X 

101 

4.000 
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TABLE  V 

Broad  Group  Structure  of  the  7  Group  Cross  Section  Set  for 
Carbon. 


Group 

Energy 

(ev) 

Interval 

1 

1.4918  x  10 ! 

to 

1.3498  x  10^ 

2 

1.3498  x  107 

to 

1.2214  x  101 

3 

1.2214  x  10? 

to 

1.1052  x  101 

4 

1.1052  x  10? 

to 

1.0000  x  10' 

5 

1.0000  x  10 6 

to 

6.0653  x  10 j? 

6 

6.0653  x  106 

to 

3.0119  x  10b 

7 

3.0119  x  10 

to 

1.0026  x  10b 

TABLE  VI 


Total  Cross  Sections  for  Carbon  Calculated  in  New  Barnyard 
and  GGC-4 


Group 

New  Barnyard  Total 
Cross  Section 
(cm-1) 

GGC-4  Total 

Cross  Section 
(cm-1) 

1 

1.0730  x  10“?- 

1.0730  x  10_1 

2 

1.1078  x  10_i 

1.1078  x  10-i 

3 

1.1223  x  10_1 

1.1223  x  10-1 

4 

9.8612  x  10"‘ 

9.8612  x  10-f 

5 

1.0006  x  10“ 1 

1.0006  x  10-1 

6 

1.4518  x  10’f 

1.4518  x  10-i 

7 

1.6053  x  10"1 

1.6053  x  10"1 
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TABLE  VII 


PO  Scattering  Transfer  Cross  Sections  for  Carbon  Calculated 
in  New  Barnyard  and  GGC-4 


New  Barnyard  PO  Total  GGC-4  PO  Total 

Group  Transfer  Cross  Section  Transfer  Cross  Section 

(cm’j  (cm"  ) 


FROM 

1 

TO 

1 

3.6984 

X 

10“2 

3.6984 

X 

10l2 

FROM 

1 

TO 

2 

1.3377 

X 

10-2 

1.3377 

X 

10-2 

FROM 

1 

TO 

3 

6.6004 

X 

10-3 

6.6004 

X 

10-3 

FROM 

1 

TO 

4 

4.0538 

X 

10-3 

4.0538 

X 

10  "3 

o 

FROM 

1 

TO 

5 

1.6083 

X 

10-2 

1.6083 

X 

10-2 

FROM 

1 

TO 

6 

3.5487 

X 

10"3 

3.5487 

X 

10-3 

o 

FROM 

1 

TO 

7 

1.0219 

X 

10“2 

1.0219 

X 

10-2 

FROM 

2 

TO 

2 

3.8940 

X 

10"2 

3.8940 

X 

10-2 

FROM 

2 

TO 

3 

1.4444 

X 

10"2 

1.4444 

X 

10-2 

FROM 

2 

TO 

4 

8.1390 

X 

10-3 

8.1390 

X 

10  “3 

o 

FROM 

2 

TO 

5 

1.8409 

X 

10-2 

A 

1.8409 

X 

10-2 

FROM 

2 

TO 

6 

5.6734 

X 

lO-3 

5.6734 

X 

10“3 

FROM 

2 

TO 

7 

7.8967 

X 

lO-3 

7.8967 

X 

10-3 

FROM 

3 

TO 

3 

3.8735 

X 

10-2 

3.8735 

X 

10-2 

FROM 

3 

TO 

4 

1.5851 

X 

10-2 

1.5851 

X 

10-2 

FROM 

3 

TO 

5 

2.3362 

X 

10-2 

2.3362 

X 

10-2 

FROM 

3 

TO 

6 

1.1659 

X 

lO"2 

1.1659 

X 

10-2 

FROM 

3 

TO 

7 

4.9461 

X 

io-^ 

4.9461 

X 

10-3 

o 

FROM 

4 

TO 

4 

3.3846 

X 

10“2 

3.3846 

X 

10  2 

FROM 

4 

TO 

5 

2.6936 

X 

10-2 

2.6936 

X 

10-2 

FROM 

4 

TO 

6 

2.6107 

X 

10-2 

2.6107 

X 

TO  -2 

FROM 

4 

TO 

7 

2.5326 

X 

10-3 

2.5326 

X 

10-3 

FROM 

5 

TO 

5 

4.2742 

X 

10-2 

4.2742 

X 

10-2 

FROM 

5 

TO 

6 

3.5086 

X 

10-2 

3.5086 

X 

10-2 

FROM 

5 

TO 

7 

1.9360 

X 

lO’2 

1.9360 

X 

10-2 

FROM 

6 

TO 

6 

9.3408 

X 

10‘2 

9.3408 

X 

10-2 

FROM 

6 

TO 

7 

5.0676 

X 

io-2 

5.0676 

X 

10"2 

FROM 

7 

TO 

7 

1.3782 

X 

10_1 

1.3782 

X 

10"1 

L 
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Six  Group  Cross  Section  Set  for  Nitrogen 

A  comparison  of  cross  sections  calculated  in  the 
ENDF/B  version  of  New  Barnyard  and  the  GGC-4  version  of 
New  Barnyard  was  made  for  several  problems.  The  6  group 
cross  section  set  for  nitrogen  that  is  presented  here  for 
comparison  is  typical  of  the  cross  section  comparisons  ob¬ 
tained  for  other  nuclides  and  mixtures.  The  cross  sections 
calculated  by  the  GGC-4  version  of  New  Barnyard  were  used 
as  the  standard  for  the  relative  difference  calculations 
between  these  two  cross  section  sets.  Table  VIII  shows  the 
6  group  structure  of  the  cross  sections.  The  absorption 
cross  sections  obtained  from  the  two  versions  of  New  Barn¬ 
yard  and  the  relative  differences  between  these  cross  sec¬ 
tions  are  shown  in  Table  IX  ;  similarly,  Table  x  shows 
the  total  cross  sections  and  their  relative  differences. 
Table  XI  and  xil show  the  comparisons  for  the  PO,  Pi,  P2, 
and  P3  scattering  transfer  cross  sections.  The  PO  cross 
section  is  the  total  transfer  cross  section  which  is  the 
sum  of  the  elastic  inelastic,  and  2  (n-2n)  transfer  cross 
section.  The  Pi  through  P3  cross  sections  are  elastic 
scattering  cross  sections  only. 

By  inspection  of  Tables  ’  IX  and  x  it  is  seen  that 
the  largest  relative  differences  are  13%  and  14%  which 
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occurred  for  the  3rd  and  4th  group  absorption  cross  sec¬ 
tions  respectively.  Even  these  relative  differences  aren't 
too  great  when  one  considers  the  following  factors: 

(1)  The  cross  section  data  tapes  are  from  two 
independent  sources,  namely  GGA  (Gulf  Gen¬ 
eral  Atom)  and  ORNL  (Oak  Ridge  National 
Lab . ) . 

(2)  Each  group  flux  is  dependent  on  the  total  and 
PO  fine  group  cross  sections.  If  there  ex¬ 
ists  a  relative  difference  for  each  group 
cross  section,  then  there  exists  a  relative 
difference  for  each  group  flux.  Since  broad 
group  cross  sections  are  obtained  by  flux 
weighting  fine  group  cross  sections,  it  is 
logical  to  argue  that  the  relative  differences 
build  up  rapidly  since  when  multiplying  or 
adding  variables  that  have  relative  differ¬ 
ences  (errors)  associated  with  them  these 
relative  differences  are  added  together. 

The  comparisons  of  the  PO  scattering  transfer  cross 
section  are  quite  good.  Table  XI  shows  that  most  relative 
differences  are  less  than  5%;  however,  the  relative  differ¬ 
ences  for  the  PO  cross  sections  for  transfer  from  1  to  3, 

1  to  4,  1  to  5,  and  1  to  6,  were  41,  99,  100,  and  100% 
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respectively,  but  these  four  cross  section  values  are  very 
small  in  magnitude  and  most  likely  ORNL  did  not  list  cross 
section  values  as  small  as  GGA's.  The  largest  relative  dif¬ 
ference  in  the  comparisons  between  the  Pi  cross  sections  was 
21%.  It  is  seen  from  Table  XII  that  the  relative  differ¬ 
ences  in  the  comparisons  between  P2  cross  sections  are  good 
except  for  cross  section  values  for  transfer  from  1  to  1 
and  1  to  2.  Note  that  in  the  P3  cross  section  comparison, 
even  larger  relative  differences  than  for  P2  are  common  and 
that  for  3  of  the  cross  section  comparisons  (denoted  by  an 
asterisk)  the  algebraic  signs  differ  between  the  cross  sec¬ 
tions.  All  of  these  large  relative  differences  and  alter¬ 
nating  algebraic  signs  can  be  explained  by  the  fact  that 
ORNL  used  an  8th  order  Legendre  polynomial  expansion  of 
the  fine  group  differential  scattering  transfer  cross  sec¬ 
tions  and  GGA  used  a  6th  order  *Legendre  polynomial  ex¬ 
pansion. 


*  GGA  only  listed  PO  through  P3  cross  sections  on  the 
GGC-4  data  tape. 
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TABLE  VIII 

Broad  Group  Structure  of  the  6  Group  Cross  Section  Set  for 
Nitrogen 


Group 

Energy  Interval 
(ev) 

1 

1.4918 

X 

10  7 

to 

3.0119 

x  106 

2 

3.0119 

X 

10  6 

to 

6.0810 

x  105 

3 

6.0810 

X 

10  5 

to 

1.2277 

x  105 

4 

1.2277 

X 

10  5 

to 

2.6126 

x  103 

5 

2.6126 

X 

103 

to 

4.7851 

x  101 

6 

4.7851 

X 

10  3 

to 

3.9279 

TABLE  IX 


Comparison  Between  Absorption  epos's  Selections  Which  Were 
Calculated  in  the  GGC-4  Version  and  the  ENDF/B  Version  of 
New  Barnyard  for  Nitrogen 


GGC-4  Absorption 

ENDF/B  Absorption 

Relative 

Group 

Cross  Section 

Cross  Section  Difference 

(cm"1) 

(cm-1) 

(%) 

1 

1.8192  x  10"5 

1.7910  x  10-5 

2 

2 

3.9154  x  10"6 

3.6110  x  10-6 

8 

3 

5.3592  x  10“7 

4.6538  x  10-7 

13 

4 

1.2364  x  10"7 

1.4085  x  10"7 

14 

5 

9.2645  x  10"7 

9.8802  x  10"7 

7 

6 

4.3170  x  10~6 

4.5740  x  10’6 

6 
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TABLE  X 

Comparisons  Between  Total  Cross  Sections  Which  Were 
Calculated  in  the  GGC-4  Version  and  the  ENDF/B  Version 
of  New  Barnyard  for  Nitrogen 


Group 

GGC-4  Total 

Cross  Section 
(cm"-*-) 

ENDF/B  Total 
Cross  Section 
(cm  (cm”1) 

Relative 

Difference 

(%) 

1 

8.7827  x  1CT5 

8.4124  x 

10-5 

4 

2 

9.9716  x  10"5 

9.6426  x 

io-5 

3 

3 

1.7344  x  IO-4 

1.5963  x 

10 ’4 

8 

4 

3.4983  x  10-4 

3.4664  x 

io-4 

1 

5 

4.9423  x  1CT4 

4.9424  x 

io'4 

0 

6 

5.3625  x  10-4 

5.3478  x 

h-4 

O 

1 

0 
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The  absolute  value  was  taken  so  that  a  comparison  could  be  made. 
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Conclusions 

The  comparison  of  the  flux  calculated  by  New  Barnyard 
with  the  flux  calculated  by  GAM-1  and  the  comparison  of  the 
group  cross  sections  calculated  by  New  Barnyard  with  the 
group  cross  sections  calculated  by  GGC-4  give  confidence  in 
the  "correctness"  or  validity  of  the  calculations  performed 
by  New  Barnyard;  that  is.  New  Barnyard  does  what  it  is  sup¬ 
posed  to  do,  and  does  it  right.  Some  thought,  however, 
should  now  be  given  to  the  usefulness  of  the  zero  moment  of 
the  neutron  flux  for  flux  weighting  cross  sections. 

Recall  that  the  zero  moment  equation  was  based  on  an  in¬ 
finite  medium  and  a  plane  source  of  infinite  dimensions  at 
the  origin.  Assumptions  such  as  these  are  valid  for  reactor 
calculations  where  one  at  least  crudely  has  a  plane  source 
(fission  neutrons  in  the  core)  and  an  infinite  medium 
(shielding  around  the  core) .  When  considering  one  dimension¬ 
al  thin  shield  transport  calculations  for  neutrons  with  ener¬ 
gies  above  thermal,  where  the  mean  free  path  of  a  neutron  is 
larger  than  the  thickness  of  the  shield,  the  zero  moment  e- 
quation  leads  to  an  inaccurate  flux  spectrum.  Thus,  for  a 
thin  shield  calculation  a  better  approximation  for  the  flux 
is  one  equal  to  the  source  spectrum  for  the  problem  since 
relatively  few  neutrons  are  slowed  down  when  passing  through 
the  shield.  It  was  for  cases  such  as  this  that 
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New  Barnyard  was  programmed  with  the  option  of  inputting  a 
flux  spectrum. 

One  should  also  recall  that  no  resonance  calculations 
are  performed  in  New  Barnyard;  however,  resonances  in  gen¬ 
eral  are  negligible  above  a  few  kev  and  it  is  above  this 
energy  range  that  New  Barnyard  is  primarily  intended. 

Finally,  recall  the  zero  moment  flux  was  used  to  flux 
weight  PO  through  P8  scattering  transfer  cross  sections. 
Cross  section  codes  such  as  GAM-1  and  GGC-4  usually  flux 
weight  the  PN  cross  sections  with  a  corresponding  (j)^  flux 
which  comes  from  the  expansion  of  the  three  dimensional 
flux,  §  (X,  E ,H),  in  a  Legendre  polynomial  series  follow¬ 
ed  by  either  a  B^  (Ref  10)  or  a  P^  (Ref  11)  solution  to  the 
Boltzmann  equation.  In  general',  as  N  becomes  large  (j>^ 
becomes  small;  thus,  flux  weighting  the  PN  cross  sections 
with  (|)^  fluxes  produces  PN  group  cross  sections  smaller 
than  the  corresponding  PN  group  cross  sections  flux  weight¬ 
ed  by  the  zero  moment  flux. 
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APPENDIX  A 

GGC-4  Source  Deck  of  New  Barnyard 

On  the  following  pages  is  listed  the  source  deck  that 
uses  the  GGC-4  cross  section  library  data  tape.  It  is  writ¬ 
ten  in  the  Fortran  IV  language  for  use  on  the  IBM  7094  digi¬ 
tal  computer.  It  consists  of  a  main  program  and  3  separate 
subroutines . 

Input  data  from  data  cards  and  from  the  GGC-4  cross 
section  data  tape  are  "read  in"  in  the  main  program.  Five 
scratch  tapes  are  used  to  store  P0,  Pi,  P2,  P3,  inelastic, 
n-2n,  and  total  P0  cross  section  arrays  until  they  are  need¬ 
ed  in  the  main  program.  The  flux  calculation  and  the  broad 
group  calculations  of  the  transport  cross  section,  diffusion 
coefficient,  and  the  average  cosine  of  the  scattering  angle 
are  performed  in  the  main  program. 

The  subroutine  named  REW  is  used  to  rewind  the  scratch 
tapes  and  to  set  a  double  subscripted  scattering  transfer 
cross  section  symbol  to  zero.  The  subroutine  named  ONE  is 
used  to  calculate  99  group  macroscopic  transfer  cross  sec¬ 
tion  sets.  The  subroutine  named  CSAV  is  called  in  the  main 

program  to  flux  weight  the  cross  sections. 

Extra  comment  cards  have  been  added  to  the  listing  that 
follows  so  that  the  program  will  be  easier  to  read.  A  glos¬ 
sary  of  computer  program  symbols  is  given  in  Appendix  C. 
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c  new  barnyard 

C  MULT  I GROUP  NEUTRON  MACROSCOPIC  CROSS  SECTION  CODE 

Cl  THIS  SOURCE  DECK  USES  THE  CGC--4  CROSS  SECTION  DATA  TAPE 

DIMENSION  QXCX(IS) »NTID(3) » AT ( 90 ) , D AD ( 21 ) »DD<3,45) *  LEN ( A )  , 

15 1  GO ( 100 ) ,TRA! 5 1 0 ) » CTH( 110 ) *ENG( 103 )  *UL( 206 ) » LLEGC ( 22 ) » 

2S IGA (100) .SIGT(IOC) *  SSI  37) .TOT (22) *A305(22) . SI  SO ( 22  *  22 ) » 
3P0(22,22) » XI  MELS (22, 22) ,XN2N<22,22! , ESS  I  40) *FTOT(22) *  AID (31) * 
4DENT (31)  >  S  S  S  5  (  9  9  ) ♦ P 1 ( 22 , 22 ) , P2 ( 22 , 22 ) , R3 ( 22 » 22  ) 

COMMON  TTT  ! 5  05C ) >  T  T ( 100»1C0) ,  LDF  , LD  » LT  ,  NBBG » L3C3 ( 99 ) *FLUX( ICO) 
EQUIVALENCE  ( ESS ( 7 ) » SS ( 1 ) ) 

MG=  1 
M  S  S  -  2 
MT  =  3 
M  T  T  "  4 
MMT-  7 
MMMT-3 

C2  MG,  MSS,  MT  »  MTT  »  MMT  ,  MMMT  ARE  THE  TAPE  NUMBERS  USED  IN 
C  THIS  SOURCE  DECK 

FPL" 1.0 
FNBT" 3 1568.00 

C3  READ  FIRST  DATA  CARD 

READ  (5*2) (3XCX(I) *1=1.19) 

C4  READ  SECOND  DATA  CARD 

READ (5*6 )NB3G  * NNUK ♦ KKK 

2  FORMAT (18A4) 

6  FORMAT (3  13) 

REWIND  MG 

C5  THE  FOLLOWING  CARDS  THROUGH  CARD  MO.  10  ARE  FOR  READING 

C  THE  GCC-4  DATA  T  APE  AND  FOR  RR INTI  NO  OUT  IMPORTANT  INFORMATION 

C  FROM  THE  DATA  TARE 

READ  (MC)  ( N T I D ( I ) » I  - 1  *  3 ) 

N3T"NT I D ( 1 ) 

NEP"NTID(2) 

NGT'MTI D ( 3  ) 

NEV=NEP- 1 
NES=NEP+1 
LNX=NEV-t  6 
NBTC-FNBT+O. 1 

I  F  (  FNBT  .LT.  0.)  NBT  C'-FNBT  ~C  •  1 
I  F ( N3T  -NBTC ) 4 , 5  *4 

4  WR I TE( 6 . 991 ) NBT .NBTC 

901  FORMAT  (1M1//5X.22HFAST  DATA  T/°E  NUMBER  I  6  935H  WAS  LOADED  INSTEAD 
1  OF  TAPE  NUMBER  I6*1X*22M  WHICH  WAS  SPECI-FIED.  // 

25X  *  20HPRC3LEM  TERMINATED. 

CALL  EXIT 

5  READ  (MG)  (AT  (I  )  ,1  1,90 
WRITE (6,1002) ( BXCX ( I ) ,1-1*18) 

Ip02  FORMAT ( 1H1 , 13A4 ) 

WRITE (6, 905) NBT, ( AT ( I ) ,1-1,90) 

905  FORMAT ( 26H0FAST  DATA  T APE  NUMBER  -  16/ 

X23H0TAPE  DESCRIPTION .  // (  IX  ,  1SA4 ) ) 

READ  (MC)NNOT 
ASSIGN  809  TO  NPR1 
IF(FPL) 300,308,750 
3p3  ASSIGN  129  TO  NPR1 
GC  TO  810 
tFO  CONTINUE 
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So4  LR INK-0 

WRITE(6,894)N3T 

894  FORMAT ( /35H0C0NTENTS  OF  FAST  DATA  TAPE  NUMBER  16// 

X5X  , 1  lHN'JCL  IDE  MO.  5X  >  20HMUCL  IDE  DESCRIPTION  /) 
glO  DO  129  IXL=1 *NNCT 

READ (MG)  ( DAD ( I )  .1  =  1*21) 

DD ( 1 > I  XL ) -  DAD ( 19 ) 

DD(2»IXL)=DAD(20) 

DD ( 2 . 1  XL ) “DAD ( 21 ) 

C  DAD ( 19)=NU SLIDE  ID  NO. $DAD ( 20 ) =N0.  OF  RESOLVED  RESONANCES 

C  DAD (21) =N0.  OF  UNRESOLVED  RESONANCES. 

I F ( IXL.GT.16.AND.IXL.LE.30)  GO  TO  129 
GO  TO  NPR1.(809,129) 

809  LR I NK=LR I NK+ 1 

IF (LR INK-46) 805 .807.305 
807  V.'RI  TE  (  6  >  895  )  NBT 
305  CONTINUE 

WRITE (6, 897) DAD (19)  .  ( DAD ( I )  ,1=1,18) 

395  FORMAT ( 3  5H1 CON TENTS  OF  FAST  DATA  TAPE  NUMBER  16// 

X  5X , 1 1HNUCL IDE  NO.  5X  » 20HNUCL I DE  DESCRIPTION  /) 

8^7  FORMAT! 1 X , FI  3. 7 ,5X  ,  18A4  ,  / ) 

1 29  CONTINUE 

C  READ  FAST  DATA  TAPE — RESONANCE  TABLES 

READ ( MG ) ( LEN ( I ) ,1=1.4) 

LBS-LEN ( 4 ) 

READ (MG)  ( TTT ( I ) , 1  =  1  ,LBS) 

C  READ  FAST  DATA  TAPE - ENERG I ES » LETHARG I ES , DELTA  U 

JMM=2*NES+NEP 

READ (MG) ( TRA ( I ) , 1=1, JMM) 

READ ( MG )  ( GTH ( I ) ,1  =  1, NOT) 

READ (MG)  ( NT  I D ( I )  ,1  =  1,1) 

DO  8  4  I = 1 ,  NEP 
ENG ( I ) "TRA(  I  ) 

I X= I +NES 

'JL  (  I  )  =TRA  (  IX) 

84  IZ=I X+NES 

ENG(NES) =TRA ( NES ) 

UL(NES) =  TRA ( 2#NES ) 

DO  8  1=1,99 

8  LBGB ! I )=  0 

C 6  READ  THIRD  DATA  CARD 

RE AD (5, 7) ( LLBGB ( I ) , I=1,NBBG) 

DO  9  I-l.NBSG 

9  L3G3 ( I ) =LLBGB( I ) 

7  FORMAT (2413) 

WRITE (6, 1002) ( B X C X ( I )  ,1  =  1,18) 

WRITE (6, 747) 

WRITE! 6 , 748 )( I  ,  ENG ( I ) , ENG ( I  +  1 ) »UL ( I ! ,UL ( I +1 ) , LBGB  (  I  )  ,  1  =  1 ,50 ) 
74?  FORMAT ( 22H0F I NE  GROUP  STRUCTURE  /80H0GROUP  ENERGY 

1 RVAL ( E . V . )  LETHARGY  INTERVAL  ,39H 

2WER  BD.  GRP  BOUND.  ) 

748  FORMAT ( I5,6X,1PE13.6,4H  TO  1PE13.6 »6X» 1PE13. 5 »4H  TO  1PE12.5, 

X 1 1  2  2  ) 

WRITE(6,i002) (BXCX< I ) ,1=1,18) 

WRITE (6, 747) 

WRI TE( 6 ,74  8 )  (  I  »ENG ( I ) , ENG ( I  +  1),UL(I),UL(I+1) , LBGB (  I ) , I =51 ,99 ) 
WRITE (6,1002) ( B  X  C  X ( I ), 1  =  1,18) 
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WRITE (6. 1009) 

1009  FORMAT (22HOBROAD  GROUP  STRUCTURE  /46H0BR0AD  GROUP  ENER 

XGY  INTERVAL ( E.V. ) 

ENG ( 101 ) =0.0 
DO  1008  1  =  1*  NBBG 
I F ( I  .NE.  1)  GO  TO  1005 
JNY=LBGB ( I )+l 

WRITE) 6 *1007) I *ENG( 1) • ENG ( JNY ) 

GO  TO  1008 
1005  K=I-1 

JNYN=LBGB(K)+1 
JNYY=LBGB ( I )+l 

WRITE) 6, 1007) I .ENG) JNYN) *  ENG ( JNYY ) 

1008  CONTINUE 

1007  FORMAT) I  5 , 10X » 1PE 1 3 . 6 , 4H  TO  1PE13.6) 

C  SKIP  FISSION  SPECTRA  SOURCE  DATA  ON  TAPE 

NSP=NT I D ( 1 ) 

DO  10  1  =  1  *  NSP 
10  READ ( MG ) DUMMY 
DO  27  I XL= 1 , NNUK 
C 7  READ  FOURTH  DATA  CARD 

READ) 5. 11) AID) I  XL) .DENT) I  XL) 
ll  FORMAT) 1F12.7*1E13. 6) 

DO  31  JZ= 1 .NNOT 
J2  =  JZ 

C8  THE  NEXT  CARD  CHECKS  TO  SEE  IF  THE  ID  NUMBERS  THAT  ARE  INPUT 

C  MATCH  AN  ID  NUMBER  ON  THE  DATA  TAPE 

I F ( ABS (AID) I  XL ) -DD ( 1,J2) ) -0. 0000 1 ) 27 , 2 7 , 3 1 
31  CONTINUE 

WRITE) 6.982 ) AID) I  XL ) , NBT , ( DD ( 1 . I ) . 1  =  1 .NNOT) 

9.82  FORMAT)  1H1//5X.  15H  NUCLIDE  NUMBER  F9.4  .31H  NOT  ON  FAST  DATA  TA 
XPE  NUMBER  I 6//5X . 27HNUCL IDES  ON  DATA  TAPE  ARE--// ( 30X.F9.4 ) ) 

27  CONTINUE 

READ ( MG )  ( NT  1 0 ( I ) . 1  =  1,1) 

C  SKIP  RESONANCE  DATA 

NMORE=2*NTID(l) 

DO  160  1=1 , NMORE 
160  READ ( MG ) DUMMY 

C  ZERO  OUT  CROSS  SECTION  ARRAYS 

DO  15  1=1,100 

SI-GQ)  I  )=0.0 
SIGA) I )=0.0 
ENG) I ) =0.0 
15  SIGT) I )=0.0 

DO  111  JSD= 1 , NNUK 
777  READ(MG) (ESS) I ) ,1=1,37) 

C9  THE  FOLLOWING  CARDS  THROUGH  CARD  NO.  66  ARE  USED  TO  FIND  THE 
C  DATA  FOR  THE  PROBLEM  ON  THE  DATA  TAPE 

IF) ABS (A  ID) JSD )-SS< 13) ) -0 . 00001 ) 2C » 20 » 21 
21  NRK=SS ( 30 ) +0 • 1 

CIO  SKIP  NRK  RECORDS  ON  THE  DATA  TAPE 
DO  666  1=1, NRK 
666  READ ( MG ) DUMMY 
GO  TO  777 

Cll  THE  NEXT  4  CARDS  ARE  FOR  WRITING  DATA  ON  THE  SCRATCH  TAPES 
2o  WRITE(MT) (ESS) I ) ,1=1,37) 

WRITE(MTT) (ESS) I ) ,1=1,37) 
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WRITE (MMT )  ( ESS  C 15,1=1,37) 

WRITE'MMMT  )  (  ESS  (  15,1  =  1,37) 

SSS=SS ( 13 ) 

IF(SS(?6) )39,39,85 
83  WR I TE ( 6 ,40 )  SSS 

^0  FORMAT ( 12H1NUCLIDE  NO.  F 10 . 7 , IX , 13HHAS  1-D  ARRAY 
NX=SS(26)+0. 1 

C  SSI 26! =  N 0 •  Oc  ONE-DIMENSIONAL  ARRAYS 

LNX=6+NEV 
NT  =  NX*LNX 

C  READ  1-D  CROSS  SECTIONS 

READ (MG  5  ( TTT  ( I  5 ,1  =  1, NT) 

02  THE  NEXT  tW0  CARDS  CALCULATES  THE  FINE  GROUP  MACROSCOPIC  CROSS 

C  SECTIONS 

DO  33  1  =  1, NEV 

33  S I GA ( I ) =S I GA ( I 5  +  TTT I 1+6 ) *DENT ( JSD ) 

I Sl=LNX+6 

IF (SSI  2  5 ) 5  8997,8997,8996 

03  SS (25)  INDICATES  WHETHER  A  FISSION  CROSS  SECTION  IS  INCLUDED  FOR 

C  THIS  NUCLIDE 

8c96  DO  8999  1=1,. NEV 

04  THE  NEXT  3  CARDS  CALCULATE  THE  FINE  GROUP  MACRO.  CROSS  SECTION 
IS=IS1+I 

8o99  SIGQ'I! =SIGQ! I 5+TTTI I S ) *DENT ( JSD ) 

I  Sl= I E+6 

05  the  NEXT  3  CARDS  ASSIGN  THE  FINE  GROUP  NU  VALUES  TO  ENG  (  I  ) 

DO  9010  1  =  1, NEV 

IS=IS1+I 

9c 10  ENG ( I ) =ENG ( I ) +TTT (IS) 

8o97  CONTINUE 

39  I F ( SS (16) )91*91*41 
91  WR I TE(6.52)SSS 

5  9  F OR M  A  T  (  1 2H  ON  U  C  L ! D  E  NO.  FlO  •  7  •  IX  * 37HD0ES  NOT  HAVE  P-0  ♦  P - 1  *  D-2  *  P-3  A 

XRRAY 
CO  TO  98 

4 1  WRITE (6.4?) S5S 

42  FORMAT  (  1 2  WON'JCL  I DE  NO.  FlO  •  7  >  IX  *  29HHAS  P-0  , P-1  *  P-2  ,  AND  P-3  ARRAY  ) 
LT  =  SS ( 1 6 ) +C •  1 

C  READ  AND  STORE  P-0  ARRAY 

READ (VC)  (TTT(  I  )  »  I  =  1  •  LT ) 

DO  400°  1  =  1 *  LT 

Cl?  THE  NEXT  CARD  CALCULATES  MACRO.  CROSS  SECTIONS 
n  Q  Q  TTT!  I  )-TTTU  )-v-DENT(  JSD) 

WRI TE ( MT ) SS ( 13 ) *SS( 17 ) • SS ( 18 ) *LT • ( TTT (  I  ) > I  =  1 *LT ) 

CI5  READ  AND  STORE  PI  ARRAY 
READ ( MG )  (TTT(  I  )  ,I  =  1,LT) 

DO  A444  1=1 >Lt 

4444  TTT ( I } =  TTT ( I } *DENT( JSD ) 

WRITE ( MMMT ) SS  t 13 )  *SS ( 17 )  *  5  S (  18)  *LT*(TTT(I  )  ,  I *1 *LT ) 

Cl7  READ  AND  STORE  D2  ARRAY 
READ ( MG )  ( TTT ( I  )  ,I  =  1,LT) 

DO  444 5  I  =  1 #  L T 
4/+45  TTT ( I } =TTT ( I  }  *DENT  (  JSD  ) 

V/R I  T E  (  MMM T  )  SS  (  1 3  )  t  SS  (  17  )  *SS  (  18  )  »  LT  *  (  TTT  (  I  )  1 1  =  1  f  LT  ) 

Clc  READ  AND  STORE  P3  ARRAY 
READ ( MG ) { TTT ( I ) » I=ltLT) 

DO  4446  I = 1 »  L T 
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4/f46  TTT  (  I  }  -TTT  (  I  5  XDEMT (  JSD ) 

WRITE  (  5  S5  (  12  )  »SS(  17  )  ♦  5S<  18  )  t  LT  t  (  TTT(  I  )  *  I  -  .1  *'.T  J 

C  READ  AMD  5 TORE  I N r l A  5*  T I C  a  R P  A Y 

98  I F ( S£ (  19  3  )  82  ,C2,73W 
32  WRI TE ( S  » 53 ) SSS 

5?  FORMAT ( 12HONUCLIDE  NO.  FlO • 7 » IX » 29HDOES  MOT  HAVE  INELASTIC  ARRAY  ) 
GC  TO  39 

12  WRITE (6, 43) SSS 

43  FORMAT ( lOHCNUCLIDE  NO.  FlO • 7  *  IX » 19HHAS  INELASTIC  ARRAY  ; 

LT=SS ( 193+C.l 
READ (MG)  (TTT ( I ! , I  =  1 » LT ) 

DO  4001  1=1, LT 

'♦CO  1  TTT  !  I  )  =  TTt  (  I  )  *DENT  (  JSD  ! 

WRITF(MTT) SS ( 13) >SS( 20 ‘SSi 21 ) »LT  » ( TTT ' I ) ,  :  =  l,LT ) 

C  READ  AND  STORE  N-2N  ARRAY 

89  I F ( SS ( 22 )  ) 92  » 9? . 8? 

92  WR I T E ( 6  *  5 4 ) S  9 S 

54  FORMAT(  12HQNUCL  IDE  NO.  c  10 •  7  *  IX  *  2 HDOES  NOT  HAVE  N-2N  ARRAY  ) 

GO  TO  74 

8  3  WRITE' 6  »  44 ) SSS 

H  FORMAT' 12H0NUCLIDE  NO.  FlO .  7 .  IX »  14HHAS  N-2N  ARRAY  ) 

LT=SS ( 22 )+0. 1 
READ (MG) ! TTT ( I ) *  I  =  1 • LT ) 

DO  4002  1  =  1 »  LT 

4^02  TTT( I ) =TTT( I )*DENT( JSD) 

WRITE (MMT)SS( 13)  »SS(23)  »SS(24)  ,LT,  <  TTT(  I )  .  I  =  1,'_T! 

C  READ  SIGMA  TOTAL 

74  READ  (MG)  (TTT'I) .I  =  1»NEV) 

0  CALCULATE  TCTAL  MACROSCOPIC  CROSS  SECTION 

DC'  45  1=1  »NEV 

43  SIGT( I) -SIGT ' I )-*-TTT( i ) *DENT ( JSD ) 

C  REA r  TOTAL  ISOTROPIC  SCATTER  ARRAY ( P-0+ 1 NEL+2+N-2N ) 

IF( SS( 27 ) )10T. 105*94 
icc  WRITE (6, 60) SSS 

6o  FORMAT ( 12|JCNUCLIDE  NO.  FlO . 7  *  IN • 29HDOES  NOT  HAVE  TOTAL  ISO. ARRAY) 
CALL  EX  T  T 
94  WRITE '6, 46) SSS 

46  FORMAT ( 12HQNVCL I DE  MO.  FlO. 7 » IX » 28HHAS  TOTAL  ISC.  SCATTER  '"RAY 
LT=SS ( 27 ) +0. 1 
READ ' MG )  ( TTT ( I )  ,I  =  1»LT  ) 

DO  4004  I  =  1  *  LT 

4004  TTT' !  )=TTT(  !  )-“-DE.NT(  JSD) 

WRITE(MSS)SS' 13) »SS<28)*SS<29),LT»<TTT(I),I=i,LT) 

C  READ  SIGMA  SCATTER  TOTAL  FOR  P-C  ARRAY 

IF(SS( 16) ) 111 ,111, 106 

l06  READ  (MG)  (  TTT  (  I  )  » I  =  l.»NEV ) 

C  READ  SIGMA  SCATTER  TOTAL  FOR  P-1  ARRAY 

READ (MG) ( TTT ( I ) ,1=1, NEV) 

111  CONTINUE 

CALL  REV)  ( MSS  ) 

C2l  SUBROUTINE  REW  REWINDS  A  TAPE  AND  SETS  TT(K,K.M  EOUAL  TO  0. 

DO  120  J J  =  1 » NNUK 

READ (MSS)SS( 13! ,SS( 28) , SS ( 29 ) , LT , ( TTT ( I ) , I = 1 , LT ) 

LDF=SS ( 28 )  +C.1 

C22  LDF  =  THE  NUMBER  OF  GROUPS  SCATTERED  FROM 
LD  =  SS (  29  )  +  0.1 

C23  LD=  NUMBER  OF  GROUPS  SCATTERED  TO 
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CALL  ONE  t 1  ) 

?'i  SUBROUTINE  ONE  SETS  UP  THE  FINE  GROUP  MACRO.  SCATTERING 
CROSS  SECTIONS  FOR  A  TWO  ARRAY  VARIABLE.  TT(K.KK) 

1 20  CONTINUE 

CALCULATE  FLUX  OR  READ  IN  FLUX 
IF!KKK- 1)125,126.126 
1 25  READ <5. 127) ( F  LUX ( I ) .1-1,99) 

1 27  FORMAT (6E 12. 6) 

GO  TO  315 

126  READ(5,127)SSSS 

?r  THE  FOLLOWING  CARDS  THROUGH  CARD  NO.  29  ARE  FOR  CALCULATING 
THE  TLUX 

FLUX ( 1 ) = ( 1 • / ( S IGT ( 1 ) -TT ( 1 » 1 ) ) )*SSSS(1) 

DO  29  LL-2,99 
SUM-C.C 
KKKK-LL-1 
DO  30  J=1,KKKK 
30  SUM  =  SUM-tFLUX!J)*TT<  J.LL) 

29  FLUX (LL)=(1./(SIGT(LL)— TT(LL.LL) ) ) * ( SSSS  (  LL ) +SUM ) 

315  MM= 1 

TT ( 100 , 100 ) =0. 0 
FLUX ( 100 ) =FLUX ( 99 ) 

C26  THE  NEXT  8  CARDS  ARE  FOR  PRINTING  OUT  THE  FLUX 
WRI TE ( 6 , 1002 ) ( BXCX ( I ) , 1=1 , 18 ) 

WRITE (6, 211) 

2 1 1  FORMAT ( 1H0/7H  GROUP, 12H  FLUX 

WRITE (6, 215) ( I .FLUX ( I ) .1=1.50) 

2 i 5  FORMAT! 116, 1PE20. 6) 

WRITE (6, 10 02) (BXCX! I), I =1,18) 

WRI TE ( 6 . 211 ) 

WRI TE ( 6 ,215  )  l I .FLUX! I) .1=51,99) 

C  CALCULATE  BROAD  GROUP  MACRO.  ABS.  AID  TOTAL  CROSS  SECTION 

C27  THE  FOLLOWING  CARDS  THROUGH  9019  CALCULATE  THE  BROAD-' GROUP  ABS. 
C  TOTAL.FISSION.AND  NU*FISSION  CROSS  SECTIONS 

DO  201  I =1 ,NB3G 
SUM=0.0 
BUM-O.C 
SUMM=0.0 
SUMMM=C  •  0 
TUM=  C  »  0 
I  I  I =LBGB ( I  ) 

DO  200  I  I = MM , I  I  I 

BUM=BUM  +  ENG ( I  I >*FLUX( I  I ) *S I GQ ( I  I ) 

SUMMM  =  SUMMM+S I GQ ( I  I  )*FLUX(  I  I ) 

SUM=SUM+S IGT ( I  I )*FLUX< I  I  ) 

SUMM  =  SUMM+FLUX ( I  I ) 

200  TUM=  TUM+S I GA (II)*FLUX(II) 

FTOT ( I ) =SUMMM/SUMM 

TOT! I )=SUM/SUMM 
ABBS! I ) =TUM/SUMM 
TRA! I ) =BUM/SUMM 

201  MM=LBGB ( I )+l 

WRITE (6, 1002 ) ( BXCX! I ) , 1=1 ,18 ) 

WRITEI6.312) 

3 1 2  FORMAT ( 48HOBROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECT  I ONS/ 1H0 ) 
WRITE(6,217) 

2 1 7  FORMAT ( 6HOGROUP , 20H  SIGMA  ABSORPT I  ON » 20H  SIGMA  TOTAL 
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XH  SIGMA  FISSION 

WRITE) 6.218) ( I  * ABBS ( I ) ,TOT< I ) ,FTOT< I ) » I  =  1,NBBG) 

2i8  FORMAT! 1 I4.1PE20.6, 1PE18.6, 1PE18.6) 

WRITE) 6, 1002 ) ( BXCX  ( I ). 1=1,18 ) 

WR I TE ( 6  *  9009 ) 

9009  FORMAT ( 6HOGROUP ,28H  NU*SIGMA  FISSION 

9ol9  FORMAT) 114, 1PE28. 6) 

WRITE (6, 90 19) { I , TRA ( I ) »I=1,NBBG> 

C  CALCULATE  TOTAL  BROAD  GROUP  ISO. TRANSFER  SCATTER  CROSS  SECTION 

CALL  CSAV(SISO) 

C28  THE  SUBROUTINE  CSAV  FLUX  WEIGHTS  FINE  GROUP  SCATTERING  TRANSFER 

C  CROSS  SECTIONS 

C  CALCULATE  P-0  MACRO.  CROSS  SECTIONS 

CALL  REW(MT) 

DO  2010  JM= 1 »NNUK 
READ ( MT ) (ESS) I ) ,1=1,37) 

IF(SS) 16) ) 2010 *2010, 3000 

3o00  READ  t  MT ) SS ( 1 3 )  »SS( 17) ,SS( 18) *LT , (TTT ( I ) , I  =  1 , LT ) 

LDF=SS (17)  +0.1 

LD=SS (18)  +0.1 

CALL  ONE ( 1 ) 

2 olO  CONTINUE 

CALL  CSAV(PO) 

C29  THE  NEXT  10  CARDS  ARE  FOR  CALCULATING  THE  BROAD  GROUP  INELLASTIC 

C  SCATTERING  TRANSFER  CROSS  SECTIONS 

CALL  REW(MTT) 

DO  3021  JM  =  1 » NNUK 

READ ( MTT ) (ESS (I), 1=1, 37) 

I F ( SS (19)) 3021 ,3021 ,3022 

3022  READ) MTT ) SS ( 13 ) ,SS< 20 ) »SS ( 21 ) ,LT , < TTT < I ) *  I =1 ,LT > 

LDF=SS( 20)  +0.1 

LD=SS (21)  +0.1 

CALL  ONE  < 1 ) 

3o21  CONTINUE 

CALL  CSAV ( X  I NELS ) 

C30  THE  NEXT  10  CARDS  ARE  FOR  CALCULATING  THE  BROAD  GROUP  N-2N 

C  SCATTERING  TRANSFER  CROSS  SECTIONS 

CALL  REW(MMT) 

DO  3027  JM= 1 , NNUK 

READ ( MMT ) (ESS) I ) ,1=1,37) 

IF (SSI  22) >3027,3027,3028 

3028  READ (MMT )SS( 13 ) ,SS<  23 ) *SS ( 24 ) ,LT » ( TTT < I ) *  I  =  1 , LT ) 

LDF=SS (23)  +0.1 

LD=SS ( 24 )  +0.1 

CALL  ONE (2) 

3027  CONTINUE 

CALL  CSAV ( XN2N ) 

C3 l  THE  NEXT  10  CARDS  ARE  FOR  CALCULATING  THE  PI  SCATTERING  TRANSFER 

C  CROSS  SECTIONS 

CALL  REW(MMMT) 

DO  7330  I  1  =  1 , NNUK 

READ ( MMMT ) ( ESS ( I ) , I = 1 , 37 ) 

I F  <  SS (16) >7330,7330.7331 

7331  READ (MMMT) SSI  13) ,SS(17) ,SS< 18) ,LT* (TTT( I ) ,I  =  1,LT) 

READ ( MMMT ) DUMMY 
READ (MMMT)  DUMMY 
LDF=SS ( 1 7 ) +0 . 1 
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LD=SS( 18 )+0. 1 
CALL  ONE ( 1 ) 

7330  CONTINUE 

CALL  CSAV(Pl) 

C32  THE  BROAD  GROUP  P2  SCATTERING  TRANSFER  CROSS  SECTION  CALC  FOLLOW 
CALL  REW(MMMT) 

DO  7332  1  1  =  1 » NNUK 

READ ( MMMT ) (ESS( I ), 1=1,37) 

I F ( SS ( 1 6 ) >7332.7332.7333 

7333  READ ( MMMT ) DUMMY 

READ ( MMMT )SS( 13) ,SS(17).SS(18)*LT.(TTT(I).I=1,LT) 

READ (MMMT) DUMMY 
LDF  =  SS ( 1 7 ) +0  •  1 
LD=SS( 18  )+0»  1 
CALL  ONE ( 1  ) 

7332  CONTINUE 

CALL  CSAV ( P2 ) 

C33  THE  P3  BROAD  GROUP  SCATTERING  TRANSFER  CROSS  SECT  CALC  FOLLOW 
CALL  REW(MMMT) 

DO  7334  11=1. NNUK 

READ (MMMT) (ESS( I ) .1=1.37) 

I F ( SS ( 1 6 ) 17334.7334.7335 
7335  READ (MMMT) DUMMY 
READ (MMMT) DUMMY 

READ  (  MMMT  )  SS  (  13)  .  SS  (  17  )  .SS  (  18  )  .  LT  .  (  TTT  (  I  )  .  I  =  1  ,  LT*) 

LDF=SS (17) +0 • 1 
LD=SS ( 1 8 ) +0 • 1 
CALL  ONE ( 1  ) 

7334  CONTINUE 
CALL  CSAV ( P3  ) 

WRITE (6, 1002) (BXCXt I ) .1=1.18) 

C34 '  THE  NEXT  33  CARDS  ARE  FOR  PRINTING  OUT  BROAD  GROUP  CROSS  SECTION 
WR I TE ( 6 , 312  ) 

WRITE(6,551  ) 

551  FORMAT ( IOOHOGROUP  P-0  INELASTIC  N-2N 

*  TOTAL  SCATTER 

LPC  =  0 

DO  4100  LL=1 »NBBG 
DO  3200  I  =LL  .NBBG 
LPC=LPC+1. 

IF(l,PC-2  5)  3200.555. 555 
555  WRITE(6»1002) (BXCX(N) .N=1.18) 

WR I TE ( 6 . 3 12  ) 

WRITE ( 6.551 ) 

LPC  =  0 

3200  WRITE) 6.38 )LL, I .P0( LL. I ) .XINELSILL* I ) .XN2N(LL, I ) »SISO(LL. I ) 

4i00  CONTINUE 

38  FORMAT(5HOFROM,1I3.3H  TO, 113.1  PE  15. 6,1  PE  15. 6,1 PE  15. 6.1 PE  15. 6) 

WRITE (6, 1002) (BXCXt I ) ,1=1,18) 

WR I TE ( 6 , 3 1 2  ) 

WRITE (6*5511  ) 

5511  FORMAT ( 54H0  P-1  P-2  P-3  ) 

LPC  =  0 

DO  1400  LL= 1 » NBBG 
DO  2300  I  =LL  .N8BG 
LPC  =  LPC  + 1 

IF (LPC-25) 2300,5555, 5555 
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5^55  WRITE! 6*1002) ( 9XCX ( N ) » N= 1 • 18 ) 

WR I TE ( 6  *  3 12  ) 

WRITE (6, 55 11 ) 

LPC=0 

2?00  WRITE (6, 3 888 )LL» I *P1 (LL* I ) »P2 (LL, I ) .P3( LL, I) 

1400  CONTINUE 

3e83  FORMAT  (  5  HO  FROM  *  1 1 3  *  3H  TO.  1 1  3 '» 1PE 15 . 6 , 1PE 1 5 . 6 . 1PE15 . 6  ) 

WRITE (6. 1002 ) (BXCX! I ) .1=1.18) 

WR I TE ( 6  *  1424  ) 

1424  FORMAT!  65H0GR0UP  AVERAGE  COS  (  THETA  )  DIFF.  COEFc.  SIC-MA  TRANS 

♦PORT  i 

C35  THE  NEXT  5  CARDS  ARE  FOR  CALCULATING  THE  TRANSPORT  CROSS  SECTION 

C  .  THE  DIFFUSION  COEFF.  AND  THE  AVG.  COSINE  OF  THE  SCATT.  ANGLE 

1425  FORMAT! 1 1 5 . 1PE17.4 . 1PE18 .4 . 1PE17.4 ) 

DC  1431  LL=1»NBBG 
SIGTR=TOT(LL)-Pl(LL»LL)/3. 

D=1./(3.»SIGTR) 

XMUBAR= (Pit  LL . LL ) /3. ) /SI SO l LL »LL ) 

1431  WRITE (6, 1425 ) LL .XMUBAR . D , S I GTR 
REWIND  MG 
STOP 
END 

Sle^TC  REWE 

SUBROUTINE  REW(N) 

C36  SUBROUTINE  REW  REWINDS  A  DATA  TAPE  AND  SETS  TT  (K.»KK.)*0.0 

COMMON  TTT ( 5050 ) .TT(IOO.IOO) .  LDF  .LD  »LT  .NBbO  .LtKitt  1 99  )  .FLUXi  100  ) 

REWIND  N 
DO  3043  K.=  1.100 
DO  3043  KK=1 .100 
3043  TT(K»KK) =0.0 
RETURN 
END 

SlBFTC  ONEE 

SUBROUTINE  ONE(JJJ) 

C3~  THIS  SUBROUTINE  SETS  UP  THE  MACRO.  SCATTERING  TRANSFER  CROSS 
C  SECTIONS  IN  TERMS  OF  TT(K.KJC) 

COMMON  TTT ( 5050 ) »TT< 100.100) »LDF  »LD . LT .NBBG  *  LBGB ( 99 ) .FLUX! 100 ) 

N  =  1 
KK=0 
NSN=LD 
BB“JJJ 

DO  122  K=1.LDF 
DO  123  L=N»LD 
K.K=I001 

123  TT(K.KK) =  BB*TTT(L)+TT!K»K.K) 

N=LD+1 

KJC=K 

NNN=X+NSN 

IF (NNN-100) 101.101. 124 

124  NDIF=NNN-100 
LD=NSN+LD-ND IF 
GO  TO  122 

lOl  LD=NSN+LD 
122  CONTINUE 
RETURN 
END 

SlBFTC  CSAVE 
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SUBROUTINE  CSAV(SIG) 

C38  THIS  SUB.  FLUX  WEIGHTS  THE  FINE  GROUP  SCATTERING  TRANSFR  CROSS  SEC 
COMMON  TTT  (5050)  tTT(  lOOt  100)  *  LDF  »LD  t  LT  *NBBG  *L6GB  (  99  )  t  FL(JX(  100) 
DIMENSION  S I G ( 22  1 22  ) 

N=  1 

NNNN= 1 

DO  41  LL= 1 tNBBG 
DO  32  I  =LL  »NBBG 
X  =  0.0 
FF=0.0 
KK=LBGB(LL) 

DO  25  J  =  N  9 KX 
F  =  0  •  0 

MMM=LBGB ( I ) 

DO  10  K=NNNN  tMMM 
IF(J.GT.X)  GO  TO  10 
F=FLUX( J)*TT ( J»K)+F 
10  CONTINUE 
FF=F+FF 

25  X=X+FLUX(J) 

NNNN=LBGB ( I ) +1 
S I G ( LL  » I  )=FF/X 
32  CONTINUE 

N=L6GB(LL)+1 
NNNN=LBGB ( LL ) +1 
41  CONTINUE 
RETURN 
END 
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APPENDIX  B 

ENDF/B  Source  Deck  of  New  Barnyard 

On  the  following  pages  is  listed  the  source  deck  that 
uses  the  two  ENDF/B  data  tapes.  It  is  written  in  the  For¬ 
tran  IV  language  for  use  on  the  IBM  7094  digital  computer. 

It  consists  of  a  main  program  and  4  separate  subroutines. 

Input  data  from  data  cards  is  "read  in"  in  the  main 
program.  The  cross  section  data  required  for  the  problem 
is  obtained  in  the  subroutine  named  FIB  which  reads  the 
data  from  the  data  tape.  Three  scratch  tapes  are  used  to 
store  the  PN  scattering  transfer  cross  sections  until  they 
are  needed  in  the  main  program.  The  flux  calculation,  and 
the  broad  group  calculations  of  the  transport  cross  section, 
diffusion  coefficient,  and  the  average  cosine  of  the  scat¬ 
tering  angle  are  performed  in  the  main  program. 

The  subroutines  named  REW,  ONE,  and  CSAV  perform  the 
same  calculations  as  the  subroutines  with  the  same  names 
in  the  GGC-4  source  deck.  These  subroutines  were  discussed 
briefly  in  Appendix  A. 

Extra  comment  cards  have  been  added  to  the  listing 
that  follows  so  that  it  will  be  easier  to  read.  A  glos¬ 
sary  of  computer  program  symbols  is  given  in  Appendix  C. 
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C  NEW  BARNYARD 

C  MULTIGROUP  NEUTRON  MACROSCOPIC  CROSS  SECTION  CODE 

C  THIS  SOURCE  DECK.  USES  THE  ENDF/B  CROSS  SECTION  DATA  TAPES 

DIMENSION  A (21)  . PCX  20*20) » PI (20 .20) ♦ P2 ( 20 . 20 ) »P3 ( 20. 20 ) * 

1  TOT l 22 ) * ABBS ( 22 ) »SSSS(99) .BXCX< 18) 

COMMON  DENT ( 29 ) .NNUK .TT ( 103 ♦ 100 ) *8(10300) • NBBG  *FLUX ( 100) » 
*LBGB ( 22 ) *  I D  »  N 
DATA  Cl / 4H  • •  •  •  / 

Cl  SETUP  SCRATCH  TAPE  NUMBERS  AND  READ  IN  DATA  FROM  DATA  CARDS 
NTCH1=3 
NTCH2=4 
NGP= 100 
NGP3=NGP+3 
N1  =  0 
I TSN=8 
M0DE=2 


CO 

COi 

1006 

1002 


9871 

2 


61 


C2 


62 

24 

C3 

1001 


67 

63 

14 

15 


C 

C4 

C5 

166 


READ (5. 1006) l BXCX ( I) ♦ I =1 » 18 ) 

CARDS  FROM  CARD  NUMBERS  1006  TO  34  ARE  FOR  READING  IN  DATA  AND 
FOR  STORING  THIS  DATA  ON  A  SCRATCH  TAPE 
FORMAT ( 18A4 ) 

FORMAT ( 1H1 * 18A4 ) 

READ (5*9871)  NBBG ,NNUK*KKK 
READ (5. 7) ( LBGB ( I ) *I=1*NBBG) 

FORMAT ( 2413 ) 

FORMAT (313) 

FORMAT ( 1X.A4.2I2.1I6.1E13.4) 

DO  34  JKX=1,NNUK 

READ (5.2)  MATNO.LORDER*  N»  NOR.  DENT ( JKX ) 

CONTINUE 

N1=N1+100 

X=MATNO 

L0R1=9 

IF(NOR.EQ.O)  GO  TO  1001 

SKIP  NOR  RECORDS  ON  THE  ENDF/B  DATA  TAPE 
DO  62  1=1 *NOR 

READ (N*24) DUMMY 
FORMAT ( 1A4) 

READ  AND  CHECK  DATA  FROM  ENDF/B  DATA  TAPE 
READ ( N  »  20 ) ( A ( I ) .1=1.21) 

WRITE(6»20) ( A  <  I)  .1  =  1.21) 

IF(A(1).NE.C1.AND.A(11).NE.X)  GO  TO  67 
GO  TO  14 
WRI TE( 6.63 ) 

GO  TO  33 

FORMAT ( 37H0  ERROR  IN  SKIPPING  DATA  - CHECK  NOR) 

L  =  0 


CONTINUE 

CALL  FI B ( NGP  »NGP3 ) 

SUBROUTINE  FIB  IS  FOR  READING  THE  ENDF/B  DATA  TAPE 
ID=N1+L 


IF(L.GT.LORDER)  GO  TO  166 

WRITE  CROSS  SECTION  DATA  FOR  THIS  PROBLEM  ON  I TSN 
WRI TE ( IT SN ) (B(K)  .K* 1.10300) 

WRITE  CROSS  SECTION  DATA  INFO.  THAT  WILL  BE  USED 
WRI TE ( 6 .23  )  MATNO.L. ID 
IF ( L.EQ. 8 )  GO  TO  34 
READIN.20) (A( I ) ,1=1.21) 


SCRATCH  TAPE 
FOR  THIS  PROBLEM 
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WRI TE ( 6  *20 ) (A( I ) .1  =  1,21) 

I F ( A ( 1 )  .EQ.  Cl)  GO  TO  16 
WRI TE( 6 ,64 ) 

64  FORMAT (27H0A11)  IS  NOT  EQUAL  TO  .... 

GO  TO  33 
16  CONTINUE 
L  =  L+1 

IF ( L.EQ.LOR1 )  GO  TO  300 
IF(MODE.NE.O)  GO  TO  15 
3q0  CONTINUE 
34  CONTINUE 

20  FORMAT ( 9A4 ♦ A 1 ,10A4,A3) 

21  FORMAT (5X*A4,4X»1HP,I1,3X,I5,2X*8H  RECORDS) 

22  FORMAT (1H1.26X.23H  ****  MATERIAL  NUMBER  *A4,6X.I3,8H  GROUPS* 

»  5X , 1HP » 1 1 , 6H  ****  /) 

23  FORMAT  (  1H  ,  28H  THE  ID  NUMBER  ^OR  MATERIAL  ,A4*4H  P,U»4H  IS*i6) 
CALL  REWIITSN) 

C  READ  P-0  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

DO  73  1  =  1  »NNUK. 

READ  (  I TSN )  ( B  (  K. )  *K=  1,10300) 

DO  75  K.=  1  *NGP 

DO  75  J= 1 »NGP3 
:<K=NGP3*(K-1)+J 

C6  CARD  75  SETS  UP  A  2  ARRAY  VARIABLE  FOR  THE  SCATTERING  TRANSFER 
C6c  CROSS  SECTIONS  AND  IT  CALCULATES  THf  MACROSCOPIC  CROSS  SECTION 
75  T T  (  J  , K)  =  6  l  K.K. ) *DENT (  I  )+TT  U*A) 

UY=LGRDER 
DO  74  JI-l.JY 

C7  S.<IP  JYN  RECORDS  ON  THE  BINARY  SCRATCH  TAPE 
74  READ ( I TSN ) DUMMY 
73  CONTINUE 

C  CALCULATE  FLUX  OR  READ  IN  FLUX 

IF  (K.K.K.-1)  125 ,126,126 
Co  READ  IN  FLUX  VALUES 

125  READ ( 5,127) ( FLUX ( I ) ,1=1,99) 

127  FORMAT (6E12. 6) 

GO  TO  315 

C 9  READ  IN* SOURCE  VALUES 

1 2  6  READ ( 5 , 127 ) SSSS 

Clo  THE  FOLLOWING  CARDS  THROUGH  CARD  NO.  29  ARE  FOR  CALCULATING  THt 
CloO  FLUX 

FLUX ( 1 ) - l 1./ ( T  T (2,1) -TT (4*1) ) ) *( SSSS ( 1 )+SUM) 

DO  29  'LL  =  2 ,99 
3UM=0.0 
£.<KR=LL-1 
DO  30  J  =  1 , KK.KR 
NYN=LL- J+4 

30  S  U M  =  S  U M + F  L  U  X ( J ) *TT { NYN.LL) 

29  FL'JX(LL)  =  (  !•  /  (  TT  ( 3, LL ) -TT  1 4 ,LL ) )  )  *  (  SSSS (  LL )  +SUK ) 

315  MM=i 

.FLUX ( 100)= FLUX (99) 

C 1 1  PRINT  OUT  THE  FLUX  SPECTRUM 

WRITE (6,1002) CBXCXt I ) ,1=1,16) 

WRI TE ( 6 ,211 ) 

2i 1  FORMAT ( 1H0/7; .  SROUF , 12H  FLUX 

WR.TE( 6,215  )  C I ,FLUX  t  I  )  » 1  =  1 ,50 ) 

2i5  FORMAT! I6.1PE20.6) 
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WRITE (6, 1002 )  tBXCXl 1 )  *  1=1,16) 

WRITE16.211) 

WRITE! 6,215) ( I ,FLUX( I ), 1=51 ,99) 

c  calculate  broad  group  macro,  ass.  and  total  CROSS  SECTION 
CI2  THE  FOLLOWING  CARDS  THROUGH  CARD  NO.  216  ARE  FOR  CALCULATING  THE 
Cl20  TOTAL,  ABSORPTION,  AND  NU«FISSION  CROSS  SECTIONS 
DO  201  1=1 »NBDG 
SUM=0.0 
BUM=0»0 
SUMM=0 • 0 
T  UM=  0. 0 
1 1  I =LBGB ( I ) 

DC  20 0  1 1 =  MM » 1 1 1 
SUM-SUM+  TT(3»II)*FLUX!II) 

BUM=BUM+TT (2, I  I )*FLUX( I  I ) 

SUMM=SUMM+FLUX(II) 

2oO  TUM  =  TUM+TT ( 1, I  I )*FLUX( I  I ) 

TOT ( I ) =SUM/SUMM 
ABBS ( I ) =TUM/SUMM 
SSSS( I ) =BUM/SUMM 
201  MM=LBGB ( I ) +1 

WRITE (6, 1002) ( BXCX  ( I ) ,1=1,18) 

WRITE!  6*»  312 ) 

3 12  FORMAT (48HOBROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECT  I ONS/ 1H0 ) 
WRITEI6.217) 

2 17  FORMAT ( 6H0GR0UP , 20H  SIGMA  ABSORPT I  ON , 2 OH  SIGMA  TOTAL  ,20 

*H  NU*SIGMA  FISSION  ) 

WRITE (6, 218) ( I , ABBS ( I ) ,TOT( I ) ,SSSS( I ) , I=1,NBBG) 

2 18  FORMAT! 114, 1PE20. 6, 1PE18.6, 1PE20.6 ) 

CALL  CSAV(PO) 

CI3  SUBROUTINE  CSAV  IS  FOR  FLUX  WEIGHTING  THE  SCATTERING  TRANSFER 
CI3O  CROSS  SECTIONS 

DO  94  LL= 1 »  NBBG 
DO  95  I =LL»NBBG 
P 1  ( LL  ,  I )=0.0 
P2!LL, I ) =0.0 
95  P3 ( LL  » I ) =0 • 0 
94  CONTINUE 

C  THE  NEXT  CARD  IS  A  CHECK  TO  SEE  IF  PI  CROSS  SECTIONS  ARE  WANTED 

IFt  LORDER.GE . 1 )  GO  TO  78 
GO  TO  86 

7e  CALL  REW(ITSN) 

14  SUBROUTINE  REW  IS  FOR  REWINDING  THE  BINARY  SCRATCH  TAPE  AND  FOR 
I4O  SETTING  T(KK,K)  EQUAL  TO  0. 

READ  P-1  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

OBTAIN  P-1  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 

FOR  THE  99  GROUPS 

JY=LORDER-l 

CALL  ONE! ITSN,1,JY) 

C15  SUBROUTINE  ONE  IS  FOR  READING  THE  DATA  FROM  THE  SCRATCH  TAPE  AND 

C150  AND  FOR  ASSIGNING  THE  CROSS  SECTION  DATA  TO  A  VARIABLE 
CI5OO  THAT  HAS  TWO  DIMENSIONS  AND  FOR  CALCULATING  MACRO.  CROSS  SECTS 
CALL  CSAV (PI) 

I F ( LORDER.GE . 2 )  GO  TO  83 
GO  TO  86 

83  CALL  REW ( ITSN) 

C  READ  P-2  DATA  FROM  THE  BINARY  SCRATCH  TAPE 
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C  OBTAIN  P-2  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 

C  FOR  THE  99  GROUPS 

JY=LORDER-2 
CALL  ONE ( ITSN.2. JY) 

CALL  CSAV(P2) 

IFILORDER.GE.3)  GO  TO  84 
GO  TO  36 

84  CALL  REW ( I TSN ) 

C  READ  P-3  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

C  OBTAIN  P-3  SCATTERING  TRANSFER  M/CROSCOPIC  CROSS  SECTIONS 

C  FOR  THE  99  GROUPS 

JY=LORDER-3 
CALL  ONE ( I TSN  *  3  1  JY ) 

CALL  CSAV ( P3 ) 

86  DO  7899  LL=1.NB3G 
DO  7870  I =LL  » NBBG 

Cl6  WRITE  OUT  THE  BROAD  GROUP  SCATTERING  TRANSFER  CROSS  SECTIONS 

Cl60  ON  A  SCRATCH  TAPE  FOR  TEMPORARY  STORAGE 

7870  WRITE! NT  CHI ) PO ( LL  » I )  ,P1 < LL » I ) »P2 < LL *  I ) .P3 < LL , I ) 

7399  CONTINUE 

I F ( LORDER.LT .4 )  GO  TO  500 

85  DO  51  LL=1»NBBG 
DO  52  I =LL  *  NBBG 
POfLL.  I  ) =0.0 

PI ( LL» I ) =0.0 
P2 ( LL  »  I  ) =0  •  0 

5 2  P3 ( LL  *  I ) =0 • 0 

5 1  CONTINUE 

CALL  REW ( ITSN) 

C  READ  P-4  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

C  OBTAIN  P-4  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 

C  TOR  THE  99  GROUPS 

JY=L0RDER-4 
CALL  ONE ( ITSN. 4. JY) 

CALL  CSAV ( PO ) 

IF ( LORDER .GE. 5 )  GO  TO  87 
GO  TO  50 

87  CALL  REW (ITSN) 

C  READ  P-5  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

C  OBTAIN  P-5  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 

C  FOR  THE  99  GROUPS 

JY=LORDER-5 
CALL  ONE ( ITSN.5.JY) 

CALL  CSAV ( PI  ) 

I F ( LORDER. GE . 6  )  GO  TO  88 
GO  TO  50 

88  CALL  REW ( ITSN) 

C  READ  P-6  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

C  OBTAIN  P-6  SCATTERING  TRANSFER  M/CROSCOPIC  CROSS  SECTIONS 

C  FOR  THE  99  GROUPS 

JY=LGRDER-6 
CALL  ONE ( ITSN.6.JY) 

CALL  CSAV ( P2 ) 

I F ( LORDER. GE . 7  )  GO  TO  89 
GO  TO  50 

89  CALL  REW (ITSN) 

C  READ  P-  7  DATA  FROM  THE  3 1 NARY  SCRATCH  TAPE 
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OBTAIN  P-7  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 
FOR  THE  99  GROUPS 
JY  =  LORDER  -7 
CALL  ONE ( I TSN  *  7  *  JY ) 

CALL  CSAV ( P3 ) 

So  DO  789  LL~ 1 »NBBG 
DO  787  I =LL  »NBBG 

17  WRITE  OUT  THE  BROAD  GROUP  SCATTERING  TRANSFER  CROSS  SECTIONS 
FOR  P4  THROUGH  P7 

787  WRITE (NTCH2 ) PO ( LL , I  )  ,P1 ( LL *  I ) » P2 ( LL , I ) , P3 ( LL , I ) 

789  CONTINUE 
SqO  REWIND  NTCH1 

IF(LORDER.LT.4)  GO  TO  SOI 
REWIND  NTCH2 

18  CARDS  FROM  CARD  NO.  501  TO  551  ARE  FOR  PRINTING  OUT  CROSS  SECTION 
VALUES 

Sol  WRITE! 6 ,1002) (BXCX( I ) » 1=1. 18) 

WRI TE ( 6  *  312 ) 

WRI TE ( 6  *  551 ) 


LPC  =  0 

DO  4100  LL  =  1 1 NBBG 
DO  3200  I  =  LL  *  NBBG 

READ (NTCH1)P0(LL*I)  ,  PI ( LL  ,  I ) . P2 ( LL *  ± ) * P3 l LL . I ) 

LPC  =  LPC+ 1 

IT  (LPC-25)3200*555*555 
5^5  WRITE (6. 1002) (BXCX(NZ) *NZ=1,18) 

WRITE (6,312) 

WRITE! 6, 551) 

LPC  =  0 

3200  WRI TE(6,38)LL, I ,P0( LL, I ) ,P1 (LL, I ) ,P2(LL,I ) ,P3(LL,I ) 

4i00  CONTINUE 

38  FORMAT (5HOFROM,1 13, 3H  T0,1I3,1PE15.6,1PE15.6,1PE15.6,1PE15.6 
3888  FORMAT ( 5HOFROM , 1 1  3 , 3H  TO , 1 1  3 , 1P4E15 . 6 )  , 

5511  FORMAT ( 90H0GROUP  P-4  P-5 

*  P-7 

551  FORMAT ( lOOHOGROUP  P-0  P-1 

*  P-3 
WRITE (6, 1424) 


P-6 

) 

P-2 

) 


1424  FORMAT ( 1H1 »65H  GROUP  AVERAGE  COS (THETA)  DIFF.  COEFF.  SIGMA  T 

TRANSPORT 

1425  FORMAT (1I5,1PE17.4,1PE18.4,1PE17.4) 

I9  BEGIN  CALCULATION  OF  THE  TRANSPORT  CROSS  SECTION,  DIFFUSION 
COEFFICIENT  AND  THEAVG.  COSINE  OF  THE  SCATTERING  ANGLE 
DO  1431  LL=  1 »NBBG 


S1GTR=T0T (LL)-Pl (LL.LL) /3. 

D=1./(3.*5IGTR) 

XKUBAR- (Pl(LL,LL)/3«)/P0(LL,LL) 

1431  WRITE (6, 1425 ) LL .XMUBAR  ,D , S I GT R 
IF(L0RDER.LT.4)  GO  TO  33 

C2o  READ  THE  SCRATCH  TAPE  AND  PRINT  OUT  THE  BROAD  GROUP  SCATTERING 
C  TRANSFER  CROSS  SECTIONS 

WRITE (6, 1002) ( BXCX ( I ) ,1=1,18) 

WRITE16.312) 

WRI TE ( 6  *  55 1 1 ) 

LPC-0 

DO  1400  LL= 1 »NBBG 
DO  2300  I  =LL  »NBBG 
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READ ( NT CH2 ) PO ( LL  1 1 ) *  PI ( LL • 1 ) t P2 ( LL 1 1 ) * P3 ( LL 1 1 ) 

LPC=LPC+1 

IF(LPC-25)2300t 5555 ,5555 
5555  WRITE(6t 1002) (BXCX(NZ) ,NZ=1,18) 

WRI TE ( 6  >312 ) 

WRI TE ( 6 » 5511 ) 

LPC-0 

2300  WRITE(6,3888)LL,I  ,PO(LL,I ) *  PI ( LL , I ) , P2 ( LL , I ) , P3 ( LL ♦ I ) 

I4OO  CONTINUE 

IF ( LORDER. EQ.8 )  GO  TO  90 
GO  TO  33 

9o  CALL  REW ( I TSN ) 

READ  P-8  DATA  FROM  THE  BINARY  SCRATCH  TAPE 

OBTAIN  P-8  SCATTERING  TRANSFER  MACROSCOPIC  CROSS  SECTIONS 
FOR  THE  99  GROUPS 
JY=LORDER-S 
CALL  ONE ( I TSN t  6  t JY ) 

CALL  CSAV(PO) 

WRITE (6, 1002) ( BXCX ( I ) ,1  =  1 >18) 

WR I TE ( 6 , 3 12  ) 

WR I TE ( 6 , 5595 ) 

5395  FORMAT ( 25HCGROUP  P-8 

LPC  =  0 

DO  5678  LL= 1 , NBBG 
DO  5679  I =LL  > NBBG 
LPC=LPC+ 1 

IF (LPC-25) 5679,4568,4566 
456C  WRI TE ( 6 , 1002 } (BXCXCNZ) ,NZ=1,16) 

WRI TE ( 6 , 312 ) 

WRI TE ( 6 , 5595  ) 

LPC  =  0 

5679  W-RITE(6,4569)LL,I  ,P0(LL,I) 

4^69  FORMAT (5 HO FROM, 113, 3H  TO , 1 1  3 , 1PE 15 . 6 ) 

5678  CONTINUE 
33  STOP 
END 

SlBFTC  F I  SB 

SUBROUTINE  F 1 5 ( NGP , NGP3 ) 

C22  THIS  SUBROUTINE  IS  FOR  READING  THE  ENDF/B  DATA  TAPE  AND  FOR 
C  ASSIGNING  THE  CROSS  SECTION  VALUES  TO  THE  VARIABLE  S ( K ) 

DIMENSION  M ( 6 ) >TEMP(6) , WORD (6,1) 

COMMON  DENT  (29)  »NNUK.  »TT  (  103,100)  ,B(  10300)  *  NBBG  *  FLUX  (  1 00  )  > 
*LBGB ( 22 )  , I D  » N 
5  FORMAT ( 6 ( I  2 , IX , E9 • 0 ) ) 

JT=NGP*NGP3 
J  =  0 

lG  READ (N, 5 ) ( M ( K ) , TEMP ( K ) ,K=1,6) 

DO  7  K= 1 , 6 
7  WORD l K , 1 ) -TEMP ( K ) 

DO  15  K  =  1 >  5 

I F ( WORD ( K , 1 )  .EG.  O.C  .AND.  WCRD(K+i,l)  .EQ.  0*0)  GO  TO  20 
NA=K+i 
GO  TO  15 
2q  NA=K 

GO  TO  21 
13  CONTINUE 
2 1  DO  35  K= 1 »NA 
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Ir(M(K)«oT»x)  GO  i  \j  j>0 
J  =  J+1 

5 ( J ; =  TEMP ( K ) 

I F ( J • GE • JT )  GO  TO  40 
GO  TO  35 

3q  I F  ( M  (  K ) • l 0 • 0 )  GO  TO  35 
M(K)=M(K)-1 
J  =  J  +  1 

3( J)=TEMP(K) 

IF  (  w. •  GE •  JT  )  GO  TO  40 
GO  TO  30 
3  5  COiM  »  I NU  E 

IF(J.LT.JT)  GO  TO  10 


4q  CONTINUE 


RETURN 

ENO 


SUBROUTINE  REW(NN) 

THIS  SUBROUTINE  IS  FOR  REWINDING 
T T ( NR » K )  EQUAL  TO  0. 

C  uMMoN  u*  »_N  i  (  2 9  )  >  N*\iG  K  ♦  T  T  (  x  C  3  >  ▲  C  0  ) 
*LBGB(22) 1 1 D » N 
REMIND  N N 


DC  3043  K  =  1 9 1  CO 
DO  3043  KR=1»1C3 

3 q43  TT  (  A  K  9  K  }  -  0  •  0 
RETURN 
*_  ND 


TAPES  AND  FOR  S«-TTINo 
B(iC^OC)  »i\BbG»Fi-UX(  100)  » 


SIbFTC  cnee 

SUBROUTINE  ONE ( ITSNtUN*UY) 

CE  4  THIS  SUBROUTINE  IS  FOR  READING  Tr.L  CROSS  SECT  I  U»\^^T/a  r  ROM 
C  BINARY  SCR/^Toi  TAF»_  AND  ASSIGNING  THE  DATrv  TO  A  TWO  m*\RmY 

C  VARIABLE t  TTCJtK) 

COMMON  DENT ( 29 ) tNNUK  >  TT ( 103  9 100 )  >B ( 1C30C ) *NBBG  >FLUX ( iOO J  f 
**  L  3  G  3  (  2  2  )  1 1 D  9  N 
DO  SO  JM  =  1 9  NNUK 
DO  79  I  =  1 » J N 
7  9  read; itsn) dummy 

READ  (  IT.SNS  (B(K)  fR^itlCSOO) 

DO  31  K  =  1 9  ICC 
DO  81  J= 1 9 103 
KK  =  1 0 3*  (  K-l  )  J 

3 1  TT ( J  9  K) =3 ( KK )  *DENT  ( JM) +TT  (  J  9 fC ) 

IF(JY.EQ.O)  GO  TO  60 
DO  82  J  =  1 9 J Y 
82  read; itsn) dummy 
80  CONT I NUE 

RETURN 
END 

SIbETC  CSAVE 

SUBROUTINE  CSAV(SIC) 

C25  THIS  SUBROUTINE  FLUX  WEIGHTS  THE  FINE  GROUP  CROSS  SECTIONS 
C  TO  OBTAIN  BROAD  GROUP- CROSS  SECTIONS 

DIMENSION  5 1 G ( 20  9 20 ) 

COMMON  DENT  (  29  )  9 NNUK.  9  TT  { 10  3  9 100  )  9  3  (  10300  )  9  NBBG  9 FLUX  (  100  )  9 
#LBGB ( 22 )  9 1 D  9  N 
NB  =  1 
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NNNN= 1 

DO  41  LL=1 »NBBG 
DO  32  I =  LL  »NBBG 
X  =  0.0 
FF=0.0 
KK=LBGB(LL) 

DO  25  J=NB  »KK 
F  =  0.0 

MMM=LBGB ( I ) 

DO  10  K  =  NNNN  *MMM 
IF(J.GT.K)  GO  TO  10 
NMY-K-J+4 

F  =  FLUX(J)*TT ( NNY  *  K ) +F 
10  CONTINUE 
FF=F+FF 

25  X=X+FLUX ( J ) 

NNNN=LBGB ( I )+l 
SIGILL* I ) =FF/X 
32  CONTINUE 

NB=LBGB ( LL ) +1 
NNNN=LBGB(LL)+1 
41  CONTINUE 
RETURN 
END 


75 


GNE/PHYS  69-8 


APPENDIX  C 

Glossary  of  Computer  Program  Symbols 
The  following  Fortran  IV  nomenclature  is  used  in  both 
source  decks  of  New  Barnyard.  The  program  themselves  are 
listed  in  Appendices  A  and  B. 


Symbol 

Meaning  or  Use 

A  (I) 

:  Alphameric  variable*  for  tempor¬ 

ary  storage  of  information  concern¬ 
ing  cross  section  data  for  a  nuclide 

ABBS  (I) 

:  Absorption  cross  section  for  the 

I-th  broad  group 

AID  (IXL) 

:  Nuclide  I.D.  number  for  the 

IXL-th  nuclide 

AT(I) 

:  Alphameric  variable  used  for 

data  tape  description 

B(K) 

:  Single  array  variable  used  for 

the  temporary  storage  of  cross 
section  data 

BUM 

:  Summing  Variable 

BXCX(I) 

:  Alphameric  variable  used  for 
problem  description 

Cl 

:  Alphameric  variable  used  to  check 

alphameric  data  from  the  data  tapes 

CSAV 

:  Subroutine  that  flux  weights 

fine  group  cross  sections 

*  Alphameric  variable  refers  to  a  variable  which  is 
"read  in"  under  an  A  format. 


76 


GNE/PHYS  69-8 


Appendix  C  (Contd') 

Symbol 

Meaning  or  use 

D  (I) 

:  Diffusion  coefficient  for  the 
I-th  broad  group 

DAD (I) 

:  Alphameric  variable  used  for 
nuclide  description 

DD (1 ,  IXL) 

:  Nuclide  I.D.  number  of  the 

IXL-th  nuclide 

DD (2 ,  IXL) 

:  Resonance  parameter  of  the 

IXL-th  nuclide 

DD (3 ,IXL) 

:  Resonance  parameter  of  the 

IXL-th  nuclide 

DENT  (IXL) 

:  Number  density  of  IXL-th 

nuclide 

DUMMY 

:  The  first  variable  of  a  data 

record 

ENG  (I) 

:  Convenience  variable  for  tempor 

ary  storage  of  the  i-th  fine  group 
boundary  and  also  the  i-th  fine 
group  average  number  of  neutrons 
emitted  per  fission  event 

ESS  (I) 

:  Alphameric  variable  used  for 

nuclide  description 

F,  FF 

:  Summing  variables 

FIB 

:  Subroutine  used  to  read  the 

ENDF/B  data  tapes 

FLUX  (I) 

:  Neutron  flux  for  the  I-th 

fine  group 

FNBT 

:  Data  tape  number  on  the  GGC-4 
data  tape 

FTOT(I) 

:  Fission  cross  section  for  the 

I-th  broad  group 
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Symbol 

FPL 

GTH(I) 

ID 

III 

IS,  IX,  I 
ISI 

ITSN 

J2 

JNVN ,  JNY , 
JY,  JT 

KK 

KKK 

KKKK 


Appendix  C  (Contd') 

Meaning  or  Use 

:  Convenience  variable  for  problem 

constant 

:  Energy  boundary  of  the  I-th  ther¬ 

mal  fine  group 

:  Convenience  variable  used  as  an 

identification  number  for  a  particu¬ 
lar  set  of  data 

:  Fine  group  number 

:  Summing  variables 

:  Convenience  variable  for  a  pro¬ 

blem  constant 

:  Logical  unit  number  of  scratch 

tape 

:  Convenience  variable  for  a  pro¬ 

blem  constant 

JNYY  :  Broad  group  numbers 

:  Convenience  variable  for  problem 

constant 

:  Convenience  variable  for  desig¬ 

nating  the  KK-th  dimension  of  a  di¬ 
mensioned  variable  and  also  used  for 
a  fine  group  number 

:  Indicator  for  performing  a  flux 

calculation  or  inputting  a  flux 
spectrum 

:  Convenience  variable  for  tempor¬ 

ary  storage  of  a  fine  group  number 

:  Convenience  variable  for  problem 

constant 
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Appendix  C  (Contd1) 

Symbol 

Meaning  or  Use 

LD 

:  Number  of  fine  groups  scattered 
to 

IDF 

:  Number  of  fine  groups  scattered 
from 

LEN(I) 

:  I-th  resonance  parameter 

LLBBG(I)  ,LBGB(I) 

:  Lower  broad  group  boundaries  of 
the  I-th  group 

LNX 

:  Length  of  I-D  cross  section  array 

LORDER 

:  Order  of  the  PN  scattering  trans¬ 

fer  cross  section 

LORI 

:  Convenience  variable  for  problem 

constant 

LPC,  LRINK 

:  Summing  variables 

LT 

:  Length  of  cross  section  array 

M(K) 

:  Convenience  variable  used  for 

temporary  storage  of  cross  section 
data 

MATNO 

:  Material  number  of  a  nuclide 

MODE 

:  Convenience  variable  for  problem 

constant 

MG 

:  Logical  unit  number  of  GGC-4 
data  tape 

MM 

:  Convenience  variable  for  tempor¬ 

ary  storage  of  a  fine  group  number 

MMMT,  MMT,  MSS,  MT, 


MTT 

:  Logical  unit  numbers  of  scratch 

tapes 
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Appendix  C  (ContdV) 

Symbol 

Meaning  or  Use 

N 

:  Logical  unit  number  of  ENDF/B 

data  tape 

NA 

:  Summing  variable 

NB 

:  Convenience  variable  used  for 

temporary  storage  of  a  fine  group 
number 

NBBG 

:  Number  of  broad  groups  for  a 
particular  problem 

NBT 

:  Data  tape  number  on  the  GGC-4 

data  tape 

NBTC 

:  Data  tape  number  on  the  GGC-4 
data  tape 

NEP 

:  Number  of  fine  group  energy 

boundaries 

NES 

:  Convenience  variable  for  pro¬ 

blem  constant 

NEV 

:  Number  of  fine  groups 

NGP,  NGP3 

:  Convenience  variable  for  pro¬ 

blem  constant 

NGT 

:  Number  of  thermal  energy  boun¬ 

daries 

NMORE 

:  Number  of  records  that  contain 

resonance  data 

NNOT 

:  Number  of  nuclides  for  which 

data  is  listed  on  the  GGC-4  data 
tape 

NNNN 

:  Convenience  variable  used  for 

temporary  storage  of  a  fine  group 
number 

80 


GNE/PHYS  69-8 


Appendix  C  (Contd') 

Symbol 

Meaning  or  Use 

NNUK 

:  Number  of  nuclides  for  a  parti¬ 

cular  problem 

NOR 

:  Number  of  data  records 

NPRl 

:  Convenience  variable  for  problem 

constant 

NRK 

:  Number  of  data  records  for  a 

nuclide 

NSP 

:  Number  of  fission  sources  on  the 

GGC-4  data  tape 

NT 

:  Convenience  variable  for  problem 

constant 

NTCH1,  NTCH2 

:  Logical  unit  number  of  scratch 

tapes 

NTID  (I) 

:  Convenience  variable  used  for 

temporary  storage  of  miscellaneous 
data  from  the  GGC-4  data  tape 

NX 

:  Number  of  I-D  cross  section 

arrays  for  a  nuclide 

ONE 

:  Subroutine  for  calculating  mac¬ 
roscopic  scattering  transfer  cross 
sections 

PO(LL,I)  Pi (LL, I)  , 

P2(LL,I),  P3(LL,I)  :  Convenience  variables  used  for 

temporary  storage  of  PN  scattering 


REW 

transfer  cross  section  for  scatter¬ 
ing  transfer  cross  section  for 
scattering  from  broad  group  LL  to 
broad  group  I 

:  Subroutine  used  for  rewinding 

scratch  tapes  and  setting  a  two 
dimensional  scattering  transfer 
cross  section  variable  to  zero 
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Symbol 
SIGA  (I) 

SIGT  (I) 

SIGQ  (I) 

SIGTR  (I) 

SISO(LL,I) 

SS(I) 

sss 

SSSS(I) 

SUM,  SUMM, 
TEMP  (K) 

TTT(I) 

TRA  (I) 


Appendix  C  (Contd') 

Meaning  or  Use 

:  Absorption  cross  section  for  the 
I-th  fine  group 

:  Total  cross  section  for  the  I-th 

fine  group 

:  Fission  cross  section  for  the 
I-th  fine  group 

:  Transport  cross  section  for  the 

I-th  broad  group 

:  Total  PO  scattering  transfer 

cross  section  for  scattering  from 
broad  group  LL  to  broad  group  I 

:  Convenience  variable  used  for 

temporary  storage  of  specific  data 
for  a  nuclide 

:  Nuclide  I.D.  number 

:  Convenience  variable  used  for 

temporary  storage  of  the  source 
spectrum  and  also  the  fission  cross 
section  times  Nu. 

SUMMM  :  Summing  variables 

:  Convenience  variable  used  for 

the  temporary  storage  of  cross  sec¬ 
tion  data 

:  Total  cross  section  for  the  I-th 

broad  group 

:  Convenience  variable  for  tempor¬ 

ary  storage  of  energies  and  lethar¬ 
gies 
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Appendix  C  (Contd1) 

Symbol 

Meaning  or  Use 

TT(K,KK) 

:  Scattering  transfer  cross  section 
for  scatter  from  fine  group  K  to  fine 
group  KK  (GGC-4  source  deck) 

TT(K,KK) 

:  Scattering  transfer  cross  section 

for  scatter  from  fine  group  J  to  fine 
group  KK  where  K  is  equal  to  KK+4-J 
(ENDF/B  deck) 

TT ( 3 , I ) 

:  Total  cross  section  for  the  I-th 
fine  group  (ENDF/B  deck) 

TT ( 2 , I ) 

:  Fission  cross  section  times  Nu  for 

the  I-th  fine  group  (ENDF/B  deck) 

TT ( 1 , I ) 

:  Absorption  cross  section  for  the 
I-th  fine  group  (ENDF/B  deck) 

TTT(I) 

:  Convenience  variable  used  for 

temporary  storage  of  cross  section 
data  and  resonance  data 

XINELAS(LL,I) 

:  Inelastic  scattering  transfer  for 
cross  section  for  scatter  from  broad 
group  LL  to  broad  group  I 

XMUBAR(I) 

:  Average  cosine  of  the  scattering 
angle  for  broad  group  I 

XN2N(LL,I) 

:  N-2N  scattering  transfer  cross 
section  for  scatter  from  broad  group 
LL  to  broad  group  I 

83 


GNE/PHYS  69-8 


APPENDIX  D 

The  99  Fine  Group  Structure 
On  the  following  pages  are  listed  the  99  fine  group 
boundaries  and  the  corresponding  99  group  lethargies  used 
in  both  versions  of  New  Barnyard.  The  99  group  structure 
is  calculated  as  follows: 

Let  E^,  the  fifth  energy  point,  be  equal  to  10  MeV 
and  let  it  be  the  reference  energy.  Therefore, 


E5  =  10'  ev 

u5  =  0 


(Dl) 


-(-.4) 


The  first  energy  point  is  taken  to  be  E^=E,_  e 
The  next  49  points  are  determined  with  a  uniform 
lethargy  mesh  U  =  0.1  i.e., 

U1  =  -.4 

U.  =  U.  +0.1  i  =  2,  3,  ...50 
l  l-l 

(D2) 

and  the  next  50  points  are  given  by 

U.  =  U.  +0.25  i  =  51,  52,  . . .100 
i  i-l 

(D3) 

Energies  associated  with  these  lethargies  are  given  by 

-U4 

E .  =  E,-  e  1. 
l  5 
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GRCLP 

_ 1_ 

2 


ENERGY  INTERVAL! E.V.  ) 
1.491825E  07  TO  1.349859E  07 
1.349859E  07  TO  1.221403E  07 


LETHARGY  INTERVAL 
-4.000CCE-01  TO  -3.00000E-01 
-3.00CCCE-01  TO  -2.00000E-01 


4 

1.105171E  07 

TO 

l.OOOOOOE 

07 

-l.OOOCCE 

-01 

TO 

0. 

v* 

5 

-  _ 1..QQO0OOE  U7  TO 

9jtQ463I4fc 

Qfe 

_  Q* 

TO 

l.OOOOOE 

-01 

6 

9.0A8374E  06 

TO 

8.1873C8E 

06 

l.OOOCCE 

-01 

TO 

_  A  9.  V  y,w  y  y  u 

2.COOOOE 

-01 

_ i 

8« 187  308 E  06 

TO 

7.408182E 

06 

2.000CCE 

-01 

TO 

3.00000E 

-01 

8 

7.4C8182E  06 

TO 

6*70320 OE 

06 

3.00CCCE 

-01 

TO 

4  •  CQOOQE' 

-01 

9 

5 

d 

9 

TO 

6.0653Q7E 

04 

4.000  CCE1 

-01 

TO 

__5.COOOOE' 

-01 

10 

6.065307E  06 

TO 

5.488116E 

06 

5.000CCE 

-31 

TO 

6  •  COOOOE' 

-01 

_ 11 

_ 5,M0116E  06 

TO 

4.965853E 

06 

_  4*  005 CCE 

-Gi  Tn 

_7*  OOOOOE' 

-01 

12 

4.965853E  06 

TO 

4.493290E 

06 

7.00CCCE 

-01 

TO 

8.  OOOOOE' 

-01 

_  11 

4 . 49  3290  E  06 

TO 

4  *0656976 

06 

-ai 

TO 

9*  OOOOOE’ 

-01 

14 

4* 06 5697 E  06 

TO 

3.6787S4E 

06 

9.000CCE-01 

TO 

l.OOOOOE 

00 

-  —15. 

3.678794E  06 

TO 

3.3287UE 

06 

l.OOOCCE 

00 

TO 

1 • 10000E 

00 

16 

3.328711E  06 

TO 

3.011942E 

06 

1. IOC CCE 

00 

TO 

1*2000 OE 

00 

_ _ 1I_ 

_ 3, Q 1194 2 E  06 

TO 

2.725318E 

06 

_  1.20CCCE 

00 

TO 

1.30000E 

00 

18 

2. 725318E  06 

TO 

2.465970E 

06 

1.30CCCE 

00 

TO 

1.40000E 

00 

_ 19 

2 • 465970E  06 

TO 

2.2313C2E 

06 

1.40CCCE 

00 

TO 

1.50000E 

00 

20 

2.231302E  06 

TO 

2 .018  96  5E 

06 

1.50CCCE 

00 

TO 

1.60000E 

00 

_ 21 

2.0189656  06 

TO 

1.826835E 

06 

1.600 CCE 

00 

TO 

1.70000E 

00 

22 

1 • 826835E  06 

TO 

1.652989E 

06 

1.70GCCE 

00 

TO 

1.800Q0E 

00 

_ 23 

_ _ 1.652989E  06 

TQ„ 

1*4954066 

Q6 

1,800  CCE 

00 

TO 

l#  90000E 

00 

24 

1.495686E  06 

TO 

1.353353E 

06 

1. 900 CCE 

00 

TO 

2.COOOOE 

00 

_ 25 

1.353353E  06 

TO 

1.224564E 

06 

2.00CCCE 

00 

TO 

2.10000E 

00 

26 

1 • 2  24364  E  06 

TO 

1.108032E 

06 

2.10  'CCE 

00 

TO 

2.20000E 

00 

.  __27 

1.108022E  06 

TO 

1.002569E 

06 

2.200CCE 

00 

TO 

2.30000E 

00 

28 

1. 002589E  06 

TO 

9.071796E 

05 

2,300  CCE 

00 

TO 

2.40000E 

00 

_ 29 

_ 9,0 717.9 6f  05 

TO 

8,20.85011 

05 

2.400CCE 

00 

TO 

2.50000E 

00 

30 

8.208501 E  05 

TO 

7.427359E 

05 

2.50CCCE 

00 

TO 

2.60000E 

00 

_ 1L. 

_  7.427359E  05 

TO 

6.720552E 

05 

2.600CCE 

00 

TO 

2.70000E 

00 

32 

6. 720552E  05 

TO 

6.081CC7E 

05 

2.700CCE 

00 

TO 

2.80000E 

00 

_ 23. 

6.081007.1  05 

TO 

5 .502323E 

05 

2.800CCE 

00 

TO 

2.90000E 

00 

34 

5. 502323E  05 

TO 

4.978708E 

05 

2 . 900  CCE 

00 

TO 

3. OOOOOE 

00 

_ 25. 

_ 4. 978708E  05 

TO 

4,5049116 

05 

_  3.000CCE 

00 

TO 

3  •  10000E 

00 

36 

4.504921E  05 

TO 

4  .07622 1 E 

05 

3. 100 CCE 

00 

TO 

3.20000E 

00 

_ 17 

4,0762216  05 

TO 

3.688317E 

05 _ 

3,200  CCE 

00 

TO 

3.30000E 

00 

38 

3. 688317E  05 

TO 

3.337327E 

05 

3.30CCCE 

00 

TO 

3.40000E 

00 

39 

3, 337327E  05 

TO 

3.019739E 

05 

3 .400  C  CE 

00 

TO 

3.50000E 

00 

40 

3.019739E  05 

TO 

2.732373E 

05 

3. 500 CCE 

00 

TO 

3.60000E 

00 

—41 

2.732273E  05 

TO 

2.472353L 

05 

3.600CCE 

00 

TO 

3.70000E 

00 

42 

2.47 2353E  05 

TO 

2.23707  8E 

05 

3.70CCCE 

00 

TO 

3.80000E 

00 

41 

2,2210181  05 

TO 

2.024192E 

05 

3.800CCE 

00 

TO 

3.90000E 

00 

44 

2.024192E  05 

TO 

1.831564E 

05 

3.90CCCE 

00 

TO 

4. OOOOOE 

00 

45 

1.831564E  05 

TO 

1.657268E 

05 

4.000CCE 

00 

TO 

4.10000E 

00 

46 

1 • 657268E  05 

TO 

1.499558E 

05 

4. IOC  CCE 

00 

TO 

4.20000E 

DO 

. . _  47 

1^4995511115. 

IQ  _ 

1 » 35685 6£ 

05  ....  . 

.4.200CCE 

00 

IQ_ 

4.30000E 

00 

48 

1.356856E  05 

TO 

1.227734E 

05 

4.300CCE 

00 

TO 

4.40000E 

00 

.49 

1.227734E  05 

TO 

1.110900E 

05 

4.400  CCE 

00 

TO 

4.50000E 

00 

50 

1 . 110900E  05 

TO 

8.651698E 

04 

4. 500  CCE 

00 

TO 

4.75000E 

00 
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GRCOP  ENERGY  INTER VAL ( E .V .  )  LETHARGY  INTERVAL 


51 

8.6516986 

04 

TO 

6.737949E 

04 

4.750CCE 

J)0 

TO 

5.0000QE  00 

52 

6  •  737949E 

04 

TO 

5.247520E 

04 

5.00CCCE 

00 

TO 

5.25000E 

00 

53 

5.247520E 

04 

TO 

4.086773E 

0<L 

5#25QCCE 

00 

TO 

5.5Q000E 

00 

54 

4.086773E 

04 

TO 

3 • 182782E 

04 

5.500CCE 

00 

TO 

5.75000E 

00 

55 

3. 182782E 

04 

TO 

2.478753E 

04 

5.750CCE 

J)  0 

TO 

6.00000E 

00 

56 

2.478753E 

04 

TO 

1.930455E 

04 

6.00C COE 

00 

TO 

6.25000E 

00 

57 

1 • 930455  E 

04 

TO 

1,503440E  04 

.  6.250CCE 

00 

TO 

6.50000E 

00 

58 

1 • 503440E 

04 

TO 

1 . 170880E 

04 

6.50CCCE 

00 

TO 

6.75000E 

00 

59 

1* 170880E 

04 

TO 

9.118021E 

01 

6.75CCCE 

00_ 

TO 

7.0QQOOE 

oo 

60 

9.118823E 

03 

TO 

7.101746E 

03 

7.00CCCE 

00 

TO 

7.25000E 

00 

61 _ 

7 » 1C1746E 

Q3 

TO 

1^2£M6£ 

Q1 

7.250CCE 

00 

TO 

7.50000E 

00 

62 

5. 530846E 

03 

TO 

4.307427E 

03 

7.50CCCE 

00 

TO 

7.75000E 

00 

63 

4.307427E 

03 

TO 

3.354627E 

03 

7.750CCE 

00 

TO 

8 .000006 

00 

64 

3.354627E 

03 

TO 

2.612587E 

03 

8.000CCE 

00 

TO 

8.25000E 

00 

65 

2.6I2587E 

03 

TO 

2.034664E 

03 

8.259CCE 

00 

TO 

e.50000E 

00 

66 

2.034684E 

03 

TO 

1 • 584614E 

03 

8.50CCCE 

00 

TO 

8.75000E 

00 

67 _ 

1.584614E 

.03 

TO 

1 «234QS8E 

03 

8.75CCCE 

00 

TO 

5.COOOOE 

90 

68 

1.23409 8E 

03 

TO 

9.611169E 

02 

9.00000E 

00 

TO 

9, 25000E 

00 

69 

9.611169E 

02 

TO 

7.485186E 

02 

9.250CCE 

00 

TO 

9.50000E 

00 

70 

7.485186E 

02 

TO 

5.829468E 

02 

9.50CCCE 

00 

TO 

9.75000E 

00 

71 

5.029468E 

02 

TO 

4.539955E 

02 

9.750CCE 

00 

TO 

l.OOOOOE 

01 

72 

4. 53 9995 E 

02 

TO 

3.535751E 

02 

l.OOOCCE 

01 

TO 

1.02500E 

01 

73 

3.535751E 

02 

TO 

2.753646E 

02 

1.025CCE 

01 

TO 

1.05000E 

01 

74 

2.753646E 

02 

TO 

2.144542E 

02 

1.050CCE 

01 

TO 

1.C7500E 

01 

75 

2 • 144542E 

02 

TO 

1 *67017  IE 

02 

1.075CCE 

01 

TO 

1.10000E 

01 

76 

1 . 67017 1 E 

02 

TO 

1.300730E 

02 

1*100  COE 

01 

TO 

1.12500E 

01 

77 

L.300730E 

02 

TO 

1.013010E 

02 

1.125CCE 

01 

TO 

1.15000E 

01 

78 

1.013010E 

02 

TO 

7.889328E 

01 

1.150CCE 

01 

TO 

1.17500E 

01 

19 _ 

7 «  889328E 

0L 

TO 

6.144214E 

01 

1.175CCE 

01 

TO 

1  *  20000E 

01 

80 

6.144214E 

01 

TO 

4 .785 1 19E 

01 

1.200CCE 

01 

Trf 

1.22500E 

01 

81 

4.785H9E 

01 

TO 

3.726654E 

01 

1.2250CE 

01 

TO 

1.25000E 

01 

82 

3. 726654E 

01 

TO 

2  .90232  IE 

01 

1.250CCE 

01 

TO 

1.27500E 

01 

83 

2. 902321E 

01 

TO 

2 . 260330E 

01 

1.275CCE 

01 

TO 

1.30000E 

01 

84 

2.260330E 

01 

TO 

1.760347E 

01 

1  • 300CCE 

01 

TO 

1.32500E 

01 

85 

1.760347E 

01 

TO 

1.370960E 

01 

1 .32  5  CCE 

01 

TO 

1.35000E 

01 

86 

1 • 370960E 

01 

TO 

1.0677C4E 

01 

1.35GCCE 

01 

TO 

1.37500E 

01 

87 

1.0677Q4E 

01 

TO 

8 .3152S0E 

00 

1.375CCE 

01 

TO 

1.40000E 

OL 

88 

8.31 5290  E 

00 

TO 

6.475955E 

00 

1.4000CE 

01 

TO 

1.42500E 

01 

89 

6.475955 E 

00 

TO 

5.043478E 

00 

1.425CCE 

01 

TO 

1.45000E 

01 

90 

5 • 04347  8  E 

00 

TO 

3 .92786  5E 

00 

1.450CCE 

01 

TO 

1.47500E 

01 

9 1 

3.927865E 

00 

TO 

3.059024E 

00 

1*475  CCE 

01 

TO 

1.50000E 

oi 

92 

3.059024E 

00 

TO 

2.382370E 

00 

1.500CCE 

01 

TO 

1.52500E 

01 

93 

2.382370E 

00 

TO 

1.8553S2E 

00 

1. 525CCE 

01 

TO 

1.55000E  01 

94 

1.855392E 

00 

TO 

1.444981E 

00 

1.550CCE 

01 

TO 

1.57500E 

01 

95 

1 • 444981 E 

00 

TO 

1.125352E 

00 

1.575CCE 

01 

TO 

1.60000E 

01 

96 

1.125352 E 

00 

TO 

8.764252E- 

01 

1.60GCCE 

01 

TO 

1.62500E 

01 

97 

8*  Z642jL2E" 

01 

IQ_ 

6.825607E- 

01 

1.625CCE, 

01 

XCL 

1.65000E  01 

98 

6.825607E- 

01 

TO 

5.315788E- 

01 

1.650CCE 

01 

TO 

1.67500E 

01 

99 

5.315788E- 

01 

TO 

4.139S40E- 

01 

1.675CCE 

01 

TO 

1 • 70000E„ 

01 
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APPENDIX  E 

Structure  of  the  Data  Tapes 

Tables  XIII  and  XIV  of  this  appendix  show  the  structure 
of  the  two  data  tapes  used  in  New  Barnyard.  Table  XV  lists 
explanations  for  some  parts  shown  in  Tables  XIII  and  XIV. 

If  one  desires  to  use  either  of  these  data  tapes  in  a  new 
computer  code,  he  should  study  these  tables  carefully. 


87 


GNE/PHYS  69-8 


TABLE  XIII 

Structure  of  the  GGC-4  Data  Tape 


Record 

Number  of 

Words 

Variable* 

Description 

i 

3 

NTID(l)  =  NBT 

Tape  identification  number. 

NTID(2)  =  NEP 

Number  of  fast  energy  bound¬ 
aries  (groups  +  1) . 

(Note  #1) 

NTID(3)  =  NGT 

Number  of  thermal  energy 
points.  (Note  #  2) 

2 

90 

AT  (1-90) 

Tape  description  (5  lines 
of  72  characters  each) . 

3 

1 

LAD  =  NNOT 

Number  of  nuclides  on  tape. 

Record  4  is 
records  4  = 

repeated  for 
LAD) 

each  nuclide  (number  of 

4 

21 

ESS  (1- 18)  = 

DAD  (1-18) 

Nuclide  description. 

ESS(19-21)= 
DAD  (19-21) 

Nuclide  I.D.  number,  number 
of  resolved  resonances, 
number  of  unresolved 

resonances . 

5 

4 

LEN(l) 

Resonance  data.  (Note  #3) 

LEN (2) 

Resonance  data.  (Note  #3) 

LEN(3) 

Resonance  data.  (Note  #3) 

LBS  =  LEN (4) 

Resonance  data.  (Note  #3) 

6 

LBS 
(max .  : 

=  5101) 

Resonance  data.  (Note  #3) 

7 

NELT  = 

NEP  + 

NEP  +  2 
1) 

TTT (1-NELT) 

Fast  energy  group  boundaries 
(NEP+1  values),  fast  leth¬ 
argy  boundaries  (NEP+1  val¬ 
ues),  fast  lethargy  inter¬ 
vals  (NEP  values) . 

8 

NGT 

TTT(l-NGT) 

Thermal  energy  points. 

(Note  #2) 

*  The  variable  names  that  appear  in  the  GGC-4  Source  deck  are  listed 
in  this  column.  This  tape  is  written  in  binary  and  no  Format  state¬ 
ments  are  required  when  'reading  in"  these  variables. 
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TableXH^Contd ' ) 


Record 

Number  of  Words 

Variable 

Description 

9 

i 

NSP 

Number  of  fission  spectra 

on  tape . 


Record  10  is  repeated  for  each  fission  spectrum  (num¬ 
ber  of  records  10  =  NSP) . 


10 

LNR=17+NEP 

ESS  (1-18) 

Description  of  fission 
source  spectrum. 

ESS  (19-LNR) 

Fission  source  spectrum 
for  each  fast  energy  group 
(Note  #4) 

11 

1 

MRESN 

Number  of  nuclides  with 
resonance  data. 

Records 

nuclide 

12  and  13  are  repeated  for  each  resonance 
(number  of  records  12  and  13  =  MRESN) . 

12 

3 

SS  (13) 

Nuclide  I.D.  number. 

(Note  #5) 

SS (14) 

Resonance  data. (Note  #3) 

SS (15) 

Resonance  data.  (Note  #3) 

13 

LT=6*SS (14)+ 
SS  (15)  +  9 

TT(l-LT) 

Resonance  data  for  nuclide 

Records 

14  through  25  are 

repeated  for  each  nuclide. 

14 

37 

ESS (1- 18) 

Nuclide  description . 

SS (13) 

Nuclide  I.D.  number 

SS(14) 

Resonance  data.  (Note  #3) 

SS (15) 

Resonance  data.  (Note  #3) 

SS(16) 

Length  of  PO  array  (same 
length  for  PI,  P2,  and  P3 
arrays)  .  (Note  7/5) 

SS (17) 

Number  of  groups  scattered 

from  PO  scattering. 
(Note  #6) 


SS (18) 
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Number  of  groups  scattered 
into  PO  scattering. 

(Note  #7) 
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Record 


TableXIIjCContd') 


Number  of  Words  Variable 


Description 


15 


16 


17 

18 

19 


SS(19-21) 

For  inelastic  array,  same 
as  SS  (16-18) . 

SS (22-24) 

For  n,  2n  array,  same  as 

SS  (16-18) .  (Note  #8) 

SS(25) 

Fission  index:  non-zero= 
fissionable  nuclide. 

SS (26)=NX 

Number  of  1-D  arrays 
(absorption,  fission,  etc.) 

SS(27-29) 

For  total  scatter  array, 
same  as  SS(16-18). 

SS (30) 

Number  of  records  for  this 
nuclide . 

SS (31) 

Mass  number  of  nuclides. 

If  number  of  1-D 
present . 

arrays  [ss(26)l 

is  zero,  record  15  not 

IMX=  (6+NEP)  NX 

TTT(l-IMX) 

One-dimensional  cross- 
section  arrays  in  form: 
type  number  (l=absorption, 
2=fission,  etc.),  descrip¬ 
tion  (5  words),  cross-sec¬ 
tion  for  each  energy  group. 
(Note  #9) 

Records  16  through  19  not  present 
SS(16)  is  zero. 

if  length  of  PO  array 

IDK=SS (16) 

TTT(l-IDK) 

PO  scattering  array  in  form 
<71-1,  <7l-*2,crl  —  3..., 
<^2-2,  cr2-3,..  .,  <73-3, 
<73-*  4,  ...  (Note  #10) 

IDK 

TT(l-IDK) 

PI  scattering  array  (same 
form  as  PO.  (Note  #11) 

IDK 

TT(l-IDK) 

P2  scattering  array  (same 
form  as  PO) .  (Note  #11) 

IDK 

TT(l-IDK) 

P3  scattering  array  (same 
form  as  PO) .  (Note  #11) 
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Record 

20 

21 

22 

23 

24 

25 

26 


Table  xilfontd') 

Number  of  Words  Variable  Description 

Record  20  not  present  if  length  of  inelastic  array  SS(19) 
is  zero. 

LIT— SS(19)  TTT(l-LIT)  Inelastic  scattering  array 

(same  form  as  P0).  (Note  #12) 

Record  21  not  present  if  length  of  n,  2n  array  jss(22)l  is 
zero, 

LIT=SS (22)  TTT(l-LIT)  n,  2n  scattering  array  (same 

form  as  P0) .  (Note  #12) 

NEV=Number  of  BST(l-NEV)  ertotal  for  each  group. 

groups  (Note  #13) 

Record  23  not  present  if  length  of  total  scatter  array 
(SS (27)]  is  zero. 

LAT— SS(27)  TTT(l-LAT)  Total  scattering  array  (same 

form  as  P0) .  (Note  #14  and  15) 

Records  24  and  25  not  present  if  P0  and  PI  arrays  omitted 
[SS(16)=0.J 

NEV  SCT(l-NEV)  ^Scatter  for  each  group. 

(Note  15) 

SPP(l-NEV)  C£l  scatter  for  each  group. 
(Note  17) 

-  End- of- file  mark. 
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TABLE  XIV 

Structure  of  the  ENDF/B  Data  Tapes 


Record  Number  of  Words 

Variable*  Description 

1  21 

A(I),I=1,21  Description  of  the  PO  data 

Format  (9A4,A1,  that  follows  (Note  #18) 
10A4,A3) 

2  12 

M(K) ,TEMP(K) ,  The  absorption,  fission, 

K=l,6  Format  total,  and  PO  scattering 

(6(I2,IX,E9 .0)  )  transfer  microscopic  cross 

sections  (Note  #19) 

Record  2  is  repeated  until  the  entire  PO  array  is 
established . 


3  21 

A(I),I=1,21  Description  of  the  PI  data 

Format  (9A4,A1,  that  follows. 

10A4,  A3) 

4  12 

M(K),TEMP(K)  The  PI  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 

(6(I2,(LX,E9.0))  sections  (Note  #20) 

Record  4  is  repeated  until  the  entire  PI  array  is  estab- 


lished . 

5  21 

A(I),I=1,21  Description  of  the  P2  data 

Format  (9A4,A1,  that  follows. 

10A4,A3) 

6  12 

M(K),TEMP(K)  The  P2  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 

(6(12,  1X,E9.0))  section  (Note  #11) 

Record  6  is  repeated  until  the  entire  P2  array  is 
established. 


7  21 

A(I),I=1,21  Description  of  the  P3  data 

Format  (9A4,A1,  that  follows 

10A4,  A3) 

8  12 

M(K) ,TEMP (K) ,  The  P3  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 

(6 (12, 1X,E9 . 0) )  sections. 

*  The  variable  names 

that  appear  in  the  ENDF/B  source  deck  are 

listed  in  this  column.  This  tape  is  written  in  BCD  and  re¬ 
quires  format  statements  when"reading  in"  these  variables. 
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Record  Number  of  Words 
9  21 


TABLE  XIV  (Contd  ' ) 

Variable 


Description 


10 


11 


12 


13 


14 


15 


16 


17 


12 


A(I),  1=1,21  Description  of  the  P4  data 
Format  (9A4,  Al,  that  follows 
10A4,  A3) 

M(K) ,TEMP (X) ,  The  P4  elastic  scattering 
K  =  1,  6  Format  transfer  microscopic  cross 
(6  (12  ,  IX,  section 
E9.0)) 


Record  10  is  repeated  until  the  entire  P4  array  is 
established . 

21  A(I),  1=1,  21  Description  of  the  P5  data 

Format  (9A4,  Al,  that  follows 
10A4,  A3) 

12  M(K) ,TEMP (K) ,  The  P5  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 
(6(12  ,1X,  E9.0))sections 

Record  12  is  repeated  until  the  entire  P5  array  is 
established . 


21 


12 


A(I),I=1,21  Description  of  the  P6  data 

Format  that  follows 

(9A4 , Al , 10A4 , A3) 

M(K) ,TEMP(K) ,  The  P6  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 

(6(12,  IX,  E9 .0) )sections 


Record  14  is  repeated  until  the  entire  P6  array  is 
established. 

21  A(I),I=1,21  Description  of  the  P7  data 

Format  that  follows 

(9A4,A1, 10A4,A3) 

12  M(K) ,  TEMP (K)  The  P7  elastic  scattering 

K=l,6  Format  transfer  microscopic  cross 

(6(12, 1X,E9.0) )  section 

Record  16  is  repeated  until  the  entire  P7  array  is 
established. 


21 


A(I),I=1,21  Description  of  the  P8  data 

Format  that  follows 

(9A4,A1, 10A4,A3) 
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Record 

18 


TABLE  XlVCContd*) 


Number  of  Words  Variable 


Description 


12 


M(K) ,TEMP (K) , 
K=l,6  Format 
(6 (12, IX, E9 . 0) ) 


The  P8  elastic  scattering 
transfer  microscopic  cross 
sections 


Record  18  is  repeated  until  the  entire  P8  array  is 
established. 


Records  1  through  18  are  repeated  for  each  nuclide. 
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TABLE  XV 


Comments  About  the  Data  Tapes 
Note  Number  Comments 

1  NEP  is  equal  to  one  more  than  the  number  of 
fine  groups.  This  number  is  100  for  a  99 
group  data  tape. 

2  This  information  is  for  a  thermal  spectrum 
calculation  and  is  not  used  in  the  GGC-4 
source  deck. 


3  The  resonance  data  included  on  this  tape  is 
not  used;  therefore,  no  effort  is  made  to 
describe  it. 

4  The  fission  spectrum  for  each  fine  group  is 
given  as  the  fractional  number  of  fission 
neutrons  bom  in  that  group,  .  Thus, 


1.000. 


5  SS (16)  is  the  total  length  of  the  P0  elastic 
transfer  array.  The  PI,  P2,  and  P3  arrays 
are  the  same  length  as  that  for  P0. 

6  SS (17)  is  the  number  of  incident  energy 
groups,  and  is  equal  to  99. 

7  SS(18)  is  the  maximum  number  of  secondary 
energy  groups  for  which  elastic  transfer 
cross  sections  are  given.  This  number  in¬ 
cludes  the  in-group  term. 

8  SS(19)  through  SS(24)  are  values  similar  to 
those  given  for  elastic  scattering.  Note: 
no  PI,  P2,  or  P3  arrays  are  given  for  in¬ 
elastic  or  (n,2n)  scattering. 
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Note  Number 


10 


TABLE  XV  (Contd') 

Comments 

The  I-D  arrays  have  been  given  a  numeric 
identification: 

I.D.  Number  Description  of  Reaction 

1  Absorption  (the  fission  cross 

section  plus  any  other  neutron 
removing  reaction) . 

2  Fission  cross  section. 

3  Nu,  \)  ,  the  average  number  of 

fission  neutrons  produced  per 
fission  event  (includes  delayed 
neutrons) . 


4 

5 

6 

7 

8 
9 

10 

11 


(n,^  )  radiative  capture 
(n,  p) 

(n,00 

(n,  d) 

(n,  t) 

(n,  n)p 
(n,  n)  d 
(n,  n)  t 


(The  1-0  arrays  are  given  in  sequential  order) 
The  elastic  P0  scattering  cross  sections  are 
given  as  a  continuous  array.  The  order  of 
giving  the  cross  sections  is  -^elas  .elas 

0  1-»1,0  l->2, 

-  elas  ^  elas  ^elas  .elas 

0 1— >3 ,  . . . ,  [  1+SS  (18)  -1]  ,  0  2— >2 , 0  2 — >3 , 

elas  .elas  .elas 

_ ,0  2— > [ 2+SS  (18)  -1]  ,  U  3  — >3 , 0  3— >4 , 

.elas 

o  3^[3+SS(18)-l]  , . 


96 


GNE/PHYS  69-8 


Note  Number 


TABLE  XV  (Contd1) 

Comments 


.elas  .  elas  ,  elas 

0  N->N,  0n->N+1,  On— >[N  +  SS(18)-1]  , 


11 


12 


13 


These  cross  sections  are  given  in  units  of 
barns.  There  are  SS(16)  total  terms  in  this 
array.  The  elastic  PI,  P2,  P3  arrays  are 
given  in  exactly  the  same  manner. 

Recall  PN  =  (2N+1)  *(XsN  (i— )  ;  therefore, 

PI  =  3*0sl(i— >j),  P2  =  5* 0 s2  (i — *j)  ,  etc. 
The  values  listed  on  the  tape  are  the  PN 
values  and  not  the  CTsn  (i — )  values. 

The  inelastic  and  (n,  2n)  scattering  arrays 
are  given  in  the  same  manner  as  the  PO  elas¬ 
tic  arrays.  However,  the  number  of  groups 
scattered  is  generally  taken  to  be  100  for 
these  reactions,  thus  for  the  inelastic 
scattering  array,  the  cross  sections  are 
given  as: 

inelas  inelas  inelas  ^inelas,  inelas 
0  Is— »1,  Ul — *2 ,  ...  0  1-^100,0  2^2, 0  2— >3, 
inelas  inelas  .inelas  inelas 

_ 0  2->100,O3— >3,  0  3->4,  _ ,0  3— >100, 

inelas  inelas 

_ 0  SS(20)^»SS(20)  ,  .  CTSS(20)— >100. 


For  a  99  group  library  tape  there  is  a  low 
energy  "dump"  group  (#100) .  This  group  re¬ 
presents  the  cross  sections  scattered  to  all 
energies  below  the  low  energy  boundary  of 
group  99  (0.414  ev) . 


BST(I),  Cf total  for  group  i  is  calculated  by 

‘  =[i  +  SS (18) -1]  , 

J.P0  (elastic) 

a... 

J  (i^>j) 


+2* 


i=100  n,2n 
Oj  Ci — 
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Note  Number 


14 


15 


16 


17 


18 


19 


TABLE  XV  (Contd') 

Comments 

SS (27)  is  the  total  size  of  the  total  trans¬ 
fer  array.  Since  the  data  tape  has  99  fine 
groups,  this  number  will  generally  be  5049. 

The  total  scattering  array  for  scattering 
from  i — >j  is  obtained  by 

total  scattering  P0  (elas)  inelas  n,2n 

Cf  =6  -Kf  +2*ff 

i— >j  i— *j  i— ►j  i— ►j 

Cf Scatter  for  group  i  is  calculated  by 


j  (i— >j) 


scatter  j=[ i+SS (18) -1]  P0 (elastic) 

ff,  ■  I  O 

j=i 

^P1  scatter  for  group  i  is  calculated  by 
PI  scatter  i= [ i+SS (18) -1] 


Cf 


Plfelastic) 


An  example  of  this  description  is 

. MATERIAL  NUMBER  1012  100 

GROUPS  P0 

These  99  group  cross  sections  are  given  as 

ABS  ,  fission  ^  total  .  nnn 

01  ,0\  99R 


^ABS  fission  total 


0  2  '1)0  2 


,  Cf, 9  ,  Cf2-*2,C/l-*2,98R 


'ABS  ...  fission  ^ total  .  . 

C fg  'y0  g  ,0  g  ,CTg-*g,tfg-l->g, 

Cf g-2-^g*  *  * 
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Note  Number 


20 


TABLE  XV  (Contd') 

Comments 


--ABS  ,  Fission 
0  99  >1/0 


99 


rf Total  rf 
*  U  99  *  ^99— >99, 


cC 


98— >99 


3R’  °*  ’  ^99— >100’  ^98 — >100 , . 


The  R  denotes  zero  and  the  number  in  front 
of  the  R  indicates  the  number  of  cross  sec 
tion  values  that  are  equal  to  zero. 

These  cross  sections  are  given  as 

3R,  3Cfl-^l,  99R 


3R,  3tfg->g,  3Cfg-l— *g,  3CTg-2->g,  ... 


3R»  3  ^99— >99*  3CT98_,9g,  ... 

3R,  0.  ,  3  C^99 — ^X00 ,  3(^98— >100’  *** 
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APPENDIX  F 

99  Group  Fission  Source  Spectrums 
99  group  fission  source  spectrums  are  presented  here 
for  the  following  nuclides: 

(1)  U-233 

(2)  U-235 

(3)  PU-239 

(4)  PU-241 

(5)  CF -252 

This  data  was  taken  from  pages  42-46  of  GA-4265,  "GAM-11, 

A  code  for  the  calculation  of  Fast-Neutron  Spectra  and 

Associated  Multi  Group  Constants"  by  G.  D.  Joanou  and 

J.  S.  Rudek.  Each  source  spectrum  has  been  normalized  to  1. 


100 


££EQ.,TRUM  .NUMBER _ 1  U-233  FISSION  SOURCE 
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APPENDIX  G 
Sample  Problems 

This  appendix  is  devoted  to  exhibiting  New  Barnyard 
output  for  selected  problems.  Further  explanation  of  some 
of  the  problems  is  given  below.  The  following  is  a  list  of 
the  sample  problems  included  in  this  appendix: 

(1)  Water:  5  group  set  using  GGC-4  source  deck 
(PO  to  P3  cross  sections  included) 

(2)  Magnesium:  4  group  set  using  ENDF/B  source 
deck.  (PO  to  P4  cross  sections  included.) 

For  both  of  these  problems  a  U-235  source  spectrum  was 
used.  Each  problem  has  a  printout  of  the  input  data  fol¬ 
lowed  by  the  cross  section  output.  The  output  cross  sec¬ 
tion  values  are  identified  and  they  need  no  explanation 
here.  The  printout  of  the  input  data  is  followed  by  the  cross 
section  output.  The  output  cross  section  values  are  identi¬ 
fied  and  they  need  no  explanation  here.  The  printout  of 
the  input  data  is  included  only  for  convenience  of  the  user. 
The  actual  source  decks  do  not  print  out  the  data  used. 
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THE  INPUT  DATA  FOR  THIS  PROBLEM  IS  ... 

WATER — 5  GKUUP  Sbl  USING  GGC-4  SOURCE  DECK - 

5  2  I 

10  22  4  I  lb  99 - - - - - 

I. 0000000  0.668000E-0I 

B.OzuuuUO  0 . 334000E-0 I  - - - 

0.412166E-04U.  113325E-U30.Z  I/033fc-030.609152E-03U7121752F^U20TZZ3Tl  1 E-02 
0.3  f9iU4L-UZU. 69 03 38E- 020.89303 7E-020. 125583E-010. 167906E-010.2I4541E-01 

0.263197E-OIO. 31 1313E-0 10. 356364E-0 10. 396 133E-0 10.42 8906E-0 10.453 567E -0 1 
O.  46962 Ot- 010.  4  T  7 1 2 5E-0 10 . 47660  IE-010. 4689O2E-75T0 .455093^-010 . 436340E-01 
0. 413815 fc-010. 3386 29E -010. 36178 7E-O10. 33415 9E—0l0. 306470E_0 10.279 30 lE-Ol 
O.Z53096E-010.228182E-010.204779E-010. 183021E-010. 162969E-010. l44&3YE-0r 
0. 12  7972E-010.  1 1Z9Z6E-010.994O8OE-O20.873164E-020.765445E-027T.  66982~3 E-02 
0. 585Z07E-020. 51054 l E-020. 444822E-020. 387106E-020. 336522E-020. 292267E-02 


0. 253612 E-020. 49 5 894 E— 020 .345783E~020.240338E-020.l66632E-020.il 5306E— 02 


0.796687E- 

■030. 549814E- 

-030. 379095E- 

-030. 

~07~ 

0. 

0. 

— rt - 

0. 

0. 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0. 

0. 

o . 

0. 

0* 

0. 

0. 

0. 

0 • 

0. 

0  • 

0. 

0  . 

0. 

o. 

0. 

0. 
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9WPHTS  WX 

WATER — 5  GROUP  SET  USING 

GGC-4  SOURCE 

DECK 

FAST  OATA  TAPE  NUMBER  = 

31568 

TAPE  DESCRIPTION . 

T  APE  I.D.  IS 

31568.0 

TH I S  T APE  I S  GGC4  VET 

304.0 

CONTENTS  OF  FAST  OAT A  TAPE  NUMBER  31568 


NUCLIDE  NU 

1.0000000 

1.2000000 

2. COOOCOO 
3.0062000 

3.0072000 

4.000UC00 

5.0000000 

5.0100000 

6.0200000 

7.0000000 

8. 02^0000 
nvoooocoo 

12.0000000 
13. COOOOOO 

14.0000000 

16.0000000 

2 0.0000000 
22".  0"0000  0  0 

24.0000000 
—  "25. 0000000 

26.0000000 
27. OOOOCOU 

28.0000000 

2^.0000000 

42.0000000 
48 • 0  COOOOO 
74.0000000 
74.  1  799994 
74. 1819992 
74.  1820993 

74.1829996 
74.1 830997 
74.1839991 

74. 1859999 
74.1 860991 

82.0000000 
~ 92.2  334995 
92.2349997 
9Z.  2 35 1999 
92.2379999 
97.2380991 

92.2381992 

94.2411995 


NUCLIDE  DESCRIPTION - 

HYDROGEN 

TEUTERIUM - 

HELIUM 

LITHIUM-6  ENDF/B  DATA AUGUST  1967 

LITHIUM-7  ENDF/B  DATA  AUGUST  1967 

BERYLLIUM - GA-5905 - 

BORON  NATURAL 
"BORON  rxr~ 

CARBON  ENDF/B  DATA  JULY  1967 

“NITROGEN 

OXYGEN  ENDF/B  OCTOBER  1967 

SODIUM 

MAGNESIUM 

ALUMINUM  GA-5884  - ^ 

SILICON 

SULFUR 

CALCIUM 

"TITANIUM - 

CHROMIUM 

MANGANESE - 

IRON 

CORALT 

NICKEL 

COPPER - - — 

MOLYBDENUM 

CADMIUM 

TUNGSTEN  GA-5885 

TUNGSTEN - 180 

TUNGSTEN  182  GA-5885 

TUNGSTEN  YW 2  RESONANCE  GA-5885 - 

TUNGSTEN  183  GA-5885 

TUNGSTEN  E83  RESONANCE  GA-5885 

TUNGSTEN  184  GA-5885 

TUNGSTEN  184  RESONAUUFTSA-5885 

TUNGSTEN  186  GA-5885 

TUNGSTEN  I"86  RESONANCE  GA-5885 - 

LEAD 

URANIUM  233  ENDF/B  JANUARY^ T  T968 

URANIUM  235  NASA  REPORT 

URANIUM  235  KAPLTENDF/B  DATA  FEB  1T67 

URANIUM  238  NASA  REPORT  JAN  1965 

URANIUM  238  RESONANCE  NASA  REPORT "JAN  1965 

URANIUM  238  ENDF/B  DATA  JULY  1967 

URANIUM  238  RESONANCE  ENDFYB  DATA  JULY  1967 

PLUTONIUM  241  ENDF/B  DATA  JANUARY  1967 


GNS/PHYS  69-8 

WATER— 5  GROUP  SET  USING  GGC-4  SOURCE  OECK 


FINE  GROUP  STRUCTURE 


GROUP 

ENERGY 

INTERVAL ( E • V. 

) 

LETHARGY  INTERVAL 

I 

1.491825E 

07 

TO 

1 • 349859E 

07 

-4.0Q000E 

-01 

TO 

-3.  OOOOOE' 

-01 

2 

1.349859E 

07 

TO 

1.221403E 

07 

-3. OOOOOE 

-01 

TO 

-2. OOOOOE 

-01 

3 

1.221403E 

07 

TO 

1. 105171E 

07 

-2. OOOOOE 

-01 

TO 

-l.OOOOOE 

-01 

4 

1.105171E 

07 

TO 

1 • OOOOOOE 

07 

-l.OOOOOE 

-01 

TO 

0. 

5 

1 • OOOOOOE 

07 

TO 

9 • 048374E 

06 

0. 

TO 

l.OOOOOE 

-01 

6 

9.048374E 

06 

TO 

8.187308E 

06 

l.OOOOOE 

-01 

TO 

2. OOOOOE 

-01 

7 

8.  187308E 

06 

TO 

7.408182E 

06 

2. OOOOOE 

-01 

TO 

3. OOOOOE 

-01 

8 

7 • 408 182E 

06 

TO 

6.703200E 

06 

3. OOOOOE 

-01 

TO 

4. OOOOOE 

-01 

9 

6.703200E 

06 

TO 

6 • 065307E 

06 

4.00000E 

-01 

TO 

5.00000E 

-01 

10 

6.065307E 

06 

TO 

5*488 1 1 6E 

06 

5. OOOOOE 

-01 

TO 

6. OOOOOE 

-01 

11 

5.488116E 

06 

TO 

4*  965853E 

06 

6. OOOOOE 

-01 

TO 

7. OOOOOE 

-01 

12 

4.965853E 

06 

TO 

4.493290E 

06 

7. OOOOOE 

-01 

TO 

8. OOOOOE 

-01 

13 

4.493290E 

06 

TO 

4. 065697E 

06 

8. OOOOOE 

-01 

TO 

9. OOOOOE 

-01 

14 

4.065697E 

06 

TO 

3.678794E 

06 

9. OOOOOE 

-01 

TO 

1*  OOOOOE 

Op 

15 

3. 678794E 

06 

TO 

3.328711E 

06 

l.OOOOOE 

00 

TO 

1 • 1 OOOOE 

00 

16 

3. 32871 IE 

06 

TO 

3.01 1942E 

06 

1.10000E 

00 

TO 

1.20000E 

00 

i? 

3.011942E 

06 

TO 

2 • 725318E 

06 

1.20000E 

00 

TO 

1.30000E 

00 

18 

2.725318E 

06 

TO 

2 • 465970E 

06 

1.30000E 

00 

TO 

1.40000E 

00 

19 

2 • 465970E 

06 

TO 

2.231302E 

06 

1.40000E 

00 

TO- 

1.50000E 

00 

20 

2.231302E 

06 

TO 

2.018965E 

06 

1.50000E 

00 

TO 

1.60000E 

00 

2.018965E 

06 

TO 

l • 826835E 

06 

1.60000E 

00 

TO 

1.70000E 

00 

22 

1  • 826835E 

06 

TO 

1 • 652989E 

06 

1.70000E 

00 

TO 

1.80000E 

00 

1.652989E 

06 

TO 

1 • 495686E 

06 

1.80000E 

00 

TO 

1.90000E 

00 

24 

1 . 495686E 

06 

TO 

1 • 353353E 

06 

1.90000E 

00 

TO 

2. OOOOOE 

00 

25 

1.353353E 

06 

TO 

1.224564E 

06 

2. OOOOOE 

00 

TO 

2.10000E 

00 

26 

1 • 224564E 

06 

TO 

1 • 108032E 

06 

2.10000E 

00 

TO 

2.20000E 

00 

27 

1. 108032E 

06 

TO 

1 • 002589E 

06 

2.20000E 

00 

TO 

2.30000E 

00 

28 

1* 002589E 

06 

TO 

9.071796E 

05 

2.30000E 

00 

TO 

2.4D000E 

00 

79 

9.071796E 

05 

TO 

8.208501E 

05 

2.40000E 

00 

TO 

2.50000E 

00 

30 

8 . 20850  IE 

05 

TO 

7.427359E 

05 

2.50000E 

00 

TO 

2.60000E 

00 

31 

7.  427359E 

05 

TO 

6.720552E 

05 

2.60000E 

00 

TO 

2.70000E 

00 

32 

6. 720552E 

05 

TO 

6.081007E 

05 

2. 70000E 

00 

TO 

2.80000E 

00 

33 

6.08 1007E 

05 

TO 

“5.502323E 

05 

2. 80000E 

00 

^ro 

2.90000E 

99^ 

34 

5. 502323E 

05 

TO 

4. 978708E 

05 

2.90000E 

00 

TO 

3. OOOOOE 

00 

35 

4. 978708E 

05 

TO 

4.504921E 

05 

3700000E 

00 

TO 

3. 10000E 

00 

36 

4. 504921E 

05 

TO 

4.076221E 

05 

3.1 OOOOE 

00 

TO 

3.20000E 

00 

37 

4. 076221E 

05 

TO 

3 • 688  3 1 7E 

05 

3.20000E 

09 

TO. 

3.30000E 

00 

38 

3.688317E 

05 

TO 

3 • 337327E 

05 

3.30000E 

00 

TO 

3.40000E 

00 

39 

3. 337327E 

05 

TO 

3.01973 9E 

05 

3. 40009E  00  TO  3.50000E 

00 

40 

3.019739E 

05 

TO 

2 • 732373E 

05 

3.50000E 

00 

TO 

3.60000E 

00 

41 

2. 732373E 

05 

TO 

2.472353E 

"05 

3.60000E 

00 

TO 

3.70000T 

00 

42 

2.472353E 

05 

TO 

2.237078E 

05 

3.70000E 

00 

TO 

3.80000E 

00 

43 

2. 237078E 

05 

TO 

2.024192E 

05 

3.B0000E 

00 

TO 

3.90000E  00 

44 

2.024192E 

05 

TO 

1 • 831564E 

05 

3.90000E 

00 

TO 

4. OOOOOE 

00 

45 

1 • 83 1564E 

05 

TO 

1 . 657268E 

"05 

4. O0OOOF 

00 

TO 

4.10000E 

00 

46 

1 • 657268E 

05 

TO 

1.499558E 

05 

4.1 OOOOE 

00 

TO 

4.20000E 

00 

4  7 

1 . 499558E 

05 

TO 

1.356856E 

"95 

4.20000E 

00 

TO 

4. 30000F 

99 

48 

1 . 356856E 

05 

TO 

1 . 227734E 

05 

4.30000E 

00 

TO 

4.40000E 

00 

49 

1. 227734E 

05 

TO 

l. 110900E 

U5 

4. 4000 OF 

00 

TO 

“4T500U0E 

00 

50 

1.  110900E 

05 

TO 

8.651698E 

04 

4.50000E 

00 

TO 

4.75000E 

00 
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GNR/PHYS  6£»8 

WATER — 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


FINE  GROUP  STRUCTURE 

GROUP 

ENERGY 

INTERVAL ( E . V . 

) 

LETHARGY  INTERVAL 

51 

8 • 651698E 

04  TO 

6.737949E 

04 

47T5000E 

00 

TO 

5 • OOOOOE 

00 

52 

6 • 7  37949  E 

04 

TO 

5 • 247520E 

04 

5.00000E 

00 

TO 

5.25000E 

00 

53 

5.Z47520E 

04 

TO 

4 • 086T73E 

04 

5.25000E 

00 

TO 

5 • 50000E 

90 

54 

4. 086773E 

04 

TO 

3.182782E 

04 

5.50000E 

00 

TO 

5.75000E 

00 

~ 55 

3.IR278ZE 

54 

TIT 

2.47875TF 

54 

5. 75000b 

00 

TO 

t>.  OOOOOE 

00 

56 

2.478753E 

04 

TO 

1 • 930455E 

04 

6.00000E 

00 

TO 

6.25000E 

00 

57 

1 • 930455E 

04 

TO 

1 • 503440E 

54 

6.55000E 

00 

TO 

6.50000E 

00 

58 

1  • 503440E 

04 

TO 

1 • 1 70880E 

04 

6.50000E 

00 

TO 

6.75000E 

00 

59 

1. 170880E 

04 

TO 

9. 118823E 

53 

67T5500E 

90 

TO 

7. OOOOOE 

00 

60 

9.  118823E 

03 

TO 

7. 101746E 

03 

7.00000E 

00 

TO 

7.25000E 

00 

61 

7.  10174  6IT 

0  3 

ITT 

5 • 530846E 

53 

7.25000E 

00 

TO 

7 . 5OO0OE 

90 

62 

5.530846E 

03 

TO 

4 • 307427E 

03 

7.50000E 

00 

TO 

7.75000E 

00 

63 

5  .397  42 7E 

03  TO 

3. 354F27E 

53 

7 .75000E 

00 

TO 

8. OOOOOE 

00 

64 

3. 354627E 

03 

TO 

2.612587E 

03 

8.00000E 

00 

TO 

8.25000E 

00 

55 

~2 • 61 2587E 

03 

TO 

Z.03?»5B4E 

“53 

BTZ5000E 

90 

TO 

8.50000E 

00 

66 

2 . 034684E 

03 

TO 

1.584614E 

03 

8.50000E 

00 

TO 

8.75000E 

00 

“67 

1. 584614"E" 

03 

1  U 

1 .234098E 

5T 

8775 DOGE 

00 

1 U 

9.0U000b  00 

68 

1 • 2  34098  E 

03 

TO 

9.611169E 

02 

9. OOOOOE 

00 

TO 

9.25000E 

00 

69 

9.611 1 69E  02 

TU 

7 • 435Y86E 

52 

972~5000E 

00 

TO 

9.50000E 

00 

70 

7.485186E 

02 

TO 

5.829468E 

02 

9.50000E 

00 

TO 

9.75000E 

00 

71 

5.  82 946 BE  7/2 

TO 

4.53  9995E 

52 

9775000E 

00 

TO 

I.OOOOOE 

01 

72 

4.539995E 

02 

TO 

3. 535751E 

02 

I.OOOOOE 

01 

TO 

1.02500E 

01 

73 

3.535751E 

02 

1  0 

2 • 75  3646E 

17  02  500E 

i . UyOOOt 

51 

74 

2. 753646E 

02 

TO 

2.144542E 

02 

1.05000E 

01 

TO 

1.07500E 

01 

75 

2.1 44542E  02 

75 

7 7  67017TE 

02 

1  • 07500E 

91 

TO 

1.10000E 

01 

76 

1 • 6701 7 1 E 

02 

TO 

1 • 300730E 

02 

1.10000E 

01 

TO 

1.12500E 

01 

77 

17  30  073 OE  02 

TO 

1 .9T39T9E 

02 

1  7T2500E 

51 

TO 

— I7T5000E 

01 

78 

1.013010E 

02 

TO 

7 • 889328E 

01 

1. 15000E 

01 

TO 

1.17500E 

01 

79 

7.889 32 8 E 

751 

TO 

6 • 1442 14E 

01 

1 • 1 7500E 

01 

TO 

1 .200005 

91 

80 

6. 1442 1 4E 

01 

TO 

4.785119E 

01 

1.20000E 

01 

TO 

1.22500E 

01 

~ 81T 

4 . 785 11 9E  01  TO 

3  .T 2  6  6  5  4  E 

01 

1.22500E  01 

TO 

1.T5000E 

01 

82 

3.726654E 

01 

TO 

2 • 90232 1 E 

01 

1  • 25000E 

01 

TO 

1.27500E 

01 

”53 

2  •  9  0232  IE 

5T 

TIT 

2.265335E 

Tir¬ 

1 . 27500E 

01 

TO 

1. 30000E 

01 

84 

2.260330E 

01 

TO 

1 • 760347E 

ol 

1.30000E 

01 

TO 

1.32500E 

01 

85 

t7560T4tf 

01 

TO 

1 .370960E 

01 

1T32505E- 

01 

T5 

T7T5000E 

91 

86 

1. 370960E 

01 

TO 

1.067704E 

01 

1.35000E 

01 

TO 

1.37500E 

01 

87  “ 

T. 56779 4E 

01  TO 

87 31 5290E 

00 

1 • 37500E 

51 

TO 

“T740000E 

01 

88 

8 •  3 15290  E 

00 

TO 

6.475955E 

00 

1.40000E 

01 

TO 

1.42500E 

01 

89 

6 • 475955E  00  TO 

5.  04 34 7 BE 

”00 

1  .4>2505E 

51 

TO r 

1 .45000 E 

01 

90 

5. 043478E 

00 

TO 

3.927865E 

00 

1.45000E 

01 

TO 

1.47500E 

01 

91 

3 . 927865E 

00 

ru 

3.059024F 

50 

1 . 47500E 

01 

TO 

I. 50000E 

51 

92 

3 • 059024  E 

00 

TO 

2. 382370E 

00 

1 . 50000E 

01 

TO 

1.52500E 

01 

93 

277823705 

00 

75“ 

1 .555T9ZE 

00 

1.52500E 

OT 

TO 

1 • 55000E 

01 

94 

1.855392E 

00 

TO 

1 .44498  IE 

00 

1.55000E 

01 

TO 

1.57500E 

01 

95 

1. 4 44 981 E 

50 

TO 

1  7175 J52E 

55 

1 • 57500E 

01 

TO 

1.60000E 

01 

96 

1. 125352E 

00 

TO 

8.764252E 

-01 

1.60000E 

01 

TO 

1.62500E 

01 

97 

8.T64252E- 

-01 

tct 

6 • 82  5607E 

-01 

1  • 62 500E 

01 

TO 

1 • 65000E 

OT 

98 

6. 825607E- 

-01 

TO 

5.315788E 

-01 

1.65000E 

01 

TO 

1.67500E 

01 

W  5. 31 578fTE-0 1  TO  4.I7W4GE-01  r."67500E  01  TO  1.70000E  01 
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ONS/PHYS  69-8 

WATER— 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


BROAD  GROUP  STRUCTURE 


BROAD  GROUP  ENERGY  INTERVAL! E.V. > 

1  1.491825E  07  TO  5.488116E  06 

2  5.488116E  06  TO  1.652989E  06 

3  1 . 652989E  06  TO  1.356856E  05 

4  I . 356856E  05  TO  1.670I71E  02 

- 5 - 1.670171E  02  TO  4.13994OE-01 


111  - 


gnVphys  69-8 

NUCLIOE  NO.  1 .0000000  HAS  1-0  ARRAY 

NUCLIDE  NO.  1.0000000  HAS  P-0 , P- 1 , P-2 , AND  P-3  ARRAY 
NUCLIDE  NO.  1.0000000  DOES  NOT  HAVE  INELASTIC  ARRAY 
NUCLIDE  NO.  1.0000000  DOES  NOT  HAVE  N-2N  ARRAY 


NUCLIDE  NO.  1.0000000  HAS  TOTAL  ISO.  SCATTER  ARRAY 


GW5/PHTS  69-8 

NUCLIOE  NO.  8.020000Q  HAS  1-0  ARRAY 


NUCLIDE  NO.  8.0200000  HAS  P-0, P- l , P-2 , AND  P-3  ARRAY 
NUCLIDE  NO.  8.0200000  HAS  INELASTIC  ARRA Y 


NUCLIDE  NO.  8 .0200000  DOES  NOT  HAVE  N-2N  ARRA Y 
NUCLIDE  NO.  8.0200000  HAS  TOTAL  ISO.  SCATTER  ARRAY 
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P  SET  USING  GGC-4  SOURCE  DECK 


GWJTJP 

tlux 

1 

5.232288E-04 

2 

1.462861E-03 

3 

3.336402E-03 

4 

7.471954E-03 

5 

1.472939E-02 

6 

2.565449E-02 

7 

4 • 33907 1 E-02 

8 

6.701 1 26E-02 

9 

1.016520E-01 

10 

1.316190E-01 

11 

1.797270E-01 

12 

2. 133987E-01 

13 

2.608718E-01 

14 

2 •  672685F-0 1 

15 

3.02031 7E-0 1 

16 

3 • 78652  5E-0 1 

17 

4 • 193974E-0 1 

18 

4.496019E-01 

19 

4.579706E-01 

20 

4. 147970E-01 

21 

3.91 1531E-01 

22 

4.168634E-01 

23 

4.063766E-01 

24 

4.044704E-01 

25 

3 .69869  1E-0 1 

26 

3 • 695762E-0 1 

27 

2 • 976449E-0 1 

28 

2.994022E-01 

29 

3.531688E-01 

30 

3.666530E-01 

31 

3. 153495E-01 

32 

2.894981E-01 

33 

2  •  707462E-01 

34 

2~.5DI992E-01 

35 

2 .040  345E-0 1 

36 

1.881376E-01 

37 

2.065525E-01 

38 

2.20I407E-01 

39 

2.07141 2E-0 1 

1.93  0225E- 0 1 

41 

1.81 1432E-0 1 

42 

1 .707197E-01 

43 

1.622364E-01 

44 

1 .55362  IE-01 

45 

1.495198E-01 

46 

1 . 4529T6E-0 1 

47 

1 . 36 14 1 8E-01 

4B 

1 . 304334E-0 1 

49 

1.252771E-01 

50 

2.931 348E-0 1 
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GHK/PHYS  69-8 

WATER — 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


OROUP - FTUX - 

51  2.709523E-01 

57  27T5T9B98E-01 

53  2.367955E-01 

54  27252021E-0I 

55  2.151739E-01 

56  2.083026E-UI 

57  2.03152  IE-01 

58  r.984380t-01 

59  1.947 144E-0 1 


KTl  1 - QO  T7  7  1 - - 

61  1 . 905074E-0 1 

"52  1.892594E-01 

63  1.886703E-01 

64  178  82078  E-OT 

65  1. 876260E-01 

66  I'.  869519E-U1 

67  1 . 862 143E-0 1 

88  1 . 8562U3E-0 1 

69  1 . 849905E-0 1 

TO  I7844139E-01 

71  1.839628E-01 


79 - 1  - 

•  l.OiOUC^t-Ul 

73  1 . 8  33264E-0 1 

74  r.831110E-01 

75  1 . 829322E-0 1 


7Z - 1 a  /m- — - - - 

'  °  I  •  O  C  I  7()7  t  “"U  i 

77  1 . 826802E-0 1 

"78  17824 064 E-OT 

79  1 . 820083E-01 

80  r.816798E-01 

81  1.814 1 70E-0 1 

"52  F78 12050E-0 1 

83  1.810321E-01 

84  1.808792E-01 

85  1 . 807562E-0 1 

88  1 . 8064 10E-01 

87  1.801672E-01 

88  1.7B9B32E-01 

89  1 . 779842E-0 1 

80  1.  /  / 1 898E-0 1 

91  1.765532E-01 

72  17  /6U2/5E-01 

93  1 . 755907E-01 


Q  A  i  r  n  i  i  n  i-  ^  , 

i .  » i^t-u  1 

95  1.748724E-01 


OTA - 1 - 1/  - frr - - - - - 

i  t-u  1 

97  1*7 19690E-0 1 


T77>84121E-0I 

1 .658085E-01 


W 
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WATER — 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 

BROAD 

GROUP  AVERAGED  MACROSCOPIC  CROSS  SECTIONS 

GROUP 

i - 

SIGMA  ABSORPTION 

SIGMA  TOTAL 

SIGMA  FISSION 

I 

2.  1 6  7554E-03 

1.250808E-01 

0. 

2 

3.959128E-04 

2 . 15531 IE-01 

0. 

3 

0  • 

5.21 8956E-0 I 

0. 

4 

6.7I5012E-05 

I.317152E  00 

0. 

5 

l . 69 1048E-03 

I . 508743E  00 

()• 

sms*! 


OUP  SET  USING  GGC-4  SOURCE  DECK 


GROUP 

- 1 - 

2 

- 3 - 

4 

- 5 - 


NU*SIGMA  FISSION 

- v~. - 

0. 

0". 

0. 

CTi - 
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WATER — 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


BROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECTIONS 


GROUP 

P-0 

INELASTIC 

N-2N 

TOTAL  SCATTER 

FROM 

1 

TO 

1 

3 .491491E— 02 

2.022698E-05 

0. 

3.493514E-02 

FROM 

1 

TO 

2 

6.29 I986E-02 

1.346142E-03 

0. 

6.426600E-02 

FROM 

1 

TO 

3 

2.01 1283E-02 

1.779472E-03 

0. 

2. 189231E-02 

FROM 

1 

TO 

4 

1 .796387E-03 

2 • 1 25040E-05 

0. 

1.817637E-03 

FROM 

1 

TO 

5 

2  •  208431E-06 

4.846245E-10 

0. 

2.2089I6E-06 

FROM 

2 

TO 

2 

1.013809E-0I 

0. 

0. 

1 . 0 I3809E-01 

FROM 

2 

TC 

3 

1  •050091E-01 

0. 

0. 

1.050091E-01 

FROM 

2 

TO 

4 

8.734505E-03 

0. 

0. _ 

8.734505E-03 

FROM 

2 

TO 

5 

1  •  073797E-05 

0. 

0. 

1.0  73797E-05 

FROM 

3 

TO 

3 

3 • 730391 E-01 

0. 

0. 

3 • 73039 1 E-0 1 

FROM 

3 

TO 

4 

1.486773E-01 

0. 

0. 

U486773E-01 

FROM 

3 

TC 

5 

1.788487E-04 

0. 

0. 

1 • 788487E-04 

FROM 

4 

TQ 

4 

U137157E  00 

0. 

0. 

1.137157E  00 

FROM 

4 

TO 

5 

1.794877E-01 

0. 

0. 

1.794877E-01 

FROM 

5 

TC 

5 

1 • 275828E  00 

0. 

0. 

I.275828E  00 

QHVPHYS  69-8 

WATER— 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


BROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECTIONS 


P-1 

P-2 

P-3 

FROM 

1 

TO 

1 

6.652408E-02 

8.839789E-02 

9.911263E-02 

FROM 

1 

TO 

2 

1.029138E-01 

8*  854647E-02 

-5.693567E-03 

FROM 

1 

TO 

3 

2. 168826E-02 

-2.922178E-02 

-5.571292E-02 

FROM 

1 

TO 

4 

5.223204E-04 

-4.348100E-03 

-I.789059E-03 

FROM 

1 

TO 

5 

2.255498E-08 

-5.520861E-06 

-7. 894038E-08 

FROM 

2 

TO 

2 

1.802595E-01 

2.221573E-01 

I.886510E-0I 

FROM 

2 

TO 

3 

1 • 669813E-01 

4.189492E-02 

-1.495191E-01 

FROM 

2 

TO 

4 

4.203339E-03 

-1.990882E-02 

-1.38 1 823E-02 

FROM 

2 

TO 

5 

1.815097E-07 

-2 • 68420 IE-05 

-6.352366E-07 

FROM 

3 

TO 

3 

6.2L1619E-01 

6.275355E-0I 

2 • 882293E-01 

FROM 

3 

TO 

4 

2.037582E-01 

-6.702652E-02 

-2.798133E-01 

FROM 

3 

TO 

5 

8.942934E-06 

-4.466725E-04 

-3. 12766 IE-05 

FROM 

4 

TO 

4 

2.166618E  00 

1 • 604800E  00 

3.325324E-01 

FROM 

4 

TO 

5 

2.356694E-01 

- 1 • 09872  5E-0 1 

- 3 • 32  3999E-0 1 

FROM 

5 

TO 

5 

2 •  474773E  00 

1 • 869701E  00 

4.282221E-01 

QNE/PHYS  69-8 

WATER— 5  GROUP  SET  USING  GGC-4  SOURCE  DECK 


GROUP 

AVERAGE  COS(THETA) 

DIFF.  COEFF. 

SIGMA  TRANSPORT 

T 

6. 3474F-0 I 

3.2392E  00 

1.029IE-01 

2 

5. 926RE-0 I 

2.L444E  00 

I  • 5544E-01 

3 

5.5505E-01 

1.0587E  00 

3  •  1484E-0 1 

4 

6.35 10  E-0 1 

5.6027E-01 

5.9495E-01 

5 

6.4658E-01 

4.8746E-01 

6.8382E-01 
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_  JA  FDR  THIS  PROBLEM  IS  ... - 

AVAGNESfUM — A  GROUP  SET  USING  ENDF/B  SOURCE  DECK 

4  I  I  - - 

13  27  54  71 

1014  4  2  12341  573900E-00 - - - 


O.41Z166E-U40.  ll33z:>E-03U.Zf  MJ33E-O30. 6091 52E-030. 12 1 752E-02P. Z23311E-02 

fl.3f9iO4E-uzo.6OO338b-uzo.893037E-020.125583E-010.  167906E-010.Z14541E-01 

c>.  Z6319/e-oiu.  31 1313E-U10. 356364E-U10.  3961 33b- 010.  428906E-010.463567E-01 
0.4696  20  E-010. 477125E-010.47660'1E- 01 0.468902  E-010.  455O93E-010.436340E-01 

O. 413815E-010. 388629E-010. 361 787E-010. 334159E-010. 3O6470E-010.279301E-01 

P. Z53O96b-oio.Z28iBZE-olo.zo4779E-010.183021E-010. 162969E-010. 144631E-01 

°*  127972E-P10. 112926E-010.994080E-020. 873164E-020. 765445E-020.669823E-02 
585Z07E-UZ0. 610541E-020.4448ZZE-020. 387106E-020. 336522 E-020.292267E-02 
^253612E-OZ0.495894E-OZ0.345783E-020.240338E-020.166632E-020.1153o6E-02 

796687E-030.  5498 14E- 030. 379095E- 030.  07 - ST - 


~0, 


0. 


0. 

07 


0. 

07 


0. 


0. 


0. 

07 


0. 

07 


0. 


0. 


0. 


0. 


0. 


0. 
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FOLLOWING  IMFO.  INDICATCS  WHETHER  OR  NOT  YOU  HAVE  THE  OATA  FOR  THE  RIGHT 

. . .  MATERIAL  NUMBER  1014  100  GROUPS  PO  moclioFS 

THE  ID  NUMBER  FOR  MATERIAL  1014  PO  IS  100 

. MATERIAL  NUMBER  1014 _ 100  GROUPS  P 1 

THE  ID  NUMBER  FOR  MATERIAL  1014  PI  IS  101 

. . MATERIAL  NUMBER  1014 _ 100  GROUPS _ P2 

THE  IU  NUMBER  FOR  MATERIAL  1014  P2  IS  102 

. MATERIAL  NUMBER  1014 _ 100  GROUPS  P3 _ 

THE  ID  NUMBER  FOR  MATERIAL  1014  P3  IS  103 

.  MATERIAL  NUMBER  1014  100  GROUPS  P4 

THE  ID  NUMBER  FOR  MATERIAL  1014  P4  IS  104 

. MATERIAL  NUMBER  1014  100  GROUPS  P5 

.  MATERIAL  NUMBER  1014  100  GROUPS  P6 

. «...  . _ MATERIAL  NUMBER  1014 _ 100  GROUPS _ P7 _ 

.  MATERIAL  NUMBER  1014  100  GROUPS  PB 
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MAGNESIUM-4 

GROUP  SET  USING  ENDF/B  SOURCE  DECK 

GROUP 

FLUX 

1 

9«  55  6081 E— 05 

2 

2.715339E-04 

3 

6.604603E-04 

4 

1.462057E-03 

5 

2.962775E-03 

6 

5.448673E-03 

7 

1.0C0503E-02 

8 

1.650977E-02 

9 

2. 829213E— 02 

10 

3.612760E-02 

11 

4.907892E-02 

12 

6.  52  3898E— 02 

13 

9* 182791 E— 02 

14 

1.420293E-01 

15 

1*  433325E-01 

16 

1*  86  5966E*-01 

17 

2.042672E-01 

18 

2.262544E-01 

19 

2*  943707E—01 

20 

2.939286E-01 

21 

3*  97  3576  E— 01 

22 

4.364193E-01 

23 

7*150182  E— 01 

24 

5.582233E-01 

25 

6.289372E-01 

26 

8.627642E-01 

27 

1* 048663E  00 

28 

1*  053914E  00 

29 

6.233043E-01 

30 

8*19 82 78 E— 01 

31 

7.346589E-01 

32 

9.804919E-01 

33 

8*60  3694E-01 

34 

8*  654983  E— 01 

35 

5.711158E-01 

36 

4.032375E-01 

37 

6*43  6866  E— 01 

38 

5* 168559E-01 

39 

3.627552E-01 

40 

2*66  8761 E-01 

41 

2.763459E-01 

42 

3.789079E-01 

43 

4 *  89  3963 E-01 

44 

5.830240E-01 

45 

6*  52  8828E-01 

46 

6*  92  52  56E-01 

47 

7.096412E-C1 

48 

6*88  5358E-01 

49 

6.048491E— 01 

50 

3.642916E-01 
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MAGNESIUM — 4 

GROUP  SET 

USING  ENDF/B  SOURCE  DECK 

GROUP 

FLUX 

51 

9.  72  6620EW01 

52 

1#  83 1392  E 

00 

53 

1.956T25E 

00 

54 

1.904309E 

00 

55 

1 • 77 1586E 

00 

56 

9.5*  6921E-01 

57 

1  •  70  43  83  E 

00 

58 

1.9481 84E 

00 

59 

1.917449E 

00 

60 

1.926907E 

00 

61 

1. 934948E 

00 

62 

1.941029E 

00 

63 

1.945448E 

00 

64 

1. 948694E 

00 

65 

1.950924E 

00 

66 

1. 952301 E 

00 

67 

1 . 95  2993  E 

00 

68 

1.953340E 

00 

69 

1.953515E 

00 

70  . 

1. 95  3345  E 

00 

71 

1.952657E 

00 

72 

1.951796E 

00 

73 

1 • 95  07  62  E 

00 

74 

1. 949555E 

00 

75 

1.948349E 

00 

76 

1.946971 E 

00 

77 

1. 945078  E 

00 

78 

1.943015E 

00 

79 

1  •  94  06 1 1 E 

00 

80 

1. 93  7868  E 

00 

81 

1.934787E 

00 

82 

1.93 1369E 

00 

83 

1 • 927447E 

00 

84 

1.923024E 

00 

85 

1. 91 8102  E 

00 

86 

1.912519E 

00 

87 

1.906110E 

00 

88 

1. 904762E 

00 

89 

1. 904762E 

00 

90 

1. 902049E 

00 

91 

1. 897636  E 

00 

92 

1  • 891695  E 

00 

93 

1.884401E 

00 

94 

1.878294E 

00 

95 

1  • 861902  E 

00 

96 

1. 842328E 

00 

97 

1. 820671 E 

00 

98 

1. 796367  E 

00 

99 

1.769217E 

00 

GHE/PHYS  69-8 

MAGNESIUM — 4  GROUP  SET  USING  ENDF/B  SOURCE  DECK 

-TJROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECTIONS 


GROUP 

SIGMA  ABSORPTION 

SIGMA  TOTAL 

NU*SIGMA  FISSION 

1 

1.389300E-02 

7.408909E-0I 

0. 

2 

3.125178E-03 

1.061269E  00 

0. 

3 

4. 178245E-03 

2.170218E  00 

0. 

—  4 

8.04I124E-04 

!• 405475E  00 

0. 
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MAGNESIUM — 4  GROUP  SET  USING  ENDF/B  SOURCE  DECK 
BROAD  lORO  UP  AVERAGED  MACROSCOPl  C  “CROSS  SECTIONS 


GROUP 

P-0 

P-T 

P-2 

P-3 

FROM 

i 

TO 

i 

4.150959E-OI 

6Y90566TE-01 

7.391705E-DT“ 

5*  947326E-01 

2 

2*912752  E-T5T 

-4.153767E-02  ““ 

2.015781E-02 

-3.5lBO52i-0Z 

FROM 

T 

TO 

3 

2.062810E-02 

0* 

0. 

“  0. 

FROM 

1 

TO 

4 

“0. 

0.  “ 

0* 

0* 

FROM  T 

TO 

2 

— 9.506294E-C1 

- 9. 948039E-01 

7.249279E-01 

2*1 604Z3E-01 

FROM 

~T 

TO 

3 

1.075208E-01 

-2.140035EH32 

-1. 247I54E-02 

-5 • 03311 81-02 

FROM 

~T 

TO 

4 

2.8903276-05 

0* 

0. 

0* 

FROM 

3 

TO 

3 

2 • 1 23 368E  00 

9.266635E-01 

9. 47551OF-0 1 

1 • 486644E-01 

FROM” 

3 

T0T 

4 

4.249124E-02 

— -4.851247E-02 

7.T76200E-0J 

-1.838741E-03 

~ From- 

~4 

TO- 

T" 

1 •377493E  00 

4.843637E-02 

“  8.T10037E-03" 

9.448551E-04 
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MAGNESIUM— 4  GROUP  SET  USING  ENDF/B  90URCE 
GROUP  AVERAGE  COS ( THETA )  OIFF.  COEFF. 

DECK 

SIGMA  TRANSPORT 

1 

5.5454E-01 

6.5270E-01 

5.1070E-01 

2 

3.4882E-01 

4.5683E-01 

7.2967E-01 

3 

I.4547E-01 

1.7908E-01 

1.8613E  00 

4 

1.1721E-02 

2.3992E-01 

1.3893E  00 
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“MAGNETS  I UH— 4  GtfOUP  SET  USING  ENDF/B  SOURCE  DECK 

BROAD  GROUP  AVERAGED  MACROSCOPIC  CROSS  SECTIONS 


GROUP 

P-4 

P-5 

P-6 

P-7 

FROM 

1 

TO 

1 

4.2387556-01 

0. 

0. 

0. 

FROM 

1 

TO 

2 

4.189399E-03 

0. 

Urn 

cr. 

FROM 

1 

TO 

3 

0. 

0. 

0. 

0. 

FROM 

1 

TO 

4 

0. 

07~ 

0. 

0. 

TROM 

2 

TO 

2 

6. 1326326-02 

0. 

0. 

0* 

FROM 

2 

TO 

3 

-1.0424306-02 

0. 

0. 

0. 

FROM 

2 

TO 

4 

0. 

0. 

0. 

0. 

FROM 

3 

TO 

3 

1.382124E-02 

0. 

0. 

0. 

FROM 

3 

TO 

4 

7.583821E-05 

0. 

0. 

0. 

FROM 

4 

TO 

4 

3 .794 49 16-04 

0. 

0. 

0. 
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